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SUMMARY 


This report describes the development of an aeroelastic analysis of a heli- 
copter rotor and its application to the simulation of helicopter vibration re- 
duction through higher harmonic control (HHC). An improved 
implementation of a finite-state, time-domain model of unsteady aerodynamics 
is developed to capture high frequency aerodynamic effects. Helicopter trim 
parameters are determined using an improved trim procedure accounting for 
flap, lag, and torsional deformations of the blade. A procedure for calculating 
4/rev rotor hub loads is developed and control algorithms are incorporated into 
the analysis to determine the optimal HHC inputs necessary to minimize these 
rotor hub loads. 

Using this analysis, the effect of unsteady aerodynamics on the aeroelastic 
response and stability of a hingeless rotor blade is studied in detail. The effect 
of time domain unsteady aerodynamics on blade aeroelastic stability and low 
frequency response is found to be small. On the other hand, the influence of 
unsteady aerodynamics on high frequency response, especially in the presence 
of high frequency pitch excitations, is found to be significant. Furthermore, 
the two aerodynamic formulations lead to significantly different response levels 
in rotor hub vibrations as a result of phasing changes introduced in open-loop 
HHC inputs. 

Several different HHC algoritms are implemented on a hingeless rotor using 
both quasisteady and unsteady aerodynamic models, and their effectiveness in 
reducing rotor hub vibratory shears is compared. A number of closed loop 
adaptive controllers are studied. All the controllers are found to be quite ef- 
fective in reducing vibrations, but very differing HHC inputs are required de- 
pending on the aerodynamic model used. Effects of HHC on rotor stability 
and power requirements are found to be quite small. Simulations of roughly 
equivalent articulated and hingeless rotors are carried out, and it is found that 
hingeless rotors can require considerably larger HHC inputs to reduce rotor 
vibratory shears. Additionally, hub moments for the hingeless rotor may in- 
crease significantly while shears are being minimized. This implies that the 
practical implementation of HHC on hingeless rotors might be considerably 
more difficult than its implementation on articulated rotors. Finally, this 
analysis is used to compare with some flight test results obtained on an OH6A 
light helicopter. 




Chapter I 

INTRODUCTION 


1.1 BACKGROUND 

A major goal of current helicopter research is to reduce the vibration levels 
experienced by crew, passengers or equipment during flight. This research has 
been driven by both commercial and military requirements. Commercial pas- 
senger acceptance would benefit greatly from the perception of the helicopter 
as having a “jet smooth ride.” Furthermore high vibration levels lead directly 
to higher maintenance costs. High vibration levels frequently cause problems 
during high speed flight and also during transition flight at 60-80 mph. De- 
creasing vibration levels and allowing higher cruise speeds would increase the 
load utilization of helicopters and so decrease relative capital costs. From a 
military point of view, increased speed leads to benefits in survivability and 
deployment response times. The same maintenance and comfort benefits as 
for civil operations apply. Reduction in vibration levels allows more accurate 
weapons deployment and more effective intelligence gathering, therefore cur- 
rent Army requirements specify vibration levels below 0.05 g. 

The traditional approach to vibration reduction in helicopters has depended 
on the use of passive means such as vibration absorbers or isolation devices. 
A comprehensive review of helicopter vibration control by Reichert[55] de- 
scribes many of these methods. The development of such methods has lead to 
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large decreases in vibration levels over the last thirty years and can now keep 
them below 0.10 g. However in the last ten years the trend has been to reduce 
acceptable vibration levels from 0.10 g to 0.05 g. This objective can be only 
rarely accomplished using a passive vibration control approach. 


1.2 REDUCTION OF HELICOPTER VIBRATIONS 

The desire to find methods to drastically reduce vibrations below the levels 
reached by these traditional means has lead to research in new areas. One 
promising approach is to design rotor blades which inherently have low levels 
of vibration. This may be done by applying optimum structural design to the 
aeroelastic tailoring of the blade. Geometry, mass, and stiffness distributions 
may be optimized to give minimum vibration levels at the rotor hub or some- 
where in the fuselage. The fuselage itself may also be tailored to reduce vi- 
brations at various points of interest such as the pilot seat, passenger 
compartment or the tail boom. Surveys of the application of structural opti- 
mization to helicopter vibration problems have been presented by 
Friedmann[20] and Miura[44]. Another approach is to optimize the use of 
conventional devices, such as vibration absorbers or isolators, via conformal 
mapping in order to reduce vibrations at critical points in the fuselage[l]. 

A different approach to vibration reduction, which is the subject matter of 
this research, is to use active controllers to reduce vibrations by eliminating 
them at their source, namely the aerodynamic excitation of the rotor. This 
concept, referred to as Higher Harmonic Control (HHC), relies on the appli- 
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cation of higher harmonic (i.e. above the 1/rev. pitch changes used for direc- 
tional control and trim) pitch changes to modify the blade airloads so as to 
minimize harmonic blade loading. For an N bladed rotor the predominant 
vibrations in the fuselage arc at N/rev. These are normally alleviated by ap- 
plying N/rev. pitch excitations superimposed on the collective (i.e. average), 
lateral (i.e. 1/rcv. sine), and longitudinal (i.e. 1/rev. cosine) pitch inputs used 
to control the helicopter attitude and velocity. This may be done by applying 
N/rev. harmonics to a standard helicopter swashplate through the use of hy- 
draulic servo-actuators. Aircraft flight tests[67, 68,43,53], wind tunnel 
tcsts[45,28,59, 37], and analytical simulations[62,46,10,34,41]. have shown 
that HHC is capable of substantial reduction in helicopter vibration levels en- 
countered in forward flight. 

The majority of all HHC studies, either analytical or experimental, have 
been based on linear, quasi-static, frequency domain representations of the 
helicopter response to control. Least-squares or Kalman filter type identifica- 
tion of the helicopter control parameters has been used along with a minimum 
variance or quadratic performance function type controller to determine opti- 
mal control harmonics for vibration alleviation. An extensive review of previ- 
ous work in this area is given in Ch. 6. The purpose of this research is not to 
present advances in the area of control theory, but rather to investigate its in- 
teraction with an aeroelastic model of the helicopter. 

Previous analytical studies have generally relied on simple analog[41] or 
frequency domain [10] models of the helicopter response. The helicopter ae- 
roelastic problem is inherently non-linear, with design parameters interacting 
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in a complex manner. Thus simplistic models of the helicopter response can 
be unreliable. Other studies[62,46] were based on using a fairly old aeroelastic 
response code, the G400 code[2,3], for simulation purposes. The aeroelastic 
model in the G400 code docs not have a consistent representation of the ge- 
ometrical nonlincarities due to moderate blade deflections. Other shortcom- 
ings of this simulation capability consist of the lack of time domain 
aerodynamics needed for capturing high frequency unsteady aerodynamic ef- 
fects combined with a step by step time integration solution technique which 
does not enable one to obtain direct information on blade stability, which is 
usually obtained by more recent methods such as numerical implementation 
of Floquet theory. 

Most previous studies involving HHC have been limited to conventional 
rotors consisting of articulated blades. In recent years the desire to decrease 
mechanical complexity and weight and minimize maintenance costs has lead 
to the development of hingeless and bearingless rotor hubs. In hingeless rotors 
the mechanical hinges in flap and lag present in articulated blades are replaced 
by a flexible cantilevered blade, where blade flexibility provides for virtual 
hinges. In such blades the mechanical pitch bearing is retained. Bearingless 
rotor blades are similar to hingeless blades except that the pitch bearing is 
eliminated and the pitch input is introduced through a torsionally flexible 
structural element. The mechanical simplicity and weight savings in hingeless 
and bearingless rotors is generally accompanied by increased vibration levels. 
The literature contains no aeroelastic simulation capability which is based on 
a hingeless or bearingless blade model and includes the effects of HHC. 
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1.3 OBJECTIVES OF THE RESEARCH 

This study has a number of objectives. The first group of objectives en- 
compass the modification of an existing acroclastic analysis, and associated 
computer code, to enable one to model HHC control effects with time domain 
unsteady aerodynamics. This is a major extension of an existing analysis and 
code because it requires the introduction of several new ingredients which are 
listed below: 

1. Adaptation of the aeroclastic analysis of Ref. 8 to allow modeling of 
rotor response to HHC. Coding must be provided for the application 
of HHC root pitch changes in the inertia and aerodynamic portions of 
the program. The Galcrkin finite element portion of the program must 
be modified to allow the modeling of articulated rotors. An improved 
procedure for calculating the N/'rev. harmonics of the six hub forces 
and moments to be used as objects of minimization must be developed. 

2. Development of an improvement on the unsteady aerodynamic theory 
of Ref. 13 overcoming its deficiencies. Algebraic manipulation soft- 
ware is used to develop implicit expressions for the aerodynamic loads 
which include previously neglected higher order terms including some 
apparent mass terms which can be important in calculating accurate 
helicopter response data. Methods for implementing these aerodyna- 
mics in an implicit manner on the aeroelastic model of this study must 
also be developed. 

3. Development of an improved trim procedure accounting for flap, lag, 
and torsional deformations of the blade and allowing modeling of real- 
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istic blade cross sectional property distributions. This procedure will 
provide for full equilibrium of the helicopter in forward flight. 

4. Provision of appropriate means to implement closed-loop HHC on the 
acroclastic model. This involves the coding of routines to calculate se- 
veral types of control and to implement filters for identification of 
problem parameters. 

These four objectives, which constitute the development of a significantly 
improved acroclastic model for the analytic study of the application of HHC 
for helicopter vibration reduction, are covered in detail in Chs. 4 through 6 of 
this research. 

With the acroclastic model developed, a second group of objectives was 
undertaken. These objectives, consisting of various comparative studies which 
are the ultimate aim of this research, are listed below: 

1. Determination of the importance of unsteady aerodynamic modeling on 
the response and stability of hingeless rotor blades in forward flight. 
Quasisteady aerodynamics, the unsteady aerodynamics of Ref. 13, and 
the improved unsteady aerodynamics of this study are evaluated. This 
will include an investigation of the importance of unsteady modeling 
when simulating application of HHC to rotors. In this simulation dif- 
ferences in high frequency response become important. 

2. Implementation of several HHC algorithms and identification ap- 
proaches in order to determine the effectiveness of different algorithms 
in reducing hub vibrations in steady flight and when step changes are 
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made in the flight condition. Differences in simulation results using 
quasisteady and unsteady aerodynamics will again be evaluated. 

3. Investigation of the effects of HHC on rotor stability. This has not been 
previously done because of a lack of appropriate simulation capacity. 
The solution procedure of this investigation is ideally suited to this since 
the Floquet theory used provides direct stability information in forward 
flight. 

4. Studying the effect of HHC on power requirements of the rotor. This 
is particularly relevant since small decreases in rotor torque or power 
requirement were observed in the first flight tests of an HHC 
systcm[67]. Such a decrease in rotor power requirements could make 
up for the power required for actuators to implement HHC. This small 
decrease has not been investigated and may in fact be an anomaly un- 
related to HHC. 

5. Simulation of the application of HHC to roughly equivalent articulated 
and hingeless rotors. Previous investigations have shown HHC to be 
very effective when applied to articulated rotors but comparisons with 
hingeless rotors have not been made. Hingeless rotors generally have 
somewhat higher vibration levels than articulated rotors and if larger 
HHC angles are needed to counter them this could decrease the attrac- 
tiveness of HHC for vibration reduction. 

6. Use of the aeroelastic model developed to simulate the application of 
HHC to an actual rotor, that of the OH6-A light helicopter. These re- 
sults will be compared with actual flight test data. Because of the limi- 
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tations of the model, particularly the lack of modeling of fuselage 
degrees of freedom, comparison will be limited to general trends in rotor 
response. 

Results of investigations toward these six objectives arc given in Chs. 7 
through 1 1 of this research. 


1.4 THE AEROELASTIC MODEL 

The research conducted in this study is carried out using an analytical mo- 
del which improves on some restrictions which were present in previous mod- 
els. A brief description of the aeroelastic model follows. A detailed description 
is presented in Ch. 4. 

The coupled flap-lag-torsional equations of motion which serve as the basis 
of this analysis are similar to those derived in Ref. 57 and are based on a 
structural model developed in Ref. 56. They contain geometrically nonlinear 
terms due to moderate blade deflections as illustrated in Figs. 2.2-2.4. These 
equations form the basis of an implicit flap-lag-torsional analysis[7,8] of a 
flexible, isotropic blade, modeled as an Euler- Bernoulli beam undergoing small 
strains and moderate deflections. Thus the equations contain geometrically 
nonlinear terms in the structural, inertia, and aerodynamic operators associ- 
ated with this aeroelastic problem. The spatial dependence in the equations is 
eliminated by using a Galerkin finite element analysis developed by Straub 
and Friedmann[61]. A modal coordinate transformation, using six rotating 
coupled modes, is performed to reduce the number of degrees of freedom. 
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These modes are calculated at a fixed value of collective pitch which depends 
only on advance ratio. For the configurations considered, the six lowest modes 
arc usually the first three flap, two lead-lag, and one torsional modes. The 
ordinary differential equations are solved using quasilincarization in an itera- 
tive manner to obtain the periodic equilibrium position in forward flight, for 
a propulsive trim type flight condition. 

The inertia loads are determined by using D'Alembert's principle. An im- 
plicit formulation for the aerodynamic loads is used. At each iteration an ap- 
proximation to the blade response is produced. This response is then used to 
generate numerical values of the modeling quantities needed in expressions for 
the aerodynamic loads to be used in the next iteration. Derivatives of the 
aerodynamic loads with respect to the generalized coordinates, required for 
stability analysis, are computed using finite difference approximations^]. The 
equations are linearized by writing perturbation equations about the nonlinear 
equilibrium position. Stability is determined by using Floquet theory. 

Use of this model for the research of this study requires major modifications 
and additions to the analysis of Ref. 8. The derivations required to make these 
changes are described in Chs. 3 through 6. 
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Chapter II 

GEOMETRIC COORDINATES AND TRANSFORMATIONS 

This chapter presents the coordinate systems and transformations used in 
deriving the equations of motion for the blade and in determining the ex- 
pressions for various forces and moments used in calculating the rotor hub 
loads. 

2.1 COORDINATE SYSTEMS 

2.1.1 Nonrotating , Hub-fixed Coordinate System 

This coordinate system, shown in Fig. 2.1, has its origin at the center of 
rotation of the rotor. The axes are x nr pointing toward the helicopter tail, y nr 
pointing to starboard, and z nr coinciding with the vector of rotor rotation. The 
corresponding unit vectors are i nr ,j nn and k nr . Hub shears and moments are 
defined in this coordinate system and are oriented along the coordinate axes 
as shown in Fig. 2.1. 


10 


2.1.2 Rotating, Hub-fixed Coordinate System 


This coordinate system, shown in Fig. 2.2, has its origin at the center of 
rotation of the rotor. The axes are x and y, which rotate in the x nn y nr plane, 
and z which is coincident with the z nr axis. The corresponding unit vectors arc 

A A A A 

i,j, and k. The coordinate system rotates at constant angular velocity Vtk 
around the z axis. 


2. 1 .3 Preconed, Undeformed Blade Coordinate System 

This coordinate system, shown in Fig. 2.2, has its origin at the blade root, 
offset from the center of rotation by a distance e , . The axes are x 0 and z 0 , in 
the x, z plane, and y 0 which is coincident with the y axis. The unit vectors e x , 

A A AAA 

e y , and e 2 are obtained by rotating the i ,j,k system by a precone angle p p 
about the y axis. This is the global structural coordinate system used for the 
finite element model. The x 0 axis represents the undeformed elastic axis of the 
blade, where the elastic axis is defined as the line connecting the shear centers 
of the blade cross-sections. A shear force applied at a point on the elastic axis 
will not cause torsional deformations. 
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2.1.4 Deformed Blade Coordinate System 

This coordinate system, shown in Figs. 2.3 and 2.4, represents the orien- 
tation of the local blade geometry after deformation. The unit vectors e x , e y , 
and e x are obtained by rotating the e x , e y , e 2 system around all three axes in 
accordance with the deformation of the blade. The e x unit vector is always 
parallel to the local deformed elastic axis. 


2.1.5 Rotated , Deformed Blade Coordinate System 

This coordinate system is the deformed blade coordinate system with the 
torsional deformation removed. The unit vectors % , fr and £ are obtained by 
rotating the e ' x , e y , % system by an angle of — 4>e’ x . This coordinate system is 
used in deriving the aerodynamic loads. 


2.2 COORDINATE TRANSFORMATIONS 


2.2.1 Rotating to Nonrotating Transformation 



cos if/ sin 0 
— sin i l> cos i p 0 
0 0 1 



(2.1) 
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2.2.2 Coned to Unconed Transformation 



ns /? 



and its inverse: 



cos P p 

0 

— sin 

0 

1 

0 

sin P p 

0 

cos 0 p 



(2.2) 


(2.3) 


2.2.3 Deformed to Undeformed Transformation 

From Ref. 56: 



1 

- ( w s*t> + V ,x) 

- (w x - <t>v x ) 



(2.4) 


and its inverse: 



1 -{v tX + <f>wj - (w* - <t>v x ) 
-(4> + v .x w J 


,x 


w 


yX 


l 

<t> 


1 



(2.5) 
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2.2.4 Deformed and Rotated to Deformed Transformation 



1 0 0 
0 cos (f> - sin (f> 
0 sin <p cos (f> 



( 2 . 6 ) 


2.2.5 Undeformed to Deformed and Rotated Transformation 
By setting <p equal to zero in Eqn. 2.5 we obtain: 



(2.7) 
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Chapter III 

METHOD OF SOLUTION OF EQUATIONS OF MOTION 

3.1 FREE VIBRATION SOLUTION AND MODAL 
TRANSFORMATION 

The natural frequencies and mode shapes of the blade are calculated using 
a local Galerkin Method of Weighted Residuals. The method, as presented in 
Ref. 8, has been modified to allow for modeling of articulated rotor blades. 
The following description uses the notation of Ref. 8. The finite clement used 
has 1 1 degrees of freedom: displacement and rotation in lag at the end nodes 
(4), displacement and rotation in flap at the end nodes (4), and torsional rota- 
tion at the end nodes and a mid-element node (3). Displacements inside the 
element are in terms of cubic Hermitian interpolation polynomials for flap and 
lag, and quadratic Hermitian interpolation polynomials for torsion. The axial 
displacements are eliminated by making the assumption that the blade is in- 
extensional. The Galerkin method requires that the sum of the weighted resi- 
duals of both the differential equations and the natural boundary conditions 
be equal to zero. These calculations are based on the linear, homogeneous, 
undamped equations of motion of the blade in a vacuum, which are obtained 
from the nonlinear partial differential equations of motion given in Eqns. 4.3 
through 4.6. The equations of motion for the elements used in modeling the 
blade, where superscript e indicates the element number and subscript 1 indi- 
cates that these are linear matrices, are: 
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C/f]{« e > + ([*?] + [A'fj + [ T^){a e } = 0 e = 1,2 ,...,£ 


(3.1) 


The E equations 3.1 are assembled and root boundary conditions are imposed 
by dropping the rows and columns of the system mass and stiffness matrices 
corresponding to the degrees of freedom restrained at the root. These degrees 
of freedom will vary depending on whether an articulated or hingeless blade 
is being modeled. 

In forward flight the blade root pitch varies periodically with azimuth. This 
causes coupling between the modes and periodicity in the mode shapes. For 
this study the modes do not vary with azimuth and are calculated for a con- 
stant root pitch angle equal to the collective pitch. The coupled mode shapes 
of the blade are normalized by dividing each mode shape by its maximum tip 
displacement in either flap, lag, or torsion. 

A modal coordinate transformation is introduced to reduce the number of 
unknowns in the problem and to assemble the various element matrices into 
the system matrices and load vectors. The modal transformation has the form: 

{a e } = LQ e m (3.2) 

The vector {y} of generalized coordinates become the new problem unknowns. 
If m modes are used to perform the modal coordinate transformation, then (y) 
is of length m. The matrix [Q e ] is of size 11 by m. The columns of this matrix 
contain the portions of the normal mode eigenvectors corresponding to the 
modal degrees of freedom for the given blade element. 
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The stiffness matrix [A\y,y)] of the complete blade is assembled by sum- 
ming the stiffness matrices of the individual blade elements, each of which has 
had the modal coordinate transformation applied to it: 

E 

[k<7 J)] = W < 3 - 3 > 

e= 1 

where [Alysy)] is the stiffness matrix of the e-th finite clement in the local 
coordinate system. The blade mass matrix [A/(y,.y)], damping matrix 
[C(y,y)], aerodynamic load vector and nonhomogencous load vec- 

tor {F}, which contains forcing terms independent of the generalized coordi- 
nates, are assembled similarly. This assembly results in the following set of 
nonlinear, coupled, ordinary differential equations which, for the case of for- 
ward flight, have periodic coefficients: 

[A/(y ,y)l{y) + ic(y J)]{y} + [*(y J)](y} + + {A} = o (3.4) 

3.2 SOLUTION METHOD FOR HOVER 

The solution method for the case of hover is similar to that of Ref. 8 but 
has been modified so that the aeroelastic response and the rotor trim condition 
are calculated simultaneously. For hover the nonlinear equations of motion 
have constant coefficients. The generalized coordinate vector {y} is written as 
the sum of a constant vector {y 0 } and a perturbation vector (Ay(r)}: 

{y} — (go) A- (Ay(/)} (3.5) 
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Substituting Eqn. 3.5 into Eqn. 3.4 results in a system of nonlinear algebraic 
equations: 


[KO'o.GXIO'o} + M(yo>°)} + {Z 7 } = 0 (3.6) 

and a system of small perturbation linear ordinary differential equations: 

[A/(y 0 >0)]{4P(')} + [C(7 0 ,O)]{Ay(r)} + [Z6T 0 , 0)]{A><r)} = 0 (3.7) 


where: 


[C(y 0 , 0)] = [C(y 0 , 0)] + 


M=y 0 J = 0 


(3.8) 


[*(vo. o )] = [*(yo> o )] + |/o/ T 

Z = I *— * — * V = LU V 


y =y^y = o 


+ 


aw 

aw 


'y = yo’ y = o 


(3.9) 


Since the aerodynamic load vector {A(y,y)} is calculated numerically, the 


~d{A}l 

and 

~d{A}~ 

d{y) J 


L d{y} J 


aerodynamic stiffness and damping matrices | | and | | must be 

computed using a numerical finite difference approximation^]. 

In this study an additional equation is appended to Eqn. 3.6 requiring that 
the rotor be trimmed in the collective so that the rotor thrust equals the given 
helicopter weight. The modified equation becomes: 


[KW.OBO'o'} + {^'(vo .0)} + i F ') = o 


(3.10) 


where: 
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(3.11) 


W - 



[*W,o)] 


K(y 0 ,O) 

0 


(3.12) 


{A’(y 0 , 0)} = 


^(vo'.O)l 

C-0 o ,O) J 


(3.13) 


{F'} = 



(3.14) 


Equation 3.10 is solved for {y' 0 }, the static nonlinear blade equilibrium po- 
sition, using a Ncwton-Raphson iterative procedure. The linearized stability 
is determined by solving the standard eigenvalue problem defined by the first- 
order, state variable form of Eqn. 3.7: 




JO] 

[A/(y 0 ,O)] -l [C(y 0 ,O)] 


jn 

-CAKy 0 .o)r 1 [^ 0 .o)] 




{AyJ 


= [*^(Toi 0 )] (Ay i ) 


(3.15) 


where: 



(3.16) 


Blade stability is determined by the eigenvalues of [S(y 0 , 0)]. The eigenva- 
lues occur in complex congugate pairs. For m modes used in the modal coor- 
dinate transformation: 
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^j— Cy ± iojj j = 1,2 


,m 


(3.17) 


The blade is stable if all £y < 0. 


3.3 TRIM PROCEDURE 
3.3. 1 Introduction 

A helicopter in straight and level flight, flying at a constant velocity, must 
satisfy force and moment equilibrium. Determining the control inputs neces- 
sary for this equilibrium condition is referred to as trim analysis. In this study 
two trim procedures are used. The first is a flap trim procedure where the 
blade flexibility included in the trim analysis is limited to the first flap mode 
and all other blade deformations are neglected. This trim procedure is used 
principally to compare with results of Ref. 8 which also used this procedure. 
The second trim procedure was developed specifically for this study and in- 
cludes the blade flap, lag, and torsional flexibility in the trim analysis and thus 
enables one to have a much more general model for the blade. The results 
from both these trim procedures are used subsequently as trim inputs to the 
full aeroelastic model. The full model produces stability and response infor- 
mation on the blade motions, including higher modes and higher harmonics 
of the solution which were neglected in the trim procedure. 
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3.3.2 Flap Trim 

The flap trim analysis used in this study was developed in Ref. 19 and it 
was used in generating all the results of Ref. 8. This flap trim is a propulsive 
trim procedure based on the following assumptions: 

1. Blade flexibility is modeled with a linear flapping equation of motion for 
the first flap mode. 

2. The blade can have a precone angle /? p , a constant offset between the 
elastic axis and aerodynamic center x A , and built-in twist with linear 
spanwise variation d B . 

3. Quasisteady aerodynamics including reverse flow are used. Stall and 
compressibility effects are neglected. 

4. Rotor shaft dynamics are neglected. Tail rotor and rotor speed vari- 
ation effects are neglected and so yawing moment equilibrium and lat- 
eral force equilibrium are not enforced. 

5. Vertical force equilibrium, longitudinal inplane force equilibrium, rotor 
pitch and roll moment equilibrium, and an inflow equation are enforced 
by with appropriate settings of the collective pitch 0 O , cyclic pitch com- 
ponents 6 U and 0 lc , the inflow ratio /, and the rotor angle of attack. a R . 
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3.3.3 Full Propulsive Trim 


3.3.3. 1 Assumptions 

A helicopter in free flight has six degrees of translational and rotational 
freedom. Consequently equilibrium is guaranteed if three force and three 
moment equilibrium equations are satisfied. In this trim procedure yawing 
moment equilibrium and lateral inplane force equilibrium are not enforced and 
consequently the tail rotor pitch setting and the main rotor shaft angle in the 
lateral plane are not considered as trim variables. The helicopter is assumed 
to be in straight, steady flight at constant speed. Dynamics of the rotor shaft 
are neglected. Details of the blade modeling are given after the next section. 

3. 3.3.2 Trim Variables and Equilibrium Equations 

The five trim variables are the inflow ratio X , collective pitch angle d 0 , cyclic 
sine component , cyclic cosine component d u , and the rotor angle of attack 
a R . 

The five trim equilibrium equations are: 

1. The inflow equation. This is a momentum theory relationship between 
the inflow ratio X, advance ratio p, rotor angle of attack a R , and thrust 
coefficient C T . 


FT( 1 ) = 0 = Cj + 2 yj p* + A 2 (p tan a R — X) 


(3.18) 



2. The rolling moment equation. Since the tail rotor and the rotor angle 
in the lateral plane are not modeled, the rotor rolling moment is just set 
to zero. 

FT(2) = 0 = M r i (3.19) 


3. The pitching moment equation. The rotor pitching moment is count- 
ered by moments due to fuselage drag and the weight of the helicopter. 
All moments are taken about the rotor hub. The drag of the fuselage 
is given by: 

D= 'Ap A V 2 f 


where f is a parasite drag area. In nondimensional form this is: 

c Df = 'V-2- 

f 

where A is the rotor area. Typically — ~0.1. If the fuselage drag is 

A 

assumed to act a nondimensional distance h from the rotor hub, then 
the drag causes a nose down moment of: 


— 2 f ~ 

M dr = — '/ip — h cos a 


The weight of the helicopter, also assumed to act a distance h below the 
hub, produces a nose up moment due to its horizantal offset from the 
hub by a distance h sin a R . The sum of these two moments and the ro- 
tor moment, as shown in Fig. 3.1, must be zero. 
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(3.20) 


— f — — 

FT(3) = 0 = M ( - Vin"-—h cos a^ + ZzC^/sina^ 

/ 1 

4. The vertical force equation. The weight of the helicopter is in equilib- 
rium with the components of the rotor thrust and rotor longitudinal in- 
planc force vectors as shown in Fig. 3.1. 

FT(4) = 0 = Cj - cos of + C f-f sin ot ^ (3.2 1 ) 

5. The longitudinal force equation. The drag of the fuselage is countered 
by components of the rotor thrust and rotor longitudinal inplane force 
vectors as shown in Fig. 3.1. 

FT(5) = 0 = C T sin ot^ — C ^ cos oc^ — Vi /j. — (3.22) 


3.3.3. 3 Aeroelastic Model for Trim 

The aeroelastic model used in the trim procedure is a simplified version of 
the model used in the actual aeroelastic model of this study. This simplified 
model retains all the geometric, stiffness, mass, and aerodynamic modeling 
capabilities of the main analysis however the number of modes used in the 
analysis is reduced to four and furthermore only the constant and first and 
second harmonics in the Fourier series representations of the generalized co- 
ordinates are retained. In addition the mode shapes are calculated using only 
two finite elements. This reduces the cost of a response calculation using 
quasilincarization by more than an order of magnitude but retains the accu- 
racy of the constant and one per rev portions of the solution which are re- 
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quircd for determining the trim state of the rotor. Harmonics of the solution 
beyond the second are not modeled and so this simplified model is not suitable 
for calculation of the vibratory loads. 

3.3.3.4 Method for Determining the Trim State 

In order to trim the rotor, values for the trim variables X, d 0 , 8 ]s , d ic , and 
a R must be found such that the five trim equations are simultaneously satis- 
fied. The problem then is to solve the set of nonlinear algebraic equations 

FT(A,0 o ,0 lj ,0 lc ,a^)=O (3.23) 

where FT is the vector of trim equations of length 5 each of the elements of 
which is implicitly dependent on the trim variables. Equation 3.23 is solved 
by a packaged 1MSL procedure called ZSPOW using a modification of the 
Powell Hybrid Method. This involves calculating the residual values of the 
components of the vector FT due to an initial estimate of the response, and the 
Jacobian of the system at this point, then using this information to make sub- 
sequent moves toward a minimum of the residuals in the five trim equations. 

The procedure ZSPOW calls an external subroutine written for this study 
which, given the five trim variables, calculates the residuals in the five trim 
equations. The computation of these five trim equations is described next. 

1. ZSPOW provides a vector T containing the trim variables, either with 
one element incremented when calculating the Jacobian (i.e. the matrix 
of first derivatives of the vector FT with respect to the five trim vari- 


25 


ables), or with all five elements changed when moving to a new solution 
estimate. 

2. The coupled mode shapes of the blade are dependent on the collective 

angle 6 0 , so if this angle has changed, the subroutine recalculates the 
mode shapes. 

3. The mode shapes and the trim settings are used to determine a con- 
verged response for the blade using quasilincarization. 

4. The new response is used to calculate the constant portion of the rotor 
hub shears and moments as described in Ch. 5. 

5. The shears and moments are used in the trim equations to determine the 
trim equation residuals for this trim condition. 

6. The trim equation residuals are returned to the algebraic equation sol- 
ver which attempts to drive them all to zero. 

A flow chart showing this procedure is given in Fig. 3.2. 

Typically the algebraic equation solver calculates six responses to determine 
the starting condition and the Jacobian of the system. Subsequently 5-10 ad- 
ditional iterations are needed to obtain a trim solution which is accurate within 
three significant figures in terms of each of the trim variables. 
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3.4 SOLUTION METHOD FOR FORWARD FLIGHT 


3.4. 1 Quasilinearization 

The quasilinearization procedure presented in Ref. 8 is used to determine 
periodic response and stability in forward flight. The equations of motion are 
written in first order form: 

Cq) = [z#)]{<7> + {<?}; t)} t)} 


The state vector {q(i/0) is defined by: 


toO/O} = 


bin) 
{ym J 


(3.25) 


and is of length n = 2m where m is the number of modes used in the modal 
coordinate transformation. The generalized coordinates (yCA)} represent the 
time varying contributions of each of the n modes to the periodic motion sol- 
ution. The matrix L contains the linear portions of the blade model, N con- 
tains the nonlinear terms and Z contains the nonhomogeneous terms. All 
aerodynamic contributions are contained in the vector A which is calculated 
numerically. 

Quasilinearization is based on a first-order Taylor series expansion of Eqn. 
3.24 about an approximate solution 




mu) 


-i k 


(W* +l -{*>*) 
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(3.26) 


~\k 


+ 


Wnl) 

L d{}} J 


({q) k+ '-{q} k ) 


where: 


“ HF NL ) " 

- mw] + [ 1 + 

~d{A) " 

d{q) 

L m J 

_ d{q] _ 


" Wm. > y _ 

[«T + 

r d{A] i 

_ _ 

L d{q) J 

i — 

cu 

s* 

1 


Equation 3.26 can be written as: 

{i>* +l -W)]*{«>* +1 +{/#)}* 


where: 



(3.27) 

(3.28) 

(3.29) 


(3.30) 


Equation 3.29 gives an iterative procedure, the sequence of solutions of 
which converge to the solution of Eqn. 3.24. Again, since the aerodynamic 
load vector ^({<7}, {<7}; i/^)} is calculated numerically, the derivatives 


d{A} 


d{q) 


and 


d{A) 


L d{<7) J 


must be computed using a numerical finite difference 

q) J 

approximation)^]. 

The response solution determined by this quasilincarization procedure is 
stored in the form of a constant and five sine and five cosine Fourier coeffi- 
cients for each of the six generalized coordinates in the vector {y} . Wherever 
in the program values of the components of {y} or {<7} are needed, these are 
calculated for the given azimuth 1 \t from this Fourier series. 

In order to solve Eqn. 3.29, the state transition matrix [<t>(27r )]* of the ho- 
mogeneous system: 


{£}*+' = 


(3.32) 


at the end of one period is needed. The matrix [0(27:)]* is computed using a 
single pass version[8] of the rc-pass algorithm[51] in which the columns of 
[<£(271)]* are obtained by solving Eqn. 3.32 n times, with initial condition vec- 
tors {<7(0)} fc+l k = 1,2 which have all their elements equal to zero, except 
for the i-th, which is equal to one. 

The initial vector (<y(0)} * +1 for Eqn. 3.29 is determined by requiring that the 
solution be periodic with period 2n. The vector {^(0)}* +1 is shown[8] to be the 
solution of the linear algebraic system: 

([/] - [<D(27t)]*){ ? (0)}* +1 = «2 k)}* (3.33) 
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where [I] is the identity matrix and {^(27 t)}* is the solution vector of Eqn. 3.29 
integrated from i/< =0 to \J/ - 2n with initial conditions {<7(0)}* = {0}. 

System stability, linearized about the time dependent equilibrium position 
{< 7 ( </<)}*, is determined from Floquet theory[19]. The stability of the system is 
governed by the characteristic exponents obtained from the transition matrix 
[<t>(2/r)]*. For the characteristic exponents 

t-cj = ± io) cj j=U2,...,m (3.34) 

the linearized system is stable if all ( cy < 0. 


3 . 4.2 Convergence 

A convergence criterion must be established to determine when the response 
solution from one iteration of quasilinearization is sufficiently close to the sol- 
ution obtained in the previous iteration so that the procedure can be termi- 
nated. The response solution at the end of an iteration is saved in the form 
of a Fourier series for each of the six generalized coordinates in the vector 
(v). The magnitude of the harmonic, i.e. the square root of the sum of the 
squares of the sine and cosine conponents of that harmonic, can be calculated 
for each harmonic. A 6x6 matrix, l M can then be constructed having as its 
six columns the magnitudes of the constant portion and first five harmonics for 
each of the six generalized coordinates. 

In Ref. 8 a convergence criterion is used which is based on a sum of the 
absolute values of all the harmonics. From the response matrices l M and 2 M 
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obtained for two consecutive response iterations, the following sums are cal- 
culated: 


6 6 



‘= 1 J= 1 


6 6 

-ZZ' 

‘= 1 h 1 


A/:; 


(3.35) 


If S 2 is not zero, convergence is assumed to have occured when: 



*100< DELTA 


(3.36) 


If S 2 is equal to zero, convergence is assumed to have occured when: 

| S 2 — Sj | * 100 < DELTA (3.37) 


In Ref. 8 the convergence control error parameter DELTA was set to 0.5. This 
required that the change in the sum of the absolute value of the harmonics 
from one iteration to the next be less than V. i% . 

The higher harmonics in an equilibrium solution have much smaller mag- 
nitudes than the constant and lower harmonics. Therefore the convergence 
criterion given above weighs more heavily the importance of the lower har- 
monics. When convergence has been achieved by this criterion, there can in 
fact still be large differences in higher harmonics of the solution. Another 
problem with this criterion is that occasionally, purely by chance, the sum of 
the harmonics from one iteration can be very near to that obtained from the 
previous iteration, without true convergence being achieved. This can lead to 
erroneous values for higher frequency phenomena such as vibratory hub shears 
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and moments. Results illustrating these problems and their solution are de- 
scribed in Ch. 7. 

In this study a more stringent convergence criterion was developed based 
on a norm of the absolute differences in harmonics between two successive it- 
erations. The matrices ] M and 2 M are calculated as before and a norm is 
defined as: 


NORM = 


6 6 


ZBH- V 

/=1 i=j 


Vz 


(3.38) 


This norm moves smoothly to its optimal value and has the additional advan- 
tage that it cannot produce a small number by chance, but only when two so- 
lutions are actually close in each harmonic. Convergence has occurred when 
NORM < N c , where N c can be varied to give different levels of accuracy in the 
harmonics. 

It is not sufficient just to have a valid convergence criterion if progress is 
not being made toward convergence. If sufficiently tight error bounds are not 
used to solve the differential equation in quasilinearization, the higher har- 
monics of the solution will not converge and will instead vary from one iter- 
ation to the next in an apparently random fashion. If a loose convergence 
criterion is being used this can create the appearance of a spurious converged 
solution. 

In this study, tighter local error bounds were used than in Ref. 8 and in 
addition bounds were always tightened systematically as the iterative process 
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continued. This led to more consistent values for the hub shears and moments 
used in the study. 
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Chapter IV 


INCORPORATION OF UNSTEADY AERODYNAMICS IN EQUATIONS 

OF MOTION 

4.1 EQUATIONS OF MOTION 

In this study finite-state arbitrary-motion aerodynamics are incorporated in 
the coupled flap-lag-torsional equations of motion which are presented in Ref. 
57 and are also given in this section for the sake of completeness. The geom- 
etry of the problem is shown in Figs. 2.2, 2.3, and 2.4. The assumptions used 
in deriving these equations are concisely summarized below. 

The blade, in its undeformed state, is straight and has no angle of sweep or 
droop and no torque offset. The blade is cantilevered or hinged at the hub and 
has a blade root offset of e, from the axis of rotation. The feathering axis, 
which coincides with the blade elastic axis, is preconed by an angle (i p which 
is assumed to be small. The blade has an angle of built in twist of d b {x) which 
is assumed to occur about the undeformed elastic axis. Collective, cyclic, and 
higher harmonic pitch controls also contribute to the total geometric pitch an- 
gle which is given by: 

dd'l')=d 0 +d lc cosiP+6\ s sm $ + 0////c#)+ 9 b (x) (4.1) 

The cross section of the blade is symmetrical with respect to its major prin- 
cipal axis and the tension center and the elastic center are coincident. The 
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cross sectional center of grav ity, aerodynamic center, and elastic axis need not 
be coincident, and values of these parameters can change over the blade span. 

The blade is made of a homogeneous, isotropic, linearly clastic material. 
The Eulcr-Bcrnoulli assumption, that cross sections remain plane and normal 
to the deformed clastic axis during deformation, is made. This neglects the 
effects of shear deformation. The blade can bend in two mutually perpendic- 
ular directions. Elastic torsional deformations occur about the deformed clas- 
tic axis while additional deformations due to the flexibility of the control 
system occur about the feathering axis. 

The blade is assumed to have moderate deflections, which implies small 
strains and finite rotations or slopes. The elastic rotations arc assumed to be 
of the order s where e < 0.20. Terms containing the squares of the slopes are 
neglected when compared to terms of order one, i.e. 

0(1) + 0(e 2 ) =0(1) (4.2) 

Structural damping forces are assumed to be of a viscous type. This is a 
reasonable assumption since modern rotor hubs are equipped with elastomeric 
dampers. No coupling exists between the blade dynamics and a fuselage. The 
rotor is rotating at a constant angular speed Q and the helicopter is in steady 
trimmed level flight. 

The final nonlinear partial differential equations of motion, including iner- 
tia loads but with aerodynamic loads left in symbolic form, are given by Eqns. 
53, 54, 55, and 56 of Ref. 57. These, with several typographical corrections, 
are given below: 
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Axial Equation 


T ,x + Pxi = 0 (4.3) 

Lag Equation 

[£(/ 2 cos 2 6 g + / 3 sin 2 0 G )v ^ + £(/ 2 -/ 3 )sin 0 G cos 0^w xx - 2v xx <f>) 

+ E(I 2 -h) ( f )W ,xx cos2d G],xx + ( GJ< t>,x w ,xx),x-( v ,x T ),x 

+ [20^ G (/ m2 -/ w3 )sin 0G cos 0G],x _ [ w ^ 2 - x /(- x o+ f ’i) cos 0G],jc 

+ mv + 2mfiw-2wQ/? p w-mx { Q^0 G sin $ G -mQv - mQxj cos 0 G 

+c lSiy~PYA = Q (4.4) 

Flap Equation 

[£(/ 2 -/ 3 ) sin 6 G cos 6^v xx +2 (f>w xx ) + E(f 2 -I^v ^ cosIQq 
+ E{1 2 sin 2 6 g +/ 3 cos^vv^]^ - (GJ4> x v xx )^ - (w^T) x 

- [mxp 2 (x 0 +^i)sin 0^ + [2D 2 S G (/ m2 sin 2 0 c +/ m2 cos 2 6 G )} X 
'f * 

+mxfl l e G cos e G + mw + 2 mQp p v + mQ 2 p^x 0 +<?j) 

+ q S F w -PZA = ° (4.5) 

Torsion Equation 
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x +v iXX w iX )] iX +E(I 2 -/ 3 )[ sin e G cos 0^v 2 xx -w 2 xx ) 

~ v ,xx w ,xx co< ' ~QS7^ -& 2 {( f m2 + + 4>) 

* * * * 

cos 0 ^[vv >jc (xq + ^j) — ( Xq +e\)v tX 4> ~ 4> v + V( P + w 
+ ^p( x 0 + e i) +2v(w >JC + (Ip)] +mxj sin 8q\_ — v^-Xq +ej) — v + v] 

+ I m 2 [( I + 2v x ) sin 6 G cos 6 G + ( - p p v x + 2v x <p + 0 +w x V r ) cos 2 0 g 

+ (2w ^ — </> ) sin 0 G ] E I — (1 + 2v x ) sin 6 G cos d G (f> sin 6 G 

+ (2 v x w x + *v> x +2 w x -20v x - 0 + P p v x ) cos 2 0 g ]} + q XA = 0 (4.6) 

4.2 DERIVATION OF UNSTEADY AERODYNAMIC LOADS 
4.2. 1 Introduction 

The finite-state arbitrary-motion aerodynamics used in this study were de- 
veloped originally in Ref. 13. In that work a generalized version of 
Greenberg's unsteady aerodynamic theory[24] was developed for the unsteady 
lift and moment on an airfoil undergoing arbitrary pitch and plunge motion 
about a steady pitch angle in the presence of time-varing oncoming velocity 
and variable inflow. Non-circulatory lift and moment were determined from 
the unsteady linearized Bernoulli's equation for the pressure on the surface of 
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the airfoil. A circulatory flow is calculated to enforce the Kutta condition at 
the trailing edge. After enforcing the Kutta condition, the circulatory lift is 
calculated in the Laplace domain. 

A generalized Theodorsen lift deficiency function is found which acts as a 
Laplace domain operator between the V* chord downwash velocity and the 
circulatory loads divided by the time-varing oncoming velocity. A second or- 
der rational approximation for the generalized Theodorsen lift deficiency 
function is used to obtain a finite-state time-domain representation of the cir- 
culatory loads. The time-domain circulatory loads are finally expressed in 
terms of the airfoil degrees of freedom and two additional state variables de- 
noted augmented states. The two state variables are governed by a system of 
ordinary differential equations driven by the V* chord downwash velocity. 


4.2.2 General Arbitrary- Motion Lift and Moment Expressions 

From Ref. 12, the general noncirculatory and circulatory arbitrary-motion 
airloads on an airfoil, such as that shown in Fig. 4.1, are given by: 

0 = \p A a fP R) 2 {^kt) + U 7 {t)M{t) + f> 7 (0[a 0 + Aa(/)J x 

- (4.7) 

M n^) = \pA a fP 'Rfv-Aty&t) + y p A afb R)\{d - Vi)U T Aa(t) 

+ a\_Ah(t) - Ufl)] + d Uj{t)loi 0 + Aot(f)] - bR{ Vs + d 2 )M(t)} (4.8) 
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Ur i 


LJt) = p 4 aib R)U-,{t)[0M6n5{-^) 2 XAt) + 0.10805(-^-)A\(')] 

bR bR 


U 


TO 


+ ±- PA aj ( bR)U l (t)Q{t ) (4.9) 

M c (t) = bR(a + Vi)L£t) - j-p A agb R) 2 U 7 (t)Q(t) (4.10) 

where the two underlined terms cancel each other and Ujq is the constant 
portion of U T . From Fig. 4.1 it can be seen that the 3 A chord downwash ve- 
locity Q(t) will have contributions due to oncoming velocity U 7 (t), vertical ve- 
locity A h(t) — U p (t), and angular rotation Aa(f): 


QHt) = Uj{t)l<XQ + Aa (/)] + [AA(0 - U p (t)-] + bR( V* - «)Aa(0 (4. 1 1 ) 


The augmented states X\ and X 2 are governed by a first order differential 
equation: 




0 

-0.01365(- 


1 


U 


to a 


u 


TO 


-r~~) — 0.3455(-=r— ) 

bR bR 



(4-12) 


4.2.3 Lift and Moment in Terms of Local Blade Velocities 

To apply Eqns. 4.7 through 4.10 to the present problem, it is convenient to 
consider, as was done in Ref. 57, a system of unit vectors , (? y , and ^ , 
where the double prime system has been rotated by an angle of — 0e' from the 
deformed blade coordinate system. The transformation between the unit vec- 
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tors is given by Eqn. 2.6. From this transformation the velocities in the two 
systems are related by: 


U”= U^cos<t>- U' z sin0 


(J" z = Uy sin </>+£/' cos <p (4. 1 3) 

The velocities used in this derivation will be velocities of the clastic axis, as 
were used in Ref. 12, rather than velocities relative to the elastic axis, as were 
used in Ref. 57. With this convention, the following relations can be identified: 

U T =Uy, <x 0 + Aa =d G + (p , U p -Ah=U z (4.14) 

By substituting the above velocity terms in Eqns. 4.7 through 4.10 and us- 
ing the following notation: 


d „ d ~ p 

Tr^^TT’ X A = ~ — > <3 = (xa — Vi) , R = — 
dt dip bR € 


U 


U = £_ n — Y — ^ X F — — 

y ns 1 z ~ ns ’ Ai ~~ x ^ 


^nc 


L 


nc 


a iPA {bR\nS) 


2 ’ 


M nc = 


M 


nc 


a iPA {bR)\nS ) 2 


L c = 


a iPA {bR\nS) 


2 ’ 


_ M. 

M=— c 


ap A {bR)\nS) 2 


(4.15) 


the noncirculatory lift and moment become: 


* 4 


4 * 44 4 


A D * __ 

L nc = —{vpc + 4 >)-ir z - bR(x A - Vz\e G + <p)+ u” y {e G + <p)} ( 4 . 1 6 ) 
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(4.17) 


* * 


U„c=~ + *)- &A - '/>)^+ - ' -K^C + *) 


* * * * 


- />/?(3/8 + x\ - x A XO G + <f > )} 


and the % chord downwash velocity becomes: 


Q(t) = im G + <t>)-u ” 2 + bR(i-x A xd G + 4>) 


The circulatory lift and moment can be expressed as: 


_ _ 2 u y0 - 

L= L^[0.006825(-=-) a ,(i^) +0.10805(-=-)A' 2 (^)] + '/ 2 U”Q 
y bR bR y 


M c = x A L c 


where the nondimensional augmented states are now governed by: 


JL 

X~> 


1 


U‘ 


-0.0 i 365( 


yb ^2 


U‘ 


-) z — 0.3455(-==r) 
bR bR . 




The circlatory lift can be written in a more convenient form if a 
associated with unsteady effects is defined as: 


u; 0 7 _ u; 0 _ 

H = [0.006825(-r=-) 2 A’i(^) + 0.10805(^=-)A' 2 (i/')] 
bR bR 


with this expression the circulatory lift becomes: 


L c = U” y [H+ */ig] 


(4.18) 

(4.19) 

(4.20) 

(4.21) 
quantity 

(4.22) 
(4-23) 
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4.2.4 A Previous Formulation for the Circulatory Lift 

In Ref. 12, the circulatory lift and moment were written in a different form. 
As will be shown, this introduced a spurious instability at higher advance ra- 
tios which was a consequence of the particular details of the formulation, and 
was not an inherent property of the time domain aerodynamics. The circula- 
tory lift and moment were expressed in terms of their quasisteady values mul- 
tiplied by a time domain lift deficiency function, leading to the following 
cxpressions[12]: 


r = — r = [ 'A +F]U"Q 

PA a,(bRxno 2 


(4.24) 



M c 

— 2 7 = v * +nu;Q 

p A af,bR)\Qi ) 2 y 


(4.25) 


Where the function F was a lift deficiency function in the time domain cal- 
culated using the V* chord downwash velocity corresponding to a typical blade 
section located at the three quarter span station of the blade. The quantity F 
was assumed to be a function of time only, and valid for the entire blade span: 
^ Vs 2 » Qts) = 


[0.006825( 


U 


TSO \2t T7 


U 


-nx TSI + mom^-nx TS2 

£/v bR 

Qts 


(4.26) 


where: 


x ts ~ nondimensional radial position of the typical section 


42 


Qrs = Q evaluated at x 0 — x TS 


U rso = x TS + e = constant portion of LT y evaluated at the typical section. 


The differential equation governing the augmented states was given by: 



° _ 1 
-0,01365 — 0.3455(— 


bR 


bR 



(4.27) 


The purpose of this formulation for the lift and moment was to simplify the 
solution of the acroclastic response and stability problem by having the aug- 
mented states governing the unsteady aerodynamics be only functions of time. 
Thus the augmented states were independent of the spanwisc locations of the 
cross sections of the blade. While the formulation was successful in accom- 
plishing this objective, it also introduced a spurious singularity in the circula- 
tory loads at high advance ratios. 

For typical rotor blades, the V* chord downwash velocity Q can become 
zero on the advancing side of the rotor at sufficiently high advance ratios. 
From Eqn. 4.26 it can be seen that when Q TS approaches zero the lift deffici- 
ency function F will approach infinity. If the circulatory lift is being calculated 
at the typical section, this produces no algebraic problem because from Eqn. 
4.24, F is multiplied by Q which cancels with the denominator in Eqn. 4.26. 
If the circulatory lift is being calculated anywhere but at the typical section, 
Q and Q rs do not cancel and as Q TS approaches zero the circulatory lift ap- 
proaches infinity. 
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This singularity in the circulatory lift leads to a spurious instability in the 
flap degree of freedom at an advance ratio of 0.45 in the flap-lag study pre- 
sented in Ref. 12. When this version of the formulation of the unsteady loads 
is used, the spurious instability is also reproduced in the present flap-lag-tor- 
sional study. With the correct formulation of the circulatory loads, given in 
Eqns. 4.18 through 4.23, no such instability is found. A comparison of results 
due to the two formulations is given in Ch. 7. 


4 . 2.5 General Expressions for Aerodynamic Forces and Moments 


Forces and moments per unit span in the deformed blade coordinate system 
must now be found. The loads will have components due to circulatory lift, 
which is assumed to act perpendicular to the resultant velocity at the elastic 
axis, noncirculatory lift, which is assumed to act perpendicular to the blade 
chord, and an additional profile drag force which is parallel to the resultant 
velocity of the elastic axis. These orientations are shown in figure 4.2. Drag 
is calculated using: 


D = Pa V 2 b rc D0 > where: U 2 = U” 2 + If 2 

»= z ? t = ^v” 2 +v*?) 

p A a-pR{M) 2 ' 


(4.28) 


Using these assumptions concerning the direction of the lift and drag forces 
(per unit span), the loads per unit length in the double primed system become: 
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(4.29) 


\J* 

F' za = L c cos[ tan (-=r-)] + L nc cos(0^ + </>) 

F’ ya =-L c sin[ tan - '(-^-)] - T„ c sin(0 G + <t>)~ D (4.30) 

^ it 


or: 


= L c cos ( 0 C7 + <t>) 

sjuf + 

= - L c 7= _ - + 0) - -%^ 2 + 

s/uf+uf 


Qxa ~ :Vf c + V/ 


rtc 


(4.31) 


(4.32) 


(4.33) 


To obtain the forces and moments in the undcformcd coned blade coordi- 
nate system the transformations given by Eqn. 2.7 must be used. Using these 
transformations the loads are: 


p = 
1 xa 

- v p* p* 

v ,x i ya YV ,x 1 za 

Qxa Qxa 

P = 
£ ya 

P" -v w F 
1 ya v ,x ,x £ za 

C 

II 

x * 

Cj 

P = 

1 za 

F 
1 za 

Qza vv ,Aa 


(4.34) 
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4.3 NUMERICAL IMPLEMENTATION OF AERODYNAMIC LOADS 


4.3.1 Introduction 

The values of the expressions for the aerodynamic loads given in Eqn. 4.34 
must be calculated at various azimuthal positions and provided as input to the 
acroelastic analysis. 

In this study two approaches to the overall problem were implemented. In 
the first approach an algebraic manipulation program was used to derive al- 
gebraic expressions for the local velocities. Numerical values for the loads were 
then computed from these implicitly derived values for the local velocities. 
This approach avoids the neccessity of introducing assumptions which would 
allow one to neglect terms to simplify the expressions. In the second approach, 
assumptions were made for the orders of magnitude of all quantities in the 
expressions. Using these assumptions, higher order terms were systematically 
neglected leading to explicit expressions for the aerodynamic loads. This sec- 
ond approach is the same as was used in Refs. 57 and 8. It was implemented 
for comparison with the implicit formulation to verify that the two approaches 
gave similar results, which they did. These two approaches are described in 
the following sections. In both cases the basic approach to finding the local 
air velocities in terms of blade displacements and the other problem parame- 
ters follows the approach used in Ref. 57. 


46 



4.3.2 Implicit Implementation 


To numerically calculate implicit values for the blade loads, v alues of the lift 
and moment expressions in Eqns. 4.16 through 4.23 arc determined and then 
these values are numerically substituted in Eqns. 4.31 through 4.33. Subse- 
quently these results are then numerically substituted in Eqn. 4.34 to give the 
final implicit loads. 

The lift and moment quantities in Eqns. 4.16 through 4.23 are in terms of 
local blade velocities in the double primed coordinate system. Expressions for 
these velocities in terms of blade slopes and displacements and other problem 
parameters must therefore be found. The local velocities are due to two com- 
ponents, \ Ea , the velocity of the elastic axis, and V, , the velocity of airflow 
due to inflow. The velocity of the elastic axis may be found from: 

V EA ~ R EA + ^ xR £/l (4.35) 

where R £/) and Q are given by: 

Rea ~ e i l + (*o +v ^y +w ?z (4.36) 

a = Qk (4.37) 

A 

Using the transformation of Eqn. 2.2, R £y4 and R £/1 can be expressed in the /, 

A A 

j , k coordinate system as: 

A 

r ea = i e i + C*o + “) cos P p ~ w sin Pp \ 1 
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(4.38) 


+ ‘7 + C(A- 0 + u) sin P p +w cos Pp] k 

* EA - C“ cos Pp ~ w sin pp] i + v j + [ u sin P p + w cos P p ] k (4.39) 

The velocity due to inflow is given by: 

V /f = pQR cos if/i- pQR sin \f/j - QRA k (4.40) 

With higher order terms included, the velocity expressions are as follows: 

y EA = - V v > sin Pp + ( v ,^ x 0 +") - v ) cos p p + v^ej + v) 

+ “ + w ,x*'$x + Q [(<^ - w) sin +{jc 0 +u + v{w iX <fr 4- V jf )} cos p p 

- u(w x 4> +V*) + 4>w + v +e{]e' y + fi[{v + w{<f> Tv^m/*)} sin P p 

+ {w x v - +w) - (f)(v x v + u + Jf 0 )} cos P p 

~ w ,x v ,A v - e \)~ w ,x u + * “*i) + w]e’ z (4.41) 

V A = klR[ - v x p sin ip + ( cos p p - w ^ sin p p )p cos if/ 

- ( sin P p + w^. cos P p W\e' x + D^?[ - p. sin \J/ -{(f) sin p p 

+ (w x 4> +v^) cos Pp)fi cos t f/ +{(w x (f> Tv^) sin p p - 0 cos p p }i \]e' y 
+ nRl(w xV>x + 4>)p sin if/ + ( - sin p p + (^v, -w^)cos P p )n cos if/ 

+ (( w s ~ <f> v pc) sin Pp - c°s P p }*]e’ z (4.42) 
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From these velocities the total velocity of the clastic axis relative to the air can 
be found from: 

U = \ EA - V A = U' x e' x + Ufa + U± (4.43) 

The velocities in the single primed system arc: 

U' z = - [w^v^ + (f>]Rn sin i p 

+ [ sin p p + cos PAw x - <pv x )]R^ cos \p 

+ [ sin p p { — vv x + 4> v x ) + cos fip\ R X 

T sin fi p {w x v x w + (f>w + v) 

+ cos p p ( - w x v x u - vT x V iX x 0 + vv x v - 4>v x v -4>u - 4>x 0 ) 

— w x v J — vv^v^-Cj — W x u + 4>v x u — (f>v — (f>e [ + vv (4.44) 

Uy = Rpi sin ip + [0 sin p p + cos P^w x 4> + v^)]/?^ cos ^ 

+ [ - sin P^w x 4> + v x ) + 4> cos pp\RX 
+ sin p^(f>v — vv) + cos p p (w x (pv + v x v + u + T 0 ) 

— vv x <fiu + (pw — v x u + v + ej (4.45) 

Since apparent mass terms are retained in the expressions for the loads, the 
derivatives of the velocities with respect to ip are also needed. These are given 
by the following expressions: 
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g * 

^'2 = t — s ‘ n + cos ~ w tX + — (^x* 7 ,* + v ^.vv^ + 4>y\R^ sin 0 

* * * 

+ [ cos [S p ( - 4>v - v x (f> + w x ) - {w >x v x + 4>Y\Rn cos 0 

+ C s >n P^4>v + v x 4> — 

+ sin P f £w pc v pc w + w x v x w + 0vv + v^iv^iv + 0iv + v) 

+ cos P p { - W x v x u - w j x u - w x V jX x 0 + w j - <j>V x V 

- 4)V <X V - 4>u - v x 4>v - v x W tX u - v^o -4>u- 4>x 0 + vv x v) 

~ H VV V _ w ,x“ ~ >vv 7 ~ - 4>*v + 4>v x u 

1 - J,- - — 1 - — ~ ■t* *, * * _ _* * 

+ V X (}>U — v^w^v - v^w^e, + vv — <f>v - <j>e j — w x u 

♦ 

^5 =[-</> sin p p - cos P^w x <p + v^Rfi sin 0 

* _ * * * _ 

+ [0 sin p p + cos ^(h'> + 0vv )JC + v x ) + l]/ty cos 0 

__ * * * * 

+ [ - sin fipiw^ + 4>w x + v x ) + 0 cos f}p]RA 

* * * * * * 

+ sin ftp(<f>v + 0v — vv) + cos P^w x (f>v + w^v + <t>w x v + v^v + v^v -f u 

— J.*— — A— , 1 — J— — — 1? * ** • * 

— w x (pu — W x (pu + 0 w — (f>W x U — V X U + 0 W + V — V x u 


(4.46) 


(4.47) 
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Given the velocities in the single primed system, the velocities in the double 
primed system can be found by using Eqn. 4.13. By diferentiating Eqn. 4.13 
the following expressions for the ip derivatives of the velocities in the double 
primed coordinates arc found. 

Uy = - 4>Uy sin <f) + U’y cos <fr ~ 4>U z cos <j> — L'' sin (f> (4.48) 

U” z = 4>Uy cos <fr + Uy sin (f> — (pU z sin 4> + U z cos (f> (4.49) 


In actual implementation, the axial displacement terms u , w, and u arc ne- 
glected in calculating the aerodynamic loads. This is because with the implicit 
nature of the aerodynamic formulation, there arc no provisions for the calcu- 
lation of axial displacements in the aerodynamic section of the program. In 
fact these are quite small terms which should have little effect on the aero- 
dynamic loads. 


4 . 3.3 Explicit Implementation 

To determine explicit expressions for the aerodynamic loads, an ordering 

scheme is adopted by which higher order terms are systematically neglected. 

This ordering scheme is consistent with that used in deriving the equations of 

motion given in Eqns. 4.3 through 4.6 and include additional assumptions on 
** 

the order of 9 G which lead to the retention of some terms which may be im- 
portant when high frequency excitations such as HHC are present. Again, 
following Eqn. 4.2, terms having an order of magnitude comparable to the 
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square of the blade slopes, i.e. 0(e 2 ), are neglected when compared to terms of 
order one. 

The orders of magnitude assigned to various parameters in this study are 
given below: 


, , ** 
0(£- /2 ): d G 


0(1): siniA, cosiA, n , x o = -jr, R = , f=- = ( ) 


7 aj 0 ' 


d_ 


0(1): ~ — = (except when applied to 6 G ) 


0(e 2 ) : d G , sin d G , d G 


v ,x = 


dv 

dXr 


X A 

u 

b 

_ *1 

— 5 IV 

bR 

» b = 

T' 

£ ' 

w „ = , 

dx 0 

w 

x l 

k m 2 

k m 3 


£ ’ 

bR ’ 

^ ’ 




0(e 3 ' 2 ): 


/v^ 2 \ - U — x f 

0(0. « = — . x, = — 


The main difference from Ref. 8 is in the magnitude assigned to d G . In this 
study HHC inputs are included in 0 G as shown in Eqn. 4.1. Because HHC is 
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input at high frequencies, the angular accelerations due to it can be quite large 
and this is reflected in a larger assigned magnitude for d G than would be oth- 
erwise expected. 

The forces and moments in the undeformed coned blade coordinate system 
arc obtained by algebraically substituting the velocities given by Eqns. 4.44 
through 4.47 in Eqn. 4.13 and substituting this result in the expressions for lift 
and moment given by Eqns. 4.16 through 4.23. Using the ordering scheme 
leads to the neglect of numerous higher order terms. With the use of a small 
angle assumption, Eqns. 4.29 through 4.33 reduce to: 

K, = L c + L„ c ( 4 . 50 ) 

P" r „ = - U-1H + i-2) - Lje c + 4,)- -~<u;) 2 (4.51 ) 

(4.52) 

The simplified algebraic expressions for lift and moment are substituted in 
Eqns. 4.50 through 4.52. Subsequently these expressions are further algebra- 
ically substituted in the expressions for forces and moments in Eqn. 4.34. 
Carrying out the various algebraic manipulations and applying the ordering 
scheme produces the following final expressions for the forces and moments: 

= « fv + vVj + F,( 4- ilf.Wc + F 4 6 g - F 2 ) 

- F x l > + (2f,F 3 0 c - F 2 F 3 )v , - F,F 3 iiF, + F 2 { 4> - F x vw x 
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+ (^f - F b v ,x w ,x +2F \ F 3 v ,x ( f > + ( 2F \ d G ~ F 2> ~ F \'*' + F \ F J> 
+ 2/ 7 ,0v - F 3 wJ - F^v x w - vv} 4 ^~{0 ( ;Fi -*0 C F S 


+ Pp F 6 + ®0 F 3 + ^ + F 2>^ ~ F 2 w r x ~ w + F \<t>} 


(4.53) 


ya 


— //{/*2 + w + 


* 

vv 


,* u u + vu; ,JC 


+ Vp 


p + v l w ,x F 3 


4 vv 




*3> 


- y{(^l 0(7 + F 4®G - F t) F 2 + Pf$G F 1 ~ 2Fi)v 4 0 G F 3 F 2 V X 


+ F i( F l^G + ^G -2^2)^ + + (¥c ~ 2F 2 - 2 Pp F l)™,x 

+ (®G F 2 + F l F 2^,x W pc ~ F l™]c + W/ + 


+ ( / 4 i U 


— 2v)F 2 w^ x 4 (f|V 4 F 2 v ,x^ w ,x ( t > + F 2^g^ 


+ l^G + ^ 4^G — 2^ 2) w + F 2 F 4 ( t > + ( F i @G W x + F 2 ( t > ) v 


4 (Fy<j> 4 F-fi G v x -2p p v -2F i w x )w 4- F 3 F 4 4 F^<f>w x v 


+ (F^ 4- F, 


j* 2 * * * 

,v w — 2vvv x )w 4 d G vw — vv + F 4 w(p 4 $vvv} 


- yi^^G “ F 5 d G + ^3 0 G) + ^G@p F 6 + M^G^i “ ^G^V 


4 0 g F 6 w x 4 (0 g F j - F 5 d G +2F ~fi G 4 P p F 6 )(f> 4 F^vv^. 
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+ F 3 <K 4- d G 6 G v + 6 G F '$v x - F$ a w x + 6 g F ,</> 4- F { cf)(f) 

c 

- Fyt> w x — 0 G W 4- 6 2 G V -<0}~ — ^-{F\ +2 F,v +2F 1 F ?> V X } (4.54) 

a i 

«xa = + V*3 + f'll + ^-{WlOG + F^G - F 2 ) 

- FrfpV + (2 F l F ?) e G - F 2 F 3 )v^ - F l F i w tX + F{4> ~ F { vw x 

+ {F\ - F 2 )v x w x +2F l F i V' X 4) + (2F l d G - F 2 )v - Fjiv 4- F|F 4 0 
+ 2F,(/>v - F 2 w J ~F 2 v x w - vw} + y {(F 3 F 5 6 G + F s F(fi p 

~ F\ F# G ~ bRF*e G )-{F 2 F*8 G + F 5 F 6 0 G )v >JC + F 5 F 6 w x 
+ F 2 F 5 4> 4- F 3 F 5 iF c Hf c - F 5 F 6 v x 4> - (F 4 6 g 4- Ftfi p )v - F { F 4 4> 

- (F 4 4> + F 5 w x )v — F 2 F 4 V x (j) + (F 1 F 5 W X 4- F 3 F 5 4> 4- F 3 F 5 d G )v x 

+ (F^F^v x — F^vF^F^w x — F$w 4- F$ G w — bRF-j<f) 4- F 5 0v } (4.55) 

where the quantities F, through F 7 are defined by: 

F ] — e x 4- x 0 4- Rn sin ip F 5 = bR{x A — Vi) 

F 2 = XR 4- nRfip cos \p F 6 = Rfi simp 

F 3 = }xR cos ip F 7 = bR(V$ 4- xp — 3^) 

F 4 = M?(1 


55 



4.3.4 Solving Numerically for the Augmented States 


The differential equation governing the augmented states, Eqn. 4.21, can 
be rewritten as: 



(4.56) 


where: 

U y0 2 U lo 

A 2l = - 0 . 01365 (-^ r -) 2 , A 22 = - 0 . 3455 (- t ^) ( 4 . 57 ) 

hR bR 

This equation is forced by the V* chord down-wash velocity Q. In quasili- 
ncarization a time history for Q from the previous iteration is known and this 
can be used to determine the augmented states for the current iteration. 


A Fourier analysis may be applied to the time history of Q from the last 
iteration of quasilinearization, then Q may be expressed as: 

NH 

Q = Qo+ ^J-Qsn sin # + Q C n cos #3 (4.58) 

n= 1 

where NH is the number of harmonics retained in the Fourier analysis. The 
quantities X x and X 2 may then be solved for in terms of their Fourier coeffi- 
cients. 

NH 

*1 = *10 + X sin + *lcn cos ] (4.59) 

t 1=1 
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(4.60) 


Ml 

X 2 = X 20 + Yjthsn sin nt P + hen cos #] 

n= 1 

Substituting Eqns. 4.58, 4.59, and 4.60 in Eqn. 4.56 and equating coefficients 
of sine and cosine leads to the following solution for the Fourier coefficients 
of the augmented states A', and X 2 : 

*10 = - * ho = 0 

A 2l 

T7 -i A 2l ±^}Qsn^h2Qcn 

1 sn j ~> 7 7 

(^21 + n ) + n ^22 

y n _hAnhhl + t l^cn 

1 cn 2 7 

(^21 + ^ ^ ^ 22 

hsn = " nX \ cn * hen = nX \sn ( 4 ' 61 ) 

If Eqn. 4.57 is substituted in Eqn. 4.22, an expression for H in terms of 
A 2] and A 22 can be found and is given below: 

H = - \*2 1*1 - (■*•62) 

The quantity H is also assumed to have a Fourier series representation. 

NH 

H=H 0 +Y j tH s „ sin mj/ + H cn cos mjj~\ (4.63) 

/ 7=1 

If Eqns. 4.59, 4.60, and 4.61 are substituted in Eqn. 4.62 and this is set equal 
to Eqn. 4.63, the following Fourier coefficients for H are found. 
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(4.64) 


B | — By 

H cn = t-Q-]Qsn + \.-g^]Qcn 

B-, _ -B, _ 

^ s n ~ ^ ] Qsn t IQcn 

where the following abreviations have been used: 
B x = - 1294 A 2 \A 22 n + 2162 A 22 n 3 

^2 = 345542i(^2l ) + 2 1 62/1 22 ^^ 

#3 = 69IO[(^2i +^~) -¥n^ A 2j\ 


(4.65) 


(4.66) 


For a given point along the blade, the augmented states are driven by the 
time varying 3 4 chord downwash at that point. To perform the Fourier anal- 
ysis described above and solve for H every time the aerodynamic loads are re- 
quired would use excessive amounts of computer time. However the Fourier 
coefficients of H are independent of the azimuth and thus have to be deter- 
mined only once for each radial coordinate where the airloads are required. 
When implementing this procedure, the airloads are only solved for at a finite 
number of radial stations where the airloads are needed for use in Gaussian 
integration. Each time an airload at a given radius is needed, computer me- 
mory is searched to see if the Fourier coefficients for H at that radius have 
already been computed. If this data is available it is read from memory. If it 
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is not available it is calculated and stored for future use. Typically for a five 
element representation of the blade and eight point Gaussian integration this 
will involve the storing of forty sets of Fourier coefficients for H. When H is 
needed for a particular value of i/', the Fourier coefficients from Eqns. 4.64 
through 4.66 arc substituted in Eqn. 4.63 to give the needed value of H at 
that azimuth. 


59 



Chapter V 

CALCULATION OF HUB FORCES AND MOMENTS 
5.1 INTRODUCTION 

Hub forces and moments arc calculated using a direct force integration 
technique. A blade response in the form of a Fourier scries representation of 
the generalized coordinates is assumed to be known and the hub loads arc 
calculated based on this. The procedure used consists of three loops, one over 
the blade azimuth, one over the four blades of the rotor, and one over the 
blade radius. A description of the sequence of calculations is provided below: 

1. For a given blade azimuth, the local aerodynamic and inertia loads are 
calculated at a number of radial stations along the blade. 

2. The forces and moments at the blade root due to the local loads are then 
calculated and are numerically integrated along the blade to give the 
root loads at a given azimuth. 

3. The root blade loads due to the four blades are combined to give the 
total rotor hub loads for a given azimuth. 

4. A Fourier analysis of the total rotor hub loads around the blade azi- 
muth is carried out. 

From the resultant Fourier series, characteristics of the rotor hub loads may 
be reconstructed. A flow chart showing this procedure is given in Fig. 5.1 and 
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a detailed mathematical treatment of these calculations is provided in this 
chapter. 

Calculation of hub loads in this study differ from that in Ref. 8 principally 
in that aerodynamic loads have been calculated in an implicit manner and the 
blade modeling is more general. Implicit load calculations imply the retention 
of various higher order terms which were previously neglected, and the elimi- 
nation of the small angle assumption. The current study allows the inclusion 
of HHC inputs in the computations of the hub loads and includes provisions 
for the modeling of a control system stiffness and articulated blades with or 
without spring restraints. 

Details of the calculations are given in the following sections. 


5.2 CALCULATION OF LOCAL LOADS 
5.2.1 Local Aerodynamic Loads 

Expressions for local aerodynamic forces and moments in the coned, unde- 
formed coordinate system, as developed in Ch. 4, are used. These forces and 
moments are: P xa , P ya , P za , q xa , q ya , and q 2a . 
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5.2.2 Local Inertia Loads 


Inertia loads per unit span of the blade are calculated using D'Alembert's 
principle and integrating over the cross section of the blade. Referring to the 
geometry of the blade shown in Figs. 2.3 through 2.4, these loads can be ex- 
pressed as: 

P/ = P x £x + ?yfy + ?zfz = ~~~ \ I P*dA x ( 5 . 1 ) 

q, = Pxfx + Pyfy + Pzk = ~TT P R p xa ^x ' (5.2) 

The acceleration of a point on the blade is given by: 


a = R + 2QxR + Qx(QxR) (5.3) 

where the position vector of a point in the blade cross section, relative to the 
axis of rotation is given by: 

R = e \i + (*0 + u $x + v *y + w *z + y$, + Z 0^z (5.4) 

the rotor angular velocity vector is given by: 

A 

& = (5.5) 

and the position vector of a point in the blade cross section relative to the 
elastic axis is given by: 


, A- 

p TO + Z Q e z 


(5.6) 
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Using the coordinate transformations given in Eqns. 2.3 through 2.4, these 
vectors may be expressed in the coned undeformed coordinate system as: 


R - [U'o +u) + c j cos P p ~y^v x + <frw x ) - z^w x - f 
+ [v + y 0 +w x v x )']e y + [w - e 1 sin p p + y 0 <p + z^z 

Rp = [ -y^ v ,x + - Z d W ,X - <}>\x)$x 

+ [>’() + Lv'o^ + Z 0~\ e z 

n = n sin p p e x 4- Q cos P p e z 

The first and second time derivatives of R are then: 

R = [« - yd v,x + <t>w,x + 4>w tX ) - yd v , x + <t> w , x ) ~ z d™,x ~ <t> V - 4> v 

- z^w x - 4>v x y\e x + [v + >0 - z d<t> + w ,x v ,x + w ,x\x) 

- zd(f> + w, x V tX i]e y + l> + y^ + + 2 0 ]e z 

R = [m — yd + (f>w x + 24 >w x + <j>Wj.) — 2j)q(v x + (pw x + 

- M v ,x + <f> W ,x) - z d™,x - <t>v ,X ~ 2( t>\x - <fr v ,x) 

-2zdw x - 0V ;X - 4>V x ) - zdw x - (f)v x )]e x 

+ [v + To - z d<t> + *x v ,x + 2v W + w ,x\x) 


(5.7) 

(5.8) 

(5.9) 


(5.10) 
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~ 2z 0 (4 + W, x V tX + w ,x v ,x) - M<t> + 


+ [v'l' + i\)(f> + 2y 0 4> +y 0 <j> + z 0 ]e z 


(5.11) 


To develop the inertia loads implicitly, Eqns. 5.7, 5.9, 5.10, and 5.11 arc 
substituted in Eqn. 5.3 and subsequently this result and Eqn. 5.8 are substi- 
tuted in Eqns. 5.1 and 5.2. The following quantities arising from integration 
over the blade cross section can be defined: 



(5.12) 

(5.13) 
(5-14) 

(5.15) 

(5.16) 

(5.17) 

(5.18) 
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Py$\) dA = ( 7 m3 - Inafic sin °G cos 6 C ( 5 - 1 9 > 

A 

py&tfiA = -I m 2 COS o (;{()}; cos e G + e G sin 0 G ) 

A 


+ !,n 3 sin 0 G< ~ 0 G sin 9 G + ^G cos 0 G) 

(5.20) 

= ( 7 m2 - ! mi) sin 9 Ci cos °G 

J A 

(5.21) 

py Q z () dA = I m 2 9 G c °s^G + ^mi 9 G s ' n ^G 

J A 

(5.22) 

= I m2 cos 0 G< - 6 G sin 0 G + #G cos #g) 


+ 7 m3 s >n d d^G cos #g + 9 g sin d G ) 

(5.23) 

P z 0^ = I m2 sin2 ^G + (m3 co ^ 9 G 
J A 

(5.24) 

pz 0 z () dA = (/ m2 - I m 3&G sin 9 G cos 9 G 

J A 

(5.25) 


pz^z^dA — I m 2 sin 6 ^ — d G sin d G + 0 G cos d G ) 

J A 

- / m3 cos d(j{d 2 G cos 0 G + d G sin 0 G ) (5.26) 
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(5.27) 


pz m> dA = -i m 2 e G s\n 2 e a - i m 3 e G cos 2 e G 

J A 

pz ( y () dA = -l m2 Sin e G (d G cos 0 G + 0 G sin 0 G ) 
J A 


I m3 cos - 0 G sin 0 G + 0 G cos 0 G ) (5.28) 

This procedure leads to implicit expressions which would be extremely 
lengthy if written out but which can be easily calculated by the numerical 
program carrying out the analysis. In order to compare with Ref. 8 and verify 
the implicit expressions, explicit inertia loads were also derived using the or- 
dering scheme detailed in Ch. 4. These expressions contain some previously 
neglected higher order terms which can be important when HHC is being ap- 
plied. These expressions are given below: 

P xi~ m l x 0 + <’i + 2v] (5.29) 

P yi =mi-v ~2 u + 2 0 p w + v 4- x r cos 0 G + xfi^ sin 0 G + 0 cos 0 G )] (5.30) 

P zi=~ m l w + 2 p p v + xfacosdf;'] ( 5 . 31 ) 

q xi = - mxj cos 0 G [w + 2p p v + p£x 0 + e,) - <0 + v0] 


- mx, sin 0 c [v - V ] - (I m2 + / mJ )[0 G + £ + 


^ A ^ 

+ ( 2,1 +c - 0) sin dfj +2vv X 0 G sin 9 Cj cos d G ] - 4i 3 [0 sin 2 0 G 

- ( I + 2v x ) sin d G cos 0 G + (w x v x +2ir x v v + w *v x +2 w x -2<pv x -20v x 

- 4> - >V,x) c °s^(7 -2ay?c Sin 0 G cos 0 G -2v X 0 G cos^0 G ] (5.32) 

V = - »«/ cos 0 G [f x (vT* + PJxq +?,))- 4>(x 0 + e, + 2v)] 

+ mx, sin 0 G [x o + c, + 2v] - (/ w2 + 7 m3 )[v J G ] 

~ hn2^^,x ~ <t>v, x + v *<£ ) cos G +2(20 - vv X )0 G sin d G cos 9 G 
+ ( V ~ V) sin d G cos e G + (»> + P p +20 - w x ) sin 2 0 G +20 G sin 2 0 G ] 

- / w3 [ cos 2 8 G (p p ~ 4>v x + +2iv^ x +2v x w x +20 

* * * * * * * 

+ 0 V* +20 + 0V X - vv x ) + (v* x - V x ) sin 0 G cos 0 G 

* * 7 * 

-2(20 - w x )0 G sin 0 G cos 0 G + 20 G cos 0^1 + v x )] (5.33) 

Qzi = “ "**/ cos 0 g[ j O + cj +2v] + + mx, sin 0 G [0(jc o + c,)] 

- f m 2^,x cos 2 0 g -20 g sin 0 G cos 0 G + sin 0 G cos 0 G (vv x -20 - p p ) 

+ «c cos 2 0 6 <2^ -20)] - (l m2 + / m3 )[iv>c] 
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- ! m}p* e G sin e G cos °G ~ sin d G cos 9 ($ ~ 2 <t> ~ Pp) +2 <pd G COS 2 ^] (5.34) 

To be consistent with the quasilincarization analysis, the axial displace- 
ments u in these expressions were neglected in the actual numerical imple- 
mentation of the shear calculations. 


5.3 CALCULATION OF ROOT LOADS 


5.3.1 Summation of Local Aerodynamic and Inertia Loads 


Aerodynamic and inertia loads are combined at each local point of appli- 
cation before any further calculations are performed. Because inertia and 
aerodynamic loads have been nondimensionalized by different factors, a pro- 
portionality constant must be used before they can be combined: 

P=P,+ TP a ; q = q, + bRVq a (5.35) 

The proportionality constant is found from the ratio of the nondimcnsionali- 
zation factors: 


a lPa bRQ 2 / 2 yl b 
2R 4 


(5.36) 


where it has been assumed that the blade semi-chord used in calculating the 
Lock number y is equal to the local blade semi-chord. 
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5.3.2 Expressions for Local Contributions to Root Loads 


The local forces and moments produce resultant forces and moments at the 
blade axis of rotation. These loads must be found in the unconed hub coordi- 
nate system. Using Eqn. 2.2 the blade root forces are found by transforming 
the local aerodynamic and inertia forces which have been calculated in the 
undeformed, coned blade coordinate system to the undeformed, unconcd co- 
ordinate system. 

P x= P x cosp p - P z s\npp 



p z = p x sin Pp + p z cos fip 


(5.37) 


The blade root moments consist of contributions due to the local moments 
transformed to the hub coordinate system and due to the local forces acting 
through a moment arm from the blade center of rotation to a given point on 
the blade. The local moments are transformed to the unconcd coordinate sys- 
tem in the same manner as the forces in Eqn. 5.37. 

D 1 _ _ 

q x = q x cos P p - q z sin fi p 
-Rl 

q y = q y 

q? 1 = q x sin P P + q z cos Pp ( 5 - 38 ) 

Moments due to the local forces can be calculated from: 
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(5.39) 


where R c is the position vector of a point on the clastic axis of the blade relative 
to the center of rotation in the unconcd coordinate system: 


R c = ( e i + *0 cos Pp ~ w sin p p ) i + vj + (x 0 sin 0 p +w cos P p ) k (5.40) 

and is the vector of local forces in the unconcd coordinate system. 


P R = pK'i + P* j + p* k 


Combining Eqns. 5.39 through 5.41 produces the following expressions for 
the total blade root moments: 

Vx = - rfa sin p p + vv cos p p ) + P*v + q£ l (5.42) 

Vy = Pfco sin P p +w cos p p ) - P*(e, + x Q cos p p -w sin p p ) + q * 1 (5.43) 

4z = “ | + *o cos Pp-w sin p p ) + q* 1 (5.44) 


These expressions are calculated based on an equilibrium blade response 
which has been calculated using the mode shapes for a given blade and so arc 
valid for either hingeless or articulated blades. 
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5.3.3 Integration of Loads Over the Blade 


The total blade root forces and moments arc calculated by integrating the 
expressions for the local blade root loads over the length of the blade. The 
integrals for these expressions are given below: 


p* 0A)=\ 

R- R 

P*{x, i P)dx 

e \ 


f R R 

q R {x,\p)dx 

J<? i 

( 5 . 45 ) 

pfo)= 

r*_o 

Py{x,4/)dx 

e \ 


[ fy(x, 4> )dx 
Je \ 

( 5 . 46 ) 


[R—v 

P?{x, <p)dx 

e, 

fz('P) = 

[ q*{x,\l>)dx 

( 5 . 47 ) 


These integrals are calculated by a spanwise integration using a 16 point 
Gauss-Legendre integration formula. 

5.4 CALCULATION OF TOTAL HUB FORCES AND MOMENTS 

5.4.1 Transformation to the Nonrotating Hub Coordinate System 

Using Eqn. 2.1, the loads in the rotating system can be converted to the 
blade root loads in the nonrotating hub fixed coordinate system. The ex- 
pressions for the various components of the hub shears and moments in the 
nonrotating coordinate system are presented below. 
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Vertical Shear: 



Lateral Shear: 


sjti'l') = PxW) sin ip + Py(\jf) cos \p 
Longitudinal Shear: 

S\h('l')= PxW) cos \p - Py(\p) sin \p 
Pitching Moment: 

Mp£tp) = qxi'P) sin tp + qy{\p) cos \p 
Rolling Moment: 

Mr lip) = £?(<A)cos <P - qy(ip)s\n ip 
Yawing Moment: 

5.4.2 Summation Over the Blades 

Since each individual blade undergoes the same periodic response around 
the azimuth, the total hub loads for a given azimuth of the reference blade may 
be found by summing the loads for a single blade calculated at azimuths ap- 
propriately phased from that of the reference blade. This procedure is carried 
out in the following manner: 


(5.49) 


(5.50) 


(5.51) 


(5.52) 


(5.53) 



(5.54) 


• S w/>A)= Yj 5 S^i) whcrc </'/=</' + 

/=1 

where NB is the number of blades. Expressions for the other hub loads, S lt , 
5 ln , M p „ M rl , and M yw are found in a similar manner. 


5.5 HARMONIC ANALYSIS OF THE HUB LOADS 

The total hub loads are functions of blade azimuth and thus are expanded 
into a Fourier series and stored in this form. For example the Fourier expan- 
sion of S V M) is: 

NHR 

s v m = j v/ o + Y ^ vicn cos + Svtsn sin ^ 5 ' 55) 

a=i 

where the number of harmonics retained in the analysis is typically NFIR = 5 
and the Fourier coefficients arc given by: 

s *«> = 

i t 2n - 

s vicn - If S V M) cos minty 
J o 

S vtsn = 4" f S ) sin mj/ dtp (5.56) 

J 0 

These coefficients are calculated using a trapezoidal integration rule. The 
constant portion of each load, required for the calculation of the trim state of 
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the helicopter and for determining thrust and power coefficients, are just the 
constant coefficients in the Fourier series for each load. 


5.6 PERFORMANCE COEFFICIENTS AND NORMALIZATION 


The constant part of the rotor thrust, power, and inplane force are conven- 
tionally expressed in nondimensional coefficient form as: 

C T= J ’ C/, = T’ Ch= ~ 7 (5-57) 

p a A(QR) P a A(QR) Pa A ^RT 

where T, P, and F ri arc average rotor thrust, power, and inplanc force, and A 
is the rotor area. These coefficients may be expressed in terms of the nondi- 
mensional loads as: 


Cr = AAs 


yi b 


Vt ' 


Cp- 


4b HI 
- Myw 


y J b 


C H = 


4 bR F 
f-^ln 

y f b 


(5.58) 


Hub shear and moment results used in this study are normalized by the 
nondimensional blade mass moment of inertia in flapping, l b . For example: 




h 


(5.59) 


where <A) has already been nondimensionalized by w 0 QV 2 . The nondi- 
mensional form of the sine and cosine components of the 4/rev. hub forces and 
moments are used in the HHC calculations. They are represented by: 

5 vo 4 * s vtc4 * Sfrs 4 ; Spc 4 ! ^lns4 > ^lnc4 
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(5.60) 


m ptsA 7 m ptc4 ' m rh 4 » m rlc\ 7 m yws4 7 m yv,cA 

Referring to Fig. 2.1, the pitching, rolling, and yawing moments are as 
shown, and the vertical, lateral, and longitudinal shears arc respectively in the 
directions of the Z nr , Y nn and X m axes. 
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Chapter VI 

HHC FOR VIBRATION REDUCTION 
6.1 PREVIOUS WORK 

In 1974 McCloud and Kretz[39] tested multicyclic control applied in the 
rotating system on a wind tunnel model of a jet-flap rotor. They investigated 
the effects of HHC on the rotor loads and blade stresses. The concept of a li- 
near, quasi-static representation of the relationship between harmonics of the 
rotor loads and harmonics of the HHC was developed. An offline weighted 
least square error technique was used to calculate the transfer matrix relating 
response harmonics to control harmonics and the notation “T” was introduced 
for this matrix. Open-loop control needed to minimize a quadratic perform- 
ance function was calculated. Further open-loop studies using this method 
were conducted on a Multicyclic Controllable Twist Rotor by McCloud and 
Wcisbrich[40] and Brown and McCloud[5]. The control needed to reduce 
both blade loads and test module accelerations was considered including the 
influence of weights in the performance index. The effect of rotor lift, propul- 
sive force and speed on open-loop control were also examined. 

In 1974 Sissingh and Donham[60] conducted wind tunnel tests on a 
hingcless rotor which had HHC introduced by oscillations of the swashplate. 
Vibratory vertical shear and hub moment due to HHC were measured and 
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from this a transfer matrix was identified off-line by a least square error tech- 
nique. The sensed vibrations and a transfer matrix calculated by direct inver- 
sion were used to calculate control to be applied to the rotor. 

In 1980 Shaw and Albion[58] measured the response to HHC of a hingcless 
rotor in a wind tunnel. The transfer matrix was identified using a Kalman 
filter and was then inverted and used with feedback of measured vibrations to 
reduce rotor loads. This controller worked very well at an advance ratio of 
/j=0.1 but was less successful at other speeds because the required HHC ex- 
ceeded an externally imposed limit of 1.5°. Shaw continued this research in 
I985[59]. Fixed gain, scheduled gain, and adaptive controllers were evaluated 
with wind tunnel tests. Response to HHC was found to be essentially linear 
up to angles of ± 3°. The most significant finding was that a fixed gain con- 
troller provided 90% multi-axis vibration suppression over a wide range of 
operating conditions from hover to 188 knots and at varying load factors. 

In 1980 Taylor et al. [62,63] conducted numerical studies using the G400 
nonlinear aeroelastic helicopter simulation[2]. A Kalman filter was used for 
on-line identification of the transfer matrix and deterministic closed-loop gains 
were found to minimize a quadratic function of the vibrations and control in- 
puts. Changes in the transfer matrix due to system nonlinearity and variation 
in the flight condition were accounted for only by the Kalman filter identifi- 
cation of the transfer matrix. The control system showed good convergence 
and significantly reduced vibrations. 
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In 1980 Hammond[30] conducted wind tunnel tests of an articualted rotor. 
The vibratory hub moments and vertical shear due to HHC inputs were mea- 
sured. Both a deterministic controller with external rate limiting and a cau- 
tious controller based on a global HHC model were considered. The transfer 
matrix and the uncontrolled vibration level were identified using a Kalman 
filter. Feedback of measured vibrations was used to minimize a quadratic 
performance function. The cautious controller was found to be less erratic 
than the deterministic controller. Vertical vibratory shears were reduced sig- 
nificantly, however hub moment reductions were smaller. In 1981 Molusis, 
Hammond, and Clinc[45] extended this investigation with an analytical de- 
velopment of six controllers, four of which were then tested in a wind tunnel. 
Vertical, longitudinal, and lateral hub shears were used as feedback to reduce 
a quadratic performance function. A cautious controller tested worked 
smoothly and reduced vertical and longitudinal vibrations but actually led to 
higher lateral vibrations at low speeds. A controller identifying only the un- 
controlled vibration level was not successful when tested in the wind tunnel. 
A dual controller incorporating probing terms designed to enhance system 
identification was derived but not tested. 

Since 1980 Wood et al. [67, 68, 25] have developed software and hardware 
for HHC implementation on an OH6-A and have conducted the first success- 
ful flight tests with active HHC. The cautious controller developed by 
Molusis[45] was used in closed-loop studies and additionally, open-loop re- 
sults showing phasing of HHC for shear minimization were presented. A 
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Kalman filter was used to identify the transfer matrix and the uncontrolled 
vibration level. This approach was successful in reducing fuselage vibrations 
at low to moderate speeds and during moderate transient maneuvers. 

In 1982 Johnson[36] wrote a comprehensive review including the analytical 
development and evaluation of essentially all the control algorithms used to 
that date. The discussion of the algorithms was based on single-input single- 
output system models and did not include actual helicopter simulations. 

In 1983 Molusis[46] and Molusis, Mookerjee, and Bar-Shalom[47] per- 
formed analytical studies of the effects of nonlinearity on controller perform- 
ance. In Ref. 46 performance of a deterministic controller and a cautious 
controller were compared using a simple, arbitrary nonlinear simulation model 
not based on an actual helicopter analysis. Reference 47 used a nonlinear 
Voltcrra scries analytical system model with parameters derived from a non- 
linear aeroelastic vibration analysis using the G400 computer simulation. This 
model was used to consider the effects of nonlinearity on controller perform- 
ance and Kalman filter stability. 

Ham[28,29] and McKillip[41 ,42] have developed the concept of Individual 
Blade Control (IBC). The pitch of each blade is controlled individually, either 
through swashplate actuators for three or less blades, or by individual actua- 
tors in the rotating system for more than three blades. Blade mounted accel- 
erometers are used to sense blade motion and feed this information back to the 
controller. These references include analytical studies using simple models not 
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based on actual helicopter aeroelastic simulations, and limited wind tunnel 
testing. Results indicate that IBC can be used to reduce helicopter vibrations 
using relatively simple controllers at the expense of a more conplicatcd physical 
implementation due to components in the rotating system. 

In 1983 Chopra and McCloud[10] developed a linear, quasi-static, fre- 
quency-domain model for the response of a rotor to HHC. This model related 
six harmonics of vibration to six harmonics of HHC input through use of a 
transfer function based on experimental data. This model was used to evaluate 
the performance of a deterministic controller for various flight speeds and in 
the presence of simulated measurement noise and inaccurate initial estimates 
of model parameters. 

In 1983 Jacob and Lehmann[34] presented an analytical study of HHC. 
Two simulations were used, one a rigid spring restrained blade, and one a fully 
clastic blade with one mode modeled in each of flap, lag, and torsion. Simu- 
lation was basically of an open-loop type but a static search algorithm was 
used to find the HHC necessary to minimize a performance criterion. In 1985 
Lehmann[37] continued this work with the wind tunnel testing of a four 
bladcd hingeless rotor. The effect of open-loop HHC consisting of 3, 4, and 
5/rev pitch changes in the rotating system was investigated. In 1988 Lehmann 
and Kube[38] successfully reduced 4/'rev. hub vibrations using closed-loop 
HHC applied to the same rotor. A frequency domain, minimum variance 
controller based on a local HHC model and using Kalman filter identification 
of the transfer matrix was used. 
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Davis[ll] in 1984 produced a computer simulation using the G400 code. 
Deterministic, cautious, and dual controllers using both local and global HHC 
models were investigated. These controllers employed Kalman filter identifi- 
cation of the vibrations and/or transfer matrix. A general controller which can 
be specialized to any of these eases by changing certain parameters was pre- 
sented. Vibration reductions of 75 to 95 percent were achieved with each of 
the controllers when they were properly tuned. Control amplitudes of less than 
one degree were needed for these reductions. 

In 1986 Miao, Kottapalli, and Fryc[43] and Walsh[64] presented results 
of HHC flight tests using an S-76A helicopter. Only open-loop HHC was im- 
plemented. Substantial vibration reductions were achieved but at high speeds 
these reductions were constrained by hardware limitations which w r ould not 
allow high enough HHC amplitudes. Various combinations of HHC inputs 
were used to reduce combinations of several vibratory quantities. 

In 1986 Polychroniadis and Achache[53] presented flight test results for 
closed-loop HHC applied to an SA 349 Gazelle. One deterministic and two 
stochastic controllers, each based on a global HHC model, were tested. All 
three controllers gave good vibration reductions of up to 90% in steady state 
flight and differed mainly in their self-adaptive performance. Since HHC 
amplitudes in this study were limited to one degree, it was speculated that by 
allowing higher values of HHC greater reductions in vibratory loads could be 
obtained. 
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In 1986 Jacklin[32,33] presented an analysis of system identification tech- 
niques for HHC. Least squares and Kalman filter type identification tech- 
niques were reviewed and a new mean least squares algorithm was developed. 
Various aspects of these algorithms were compared. The mean least square 
algorithm uses equations similar to those of a Kalman filter however these 
equations turned out to be simpler. It was found that this algorithm converged 
to the correct estimate only when a signal averaging method over several iter- 
ations was employed. 

In 1986 Hanagud et al.[31] produced an analytical, coupled rotor-airframe 
analysis of HHC using a one degree of freedom, flap only, rotor model. A 
global HHC model was assumed and a deterministic control law was used. 
The transfer matrix was calculated by a method which required the matrix to 
be consistent with measured vibration data while at the same time being as 
close as possible to some a priori estimate of the transfer matrix. This ap- 
proach led to reductions in vibration levels. 

Gupta[26], Gupta and DuVal[27], and DuVal, Gregory, and Gupta[15] 
have developed an extension of linear-quadratic-Gaussian (LQG) design 
methods using frequency-shaped cost functionals. A vibration controller was 
obtained by minimizing a cost functional which places a large penalty on fu- 
selage accelerations at set vibration frequencies. The optimal control solution 
involved feedback of fuselage accelerations through undamped oscillators, 
tuned to the frequency at which vibrations are to be suppressed. This ap- 
proach has the advantage that on-line harmonic analysis of the vibrations is 
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not required and the resultant controller is simple to implement because it is 
a constant gain regulator with filters in the feedback loops. A dynamic model 
of the rotor-fuselage combination is needed in the LQG design procedure. 
Since this model will change with flight condition, gain scheduling is required 
to account for these different conditions. This control law was implemented 
using a blade-element simulation of the Rotor Systems Research Aircraft 
(RSRA). Accelerations were reduced by at least 80% in all channels except 
the vertical, in which the initial vibration level was two orders of magnitude 
below that of the other channels. 

In the appendices to Murphy's 1987 report[48], Motyka et al. developed 
an approach to HHC implementation using a Linear Quadratic Guassian with 
Loop Transfer Recovery (LQG/'LTR) design methodology. Baseline control- 
lers for various operating conditions were designed using LQG/LTR to ensure 
robustness. These model based compensators contained elements of a linear 
system model of the helicopter and so had to be implemented with a gain 
scheduling or table look-up approach for implementaion at different flight 
conditions. Attempts were made to enhance these controllers using various 
adaptive control techniques. Numerical simulations of the control algorithms 
were carried out on an elastic model of a helicopter. The basic model based 
compensators were very effective in reducing vibrations with vibration re- 
ductions of 99% in vertical, 94% in lateral, and 68% in longitudinal directions. 
Various adaptive enhancement approaches did not improve control perform- 
ance. 
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6.2 THE HHC MODELS 


6.2.1 Linearity of HHC 

It is generally assumed that the relationship between the vector of HHC 
input harmonics 6 and the vector of vibration output harmonics Z can be re- 
presented by a linear, quasi-static, frequency-domain model. A discrete time 
model of the helicopter is assumed such that the sampling time step At is suf- 
ficient for all transients to die out and for the vibration harmonics to be mea- 
sured. This time step would typically be one rotor revolution. The actual 
helicopter being modeled, and the mechanisms being used to control it, are of 
course complicated continuous time systems. A schematic of such a helicopter 
model is shown in Fig. 6.1. This figure shows a helicopter, with given flight 
conditions, which produces some time history of vibrations. A harmonic 
analysis of these vibrations is carried out and a discrete set of vibration har- 
monics is fed to the controller. These harmonics are used, typically in a Kal- 
man filter, to identify problem parameters. From these parameters and stored 
information from previous samplings, gain calculations are made and har- 
monics of an HHC feedback to minimize vibrations are calculated. From these 
harmonics a time history of HHC inputs is calculated and this is fed to elec- 
tro-hydraulic servo-actuators in the helicopter which set the actual blade pitch 
angles. These HHC inputs modify the helicopter dynamics producing new vi- 
bration levels which are used to begin another iteration of the process. 
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The baseline uncontrolled vibration level of the helicopter, Z 0 , is a nonlinear 
function of the helicopter operating condition and is dependent on nonlinear 
phenomena associated with aerodynamic, structural, and inertial aspects of the 
helicopter dynamics. Linearity is being assumed only in the response to HHC 
inputs. Previous experimental studies[30,59] have indicated that only 0°-3° 
of HHC inputs are needed to alleviate vibration, so the assumption of linearity 
is thought to be reasonable. Still, some studics[46,47] have found evidence of 
nonlinear response to HHC over this amplitude range, in which case the locally 
linearized model explained in the next section may be needed. 


6.2.2 Global Model 

The globalmodel of the helicopter response to HHC assumes linearity over 
the entire range of control application: 

Z(/ + l)=Z 0 + T6(i) (6.1) 

The vibration vector Z at time step / + 1 is equal to the baseline vibration Z 0 
plus the product of the transfer matrix T and the HHC vector 6 at time step 
i. This implies that T, the transfer matrix relating HHC input harmonics to 
vibration output harmonics, is independent of 9(i). This concept of a linear 
quasi-static representation of the relationship between harmonics of helicopter 
response and harmonics of HHC was introduced in Ref. 39 where the notation 
T for this matrix was first used. 
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6.2.3 Local Model 


The local model of the helicopter response to HHC is based on a lineari- 
zation of the response about the response to the current value of the control 
vector: 

Z(i + 1 ) = Z(i) + T [0(i + 1 ) - 0(/)] (6.2) 

The vibration vector Z at time step i + 1 is equal to the vibration vector at 
time step i plus the product of the transfer matrix and the difference in the 
control vector from time step i to time step / + 1. The transfer matrix T now 
relates changes in the HHC input harmonics to changes in the vibration output 
harmonics. This allows for variation of the transfer matrix T with input 0(/). 
Alternately this may be written as: 

AZ(/ + 1) = TAd(i + 1) (6.3) 


6.3 IDENTIFICATION 
6.3.1 A ssumptions 

In applying HHC algorithms for vibration reduction, it is assumed that the 
HHC inputs 0(t) are known without error. Measurements^/) of the hub vi- 
bratory loads Z(/) are made and are assumed to include zero mean white 
Gaussian measurement noise V(i): 
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y<i)=Z(i)+ V{i) 


(6.4) 


Based on the measurements^/), different parameters may be identified for the 
following cascs[36]: 

1. Local Model, Identify transfer matrix T 

2. Global Model, Identify baseline vibrations Z 0 , T assumed known 

3. Global Model, Identify T, Z 0 assumed known 

4. Global Model, Identify T and Z 0 

If the transfer matrix T or the baseline vibration vector Z 0 is to be identified, 
it is assumed to vary with flight condition and is represented by a random walk 
model[45]: 


7\i + 1 ) = 710+ W{Q 

(6.5) 

Z^i+ i)=Zo(/)+ iv 0 (0 

(6.6) 


where W(i) and W Q {i) are zero mean white Gaussian process noise. 


6.3.2 Kalman Filter Identification 
6.3.2. 1 General Kalman Filter 

A discrete Kalman filter[6] may be used to estimate a discrete random 
process x which can be modeled in the form: 

x k+\ = < Pk x k+ w k ( 6 - 7 ) 

Observation of the process occurs at discrete points in time according to: 
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z k = H k x k + v k 


( 6 . 8 ) 


where: 

*k = (n x 1 ) process state vector at time t k 

4>k =(nx n ) matrix relating x k to * fc+I 

in the absence of a forcing function 
w k = (« x 1) vector - white (uncorrelated) sequence 

with known covariance structure 
=(mx 1 ) vector of measurements at time t k 
Hk — (m x n) matrix giving the ideal (noiseless) connection between 

the measurement and the state vector at time t k 
Vk = (m x 1 ) vector of measurement errors 

- white sequence with known covariance structure 
and uncorrelated with the w k sequence 
The covariance matrices for the w k and v k vectors are: 

/;** 

^={0* /;** 

T 

E\_w k v i ] = 0 for all / and k (6.11) 

The Kalman filter consists of the following algorithm to produce an updated 
estimate of the state vector x with the assumption of some prior estimate x~. 
Update equation: 


(6.9) 

( 6 . 10 ) 
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*k = *k + K k( z k ~ H k x k ) 


( 6 . 12 ) 


Kalman gain equation: 

K k =P k H T k (H k P k H T k +R k V' (6.13) 

Error variance matrix equation: 

P k = (I-K k H t )P k (6.14) 


State vector projection equation: 

A— , A 

x k + 1 - V/Sk 


(6.15) 


Error variance matrix projection equation: 

P k + 1 = 0 k P k^k + Qk ( 6 - 1 


This algorithm is used to identify the transfer matrix T and/or the baseline 
vibrations Z 0 . From the assumptions of Eqns. 6.5 and 6.6 the matrix <f> k will 
always be the identity matrix I. Therefore Eqn. 6.15 becomes: 


A— A 

*A+l = *k 


(6.17) 


and Eqn. 6.14 and 6.16 can be combined to give: 


P k + 1 = P k ~ K k H k P k + Qk 


(6.18) 


If there are no process dynamics, i.e. Q k = 0, then the Kalman filter is 


equivalent to a recursive weighted least-squares algorithm in which the 
weighting on the measurement error is equal to . 
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6.3.2.2 Kalman Filter for Case 1 

For the local HHC model only T is identified. When Eqn. 6.5 is considered 
one row at a time, one has for the yth row: 

Tp+ 1)= Tp)+wjtQ (6.19) 

From Eqns. 6.4 and 6.3 an observation equation for one row of the T matrix, 
equivalent to Eqn. 6.8, can be written as: 

Ayp)=Ad T (i)T J T (i)+v(i) . (6.20) 

where Ay/;) is the yth componemt of A y(i) and Tp) is the yth row of 71 (/'). 

By comparing Eqns. 6.20 and 6.5 with Eqns. 6.7 and 6.8, the Kalman filter 
equations for one row of the T matrix can be obtained as: 

Tj{i + 1)= t[(i) + K(i + 1 )[A yp + 1 ) - A6 T {i)fJ{i)-} (6.2 1 ) 

K{i + 1 ) = Ppmi)lA6 T (i)Pp)AdU) + RJLOT 1 ( 6 . 22 ) 

Pp + 1 ) = Pp) - K(i + 1 )A6 T (i)Pp) +Qp) (6.23) 

If it is assumed that the ratio of parameter and measurement noise vari- 
ances, Qp)IRp ) is the same for every measurement, then Eqn. 6.21 can be 
written[36] in matrix form as: 

T(i + 1) = T(i) + [Ay(i) - T(i)Ad(i)-]K T (i+ 1) (6.24) 
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and K(i 4- 1) and P(i + 1) only have to be calculated once. There is no reason 
to expect this assumption to be true, but a great reduction in computation is 
gained by its acccptancc[36]. 

6.3. 2. 3 Kalman Filter for Case 2 

Using the global HHC model, identifying only the T matrix, and assuming 
Zq/ is known, requires only the substitution in Eqn. 6.21 of yp +1)— Z 0j for 
A yp + 1) and 6{i) for A0(/), thus one obtains: 

T]\i + 1 ) - Tj{i) 4- K(i + 1 )[yp 4- 1 ) -Zqj - 0 T {i)fJ{i)] (6.25) 

K(i+ 1)= Pp)d T (i)ie T (i)Ppm+ Rp)T l (6-26) 

Pp + 1 ) = Pp) ~ K{i + 1 )d T (i)Pp) +Qp) (6.27) 

and the matrix form of these equations may be written similarly to Eqn. 6.24. 

6.3.2. 4 Kalman Filter for Case 3 

Using the global model, identifying only Z 0 , and assuming the T matrix is 
known, requires an observer for Z 0 in the form of Eqn. 6.8. From Eqns. 6.1 
and 6.4 this can be written as: 

yp) - e T {i)rJ{i) = z 0 p) + *0 (6.28) 
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Comparing Eqns. 6.6 and 6.28 with Eqns. 6.8 and 6,7 the Kalman filter for 
this case can be seen to be: 


% + 1 ) = A)/') + K (> + I )k’jd+ 1 ) - e\i)Tj{i) - Z 0y (0] (6.29) 

K(i+ D= Pp)lPfi)+ R k T l (6.30) 

Pp+ 1 ) = Pp) — K{i + \)Pp)+Q k (6.31) 

where now all quantities arc scalars. In vector form Eqn. 6.29 can be written 
as: 

Zp + 1 ) = zp) + 0(1 + 1 ) - 7 \i)d(i) - Zp)~\K{i + 1 ) ( 6 . 32 ) 


6.3.2.5 Kalman Filter for Case 4 


Using the global HHC model and identifying both Z 0 and T requires that 
elements of Z 0 and T be incorporated in a single vector. From Eqns. 6.1 and 
6.4 an observation equation for an element of Zq and a row of the T matrix 
may be written as: 


yJL0=U 


d T (i)l 


Zq P) 

rj{i) 


+ K0 


(6.33) 


or by defining a new vector and matrix as: 


0(0 = 



C(0 = [zp) 7X0] 


(6.34) 
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this becomes: 


yfi) = [ 2 r ( 0 ][ c /( 0 ] + ><0 


(6.35) 


Comparing Eqns. 6.6, 6.5, and 6.33 with Eqns. 6.7 and 6.8 the Kalman fil- 
ter for this case can be seen to be: 


Cj{i + l)= tj(i) + K(i + 1 )[yp + 1 ) - Zq/0 - (6.36) 

K(i+ 1 ) = Ppft T (i)lO T (i)PpMi) + Rp)T 1 ( 6 - 37 ) 

pp +\)=pp)- Ku + 1 n T mp) + Qp) < 6 - 38 ) 

where now Pp) is the covariance matrix for Cj{i) and is composed in the fol- 
lowing manner: 


Pp) 


P \\j 75 12/ 

P 2 1/ P 22j 


(6.39) 


where P Uj is the covariance of Z 0y , P 22j is the covariance of theyth row of the 
T matrix, and P nj is the cross covariance of Z 0y and theyth row of the T ma- 
trix. 


Again Eqn. 6.36 may be written in matrix form as: 

C(i + 1)= C(0 + [y{i + I) - Zp) - T(i)d(i)]K T (i + 1) (6.40) 

In summary the inputs to the Kalman filter given in Eqns. 6.12 through 6.18 
for the four cases are given in table 6.1. 
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TABLE 6.1 


Kalman Filter Inputs 


Case 

X 

z 

H 

1 

Tj 

A yj 

A0 

2 

Tj 

yj ~ z v 

0 

3 

Z 0J 

yj-Tp 

1 

4 

1 1 

1 1 

yj 



6.3.3 Least Mean Squares ( LMS) Identification 

Jacklin[32,33] has developed an extension of the LMS algorithm of Widrow 
and Hoff[65,66] to handle multi-input, multi-output helicopter parameter 
identification problems. He uses this approach with “inverse control,” i.e. 
application of the inverse of the local transfer matrix as a controller to reduce 
vibrations. Denoting the inverse of T by C , the vibration control commands 
are produced by: 

*A 6{i)=-Cy(i) (6.41) 

*0(i+ l)=0(/)+*A0 (6.42) 

The matrix C need not be a true inverse of T but is a general matrix which 
performs the inverse function of the T matrix. 
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The LMS algorithm forms an adaptation error vector consisting of the dif- 
ferences between the actual and estimated changes in HHC input harmonics. 


Referring to Fig. 6.2: 


£ = Ad - Ad 

= A6 - CAZ (6.43) 

= A6 - CTAd 


If C were the exact inverse of T , the error vector would be zero. A steepest 
descent method is used to update C: 


C(/+ 1 ) = C(i) - K s \ 



(6.44) 


Here K s is a gain term governing the amount of correction being made. 


The error squared term for one row of C is given by: 


«/ - ( £ yX £ y) 7 (A Oj - CjAZXAdj - A Z ' c; ) 




(6.45) 


Expanding this and differentiating with respect to Cj leads to: 


de 


dCj 


7 = 2{CjAZ- A6j)AZ T 


T 


(6.46) 


The update equation for a row of C then becomes: 


Cp + 1 ) = Cp) - 2kpCpZ - A0j)AZ 1 


(6.47) 


All the row equations may be combined in matrix form to give: 


C(i + 1 ) = C(i) + 2[A d(i) - C(/)AZ(/)][AZ ; ] 


(6.48) 
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where A Z T is a matrix whose yth row is the A Z T (i) vector multiplied by k 2 : 


k l AZ T (i) 

k 2 AZ T {i) 

kjAZ T {i) 


(6.49) 


In the investigations of Refs. 32 and 33 it was found that the values of k } 
needed to be tuned for optimal performance. 


6.4 CONTROL ALGORITHMS 

/ 

6.4.1 Minimum Varaince Control Algorithms 

The vast majority of all HHC investigations[67,68,53,45,59,37,62,10,34] to 
date have used linear optimal control solutions based on a quadratic cost 
functional. Minimum variance control is based on the minimization of a cost 
functional which is the expected value of a weighted sum of the mean squares 
of the control and vibration variables. 

Minimum variance controllers are obtained by minimization of the cost 
functional: 

J = E{z T (i)W z Z(i ) + d T (i)W e d(i) + A6 T (i)W Ae A9{ij\ ( 6 . 50 ) 

where: 
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E{) = indicates expected value, 

accounting for uncertainty in the system 

A 0(/) =0(0-00-1) 

\V 2 = diagonal weighting matrix on vibrations 

lV g = diagonal weighting on control amplitudes 

W s q = diagonal weighting matrix on rate of change 

of the control amplitudes 

Typically Z, 0, and A 0 consist of the sine and cosine components of the 
N/rcv vibrations and HHC inputs. The weightings of each of these parameters 
may be changed to make it more or less important than the other components. 

The minimum variance controllers are obtained by taking the partial de- 
rivative of J with respect to 6{i) and setting this equal to zero: 


dJ 

d6(i) 


= 0 


( 6 . 51 ) 


The resulting set of equations may be solved for the optimal HHC input *0(0- 


The form of the resulting algorithm will depend on whether the global or 
local system model is used and on whether a deterministic, cautious, or dual 
controller is desired. For the dual controller the cost functional must be mod- 
ified to include system probing terms. Detailed descriptions of these algo- 
rithms including the terminology for deterministic, cautious, and dual 
controllers are given in the following sections. 
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6.4. 1.1 Deterministic Controller 


For deterministic control, uncertainties in T and Z 0 are ignored and all 
system parameters arc assumed to be known exactly. The cost functional then 
becomes: 

J = Z T (i)W^Z(i) + 6 T {iW e 6(i) + A0 T (iW A9 Ad(i) (6.52) 

For the local model Eqn. 6.2 is substituted for Z(i) and setting the derivative 
equal to zero in Eqn. 6.51 leads to an expression for the optimal control. The 
algebra involved is quite lengthy but will be outlined for this one case. First 
substituting Eqn. 6. 2 in Eqn. 6.52 gives: 

J = [Z(/ - 1 ) + 7(0(0 - d(i - 1 ))] 7 ’lF z [Z(/ - 1 ) + 7(0(0 - 6(i - 1 ))] 

+ d T (iW e d(i)+ [0(0 -0(1- l)] r W / Ae [0(O-0(/- 1)] (6.53) 

Expanding the transposes gives: 

/= \Z T {i- \)W z T-e T {i- 1 )T T W z T-6 T (i- 1)^ A0 ]0(O 
+ d T (i)[T T W z Z(i - 1)- T T W z T6{i- 1)- W A9 d(i- 1)] 

+ d r (i)[T T W z T+ W e + W Ae nO+Z T (i- 1 WIZKi- 1 )- Td(i — 1 )) 

+ e T {i- l)[lV Ad d(i- 1)- T r w/Z(i- I)- T8(i — 1))] (6.54) 

Taking the derivative with respect to 0(0 gives: 


98 



-^L-=[Z T (1- \)W Z T — d T (i — \)T T lV z T-e T (i- 1 )W &e 2 
+ [T T W^{i- 1)- T T W Z T6(1 - \)-n A9 d(i- 1)] 

+ 2\T T W : T + W Q + W Ae -]d(i) = 0 (6.55) 

Solving for the optimal control *0(/) gives: 

‘0(0 = D[ - T r lV z Z(i) + W Ae 0{i) + T T lV z Td(i)] (6.56) 

where: 

D = [T T W Z T+ W e + W A6 ]~ l (6.57) 

Subtracting 6(i) from both sides gives: 

*A0(i + 1 ) = D[ - ^W^i) - W e d{i)\ (6.58) 

For the global model, using Eqn. 6.1, Z 0 + 70(0 may be substituted for 
Z(0 giving: 

*0(/+ l)=Z)[-7'V z Z 0 + W Ad d(i)] (6.59) 

or: 

*A0(t + !)=£>[- r 7 ’w / z Z 0 + iv d d(i) - T r W z Td{i)\ (6.60) 
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6.4. 1.2 Cautious Controller 

A cautious controller takes account for system parameter uncertainties by 
taking the expected value of the cost functional. The control 9 is assumed to 
be explicitly known so that all uncertainty is in the vibrations Z. The cost 
functional then becomes: 

J = £ WtfZji ol + 0 T UW e d(i) + Ad T (i)lV A8 Ad(i) (6.61 ) 


For the local HHC model, Eqn. 6.2 is used for Z{i) and the covariance of 
theyth row of the T matrix, Pp ) , is taken from the Kalman filter for the local 
case, giving: 




(6.62) 


The solution may be found by defining an effective W &e as: 




(6.63) 


and substituting in Eqn. 6.56. The optimal control is then: 


*0(i+ 1)= D 


T T W^.(i) +{^V a6 + T T W z T + W ZJj Pfi)J{i) ^ 


‘A d(i +!)=£>[- T T W^0 - W e 9{ij\ 


(6.64) 


(6.65) 


where: 
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-1 


( 6 . 66 ) 


D = 


rV z r+ iVq + w'** + 

j 


The constant X c has been introduced to allow modification of the amount 
of caution in the controller. For X c = 0 this reduces to the deterministic con- 
troller. This controller introduces a constraint on the rate of change of control 
similar to W Ad . This control is proportional to the uncertainty in the T matrix. 
As the covariance of the T matrix increases, the rate of change allowed by the 
cautious controller is reduced. 


For the global HHC model, Eqn. 6.1 is used for Z(i). Uncertainties in Z(i) 
arc now due to both Z 0 and T. This gives: 

E^W zjj zf{i)^ = Z r (i)H^Z(0 + e\i)(^W zjj P^Q{i) (6.67) 

where 6 and Pj are defined in Eqns. 6.34 and 6.39. The optimal control is 
found by using J as defined in Eqns. 6.61 and 6.67 in Eqn. 6.51. After al- 
gebraic manipulations similar to Eqns. 6.53 through 6.55 the optimal solution 
is: 


*0(i+ 1)= D\ 


-t t w z z 0 + 

j 


where now: 


D = 


t t w z t + U'g + + X^£w zJ1 P 2 ^i) 


oT l 


( 6 . 68 ) 


(6.69) 
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The constant X c has been added in again as in the local cautious controller. 
This controller introduces a constraint on the magnitude of the total control 
angles proportional to the covariance of the T matrix. It also adds a constant 
term to the solution which is proportional to the cross covariance P n - 


6.4. 1.3 Dual Controller 


A dual controller attempts to improve long term system identification by 
actively probing the system while at the same time maintaining good control. 
Optimal dual controllers generally are too conplex for practical implementa- 
tion. the sub-optimal dual controller presented here is taken from Ref. 23 and 
is presented in Refs. 36 and 11. In this approach, the cost functional J of Eqn. 
6.52 is modified by having a term added which is a function of the estimation 
error and acts to probe the system. The new cost functional is: 


J D 


= J~*D 


I API 
I f\i + 1)1 


(6.70) 


Here | | indicates the determinant of a matrix. The controller therefore at- 
tempts to provide good control by minimizing J while furthermore attempting 
to reduce P{i) in comparison to P{i + 1), i.e. improve identification, in order to 
reduce the second term in Eqn. 6.70. The constant X D allows tuning of the 
degree of probing by the algorithm. In Ref. 23 it is shown that: 


\f\0\ 

I + 1)1 


= X 


D\ 


1 + 


A6 T {i)P(i)Ad(i) 

R 


(local) 


(6.71) 
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(global) 


(6.72) 


e T (i)Pumi) 


Because this gives a quadratic term in J, it allows the minumum variance sol- 
ution to be easily found. 


IPO) I 




I P(/+ 1)1 


= X 


D 


For the local HHC model the cost functional is given by Eqns. 6.70 and 
6.71. The dual control solution may be found by replacing IV A6 by an effective 
W &0 in Eqns. 6.56 through 6.58. 

IV A 0, # - "'’aS-Vx- < 6 ”> 


This leads to the following solution: 


*0(/ + 1 )=D 


r T w z z(i) + 



+ T T W z T- 


^D' 



(6.74) 


where: 


D = 


W 7 T + W r + 


w 1 
W A9 ~ d~ 


Subtracting 0(/) from both sides gives: 


*A0(i +!)=£>[- 7' 7 'fT 2 Z(«) - W d 0(ij] 


(6.75) 


(6.76) 


The probing term in the local model control leads to a reduction in the 
constraint on the rate of change of the control proportional to the covariance 
of the T matrix. 
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For the global HHC model the cost functional is given by Eqns. 6.70 and 


6.72. This gives: 

J D = Z T {i)\V z Z(i) + d T (i)lV e d(i) + Ad T UW A6 Ad(i) 

T PU ) 

- 0 (6.77) 


Comparing this with the cautious controller in Eqns. 6.61, 6.62, and 6.64 
through 6.66 leads to the following dual controller: 


*6{i + 1)= D\ 


-T T W^W^i)-X 



(6.78) 


where now: 


D = [t T W 2 T + W t + W &e - 1 (6.79) 

Again P u and P 22 are as defined in Eqn. 6.39. 


The probing term in the global HHC model dual controller reduces the ef- 
fective constraint on control magnitude in proportion to the covariance of the 
T matrix. In addition a constant term is added to the control which is pro- 
portional to the crosscovariance P n . 


6.4. 1.4 General Minimum Variance Controller 

Reference 1 1 suggests that the deterministic, cautious, and dual controllers 
just outlined can be programmed into two combined algorithms, one for the 
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local and one for the global HHC model. This form depends on the assump- 
tion that the noise characteristics for each row of Z 0 and the T matrix are the 
same so that P is the same for each row. In these algorithms the appropriate 
controller can be implemented by changing the value of a parameter /?. The 
algorithms arc as follows: 

Local Model 

* A 0(i +!)=-£>[ W e d{i) + r T W^Z(i)] (6.80) 


D = 


T r lVJ + U' e + IV AS + /Up]T(*C| ' 


(6.81) 


Global Model 


*A0(/+ 1 )=- D 


t t w z t + 


W a + 


fX + T T W^ + fXP^Z W-. 

j ) j 


ZJJ 


(6.82) 


D = 


t t wj + w, + IV Ae + 

j 


In each case the value of ft is given by: 

Controller 

Deterministic 

Cautious 

Dual 


Value of P 
0 
1 

-I 

j 


(6.83) 
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6 . 4.2 LMS Direct Inverse Control 


Direct inverse control with LMS identification of the inverse matrix was 
implemented by Jacklin[32,33]. From Eqn. 6.2, if 7A0 is equal to -Z{i) then 
vibration will be eliminated. This leads to the inverse control A0 given by Eqn. 
6.41. As in Eqns. 6.41 and 6.42, the optimal change in HHC, *A0(/), is found 
by direct multiplication of the C matrix and the total vibration vector >{/) . 
The C matrix performs the inverse of the operation carried out by the T matrix 
however it docs not have to be the formal inverse, in a mathematical sense, of 
the T matrix. The inverse matrix C is identified using the LMS identification 
algorithm explained previously. In Ref. 32 it was found that even when vi- 
brations were nearly eliminated, a steady-state error in inverse matrix identifi- 
cation remained because as vibrations went to zero the update term for the 
inverse matrix also went to zero. To improve identification, a control relaxa- 
tion parameter was used allowing the implemented control to be only a frac- 
tion of commanded control: 

*Q{i + 1)= 0(0+ K C *A6 (6.84) 

The use of control relaxation was shown to produce lower steady-state iden- 
tification errors at the expense of slower vibration reduction. 
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6.5 HHC IMPLEMENTATION 


6.5.1 HHC Input Formulation 

Required HHC pitch commands arc introduced into the structural model 
as time varying changes in the blade root pitch angle. Two methods of for- 
mulating the HHC inputs will be presented. The first method expresses HHC 
components in terms of the magnitude and phase of a sine input. This is useful 
when dealing with open loop applications of HHC but docs not allow each 
degree of freedom to be varied while all others remain zero. The second 
method of formulation expresses HHC components in terms of the amplitudes 
of sine and cosine inputs. This allows the degrees of freedom to be individually 
incremented and lends itself to implementation of automatic controllers. 

First Formulation 

0//// ~ [00 sin ( 5J ///W' + ^o)3 + [0 C s > n (5W' + cos ^ 

+ [0^ sin(aJ///y0 + sin ij/ (6.85) 

where a> HH — — and 0 O , 6 C , 6 S , 4> a , (f> c , (f> s , and io HH are constant with 
respect to t/^. 

In addition, expressions for the first and second derivatives of d HH with re- 
spect to i jj are needed. They are as follows. 

* 

0//// = [0Q ™HH cos ( a> //// 1 /' + $<>)] 
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+ [_Q&>HH cos(o> j / 1 fij/ + <f> s ) — 0^>sin + </> c )] sin \}/ 

+ [flea?/,// cos (ro / //W' + 4> c ) +9 S sin(to /7// ^ + (f> s )] cos \J/ (6.86) 

* * —2 • — 

^//// = C _ "//// sin ("///W' + <£ 0 )] 

+ [2As"//// cos((o //// <A + 4- o5 2 [iif)s\n{(5ijfi^ + $ c )] cos tj/ 

+ [ ~29(<Ofjj { cos + <f> c ) — 1 + o7////) s\n(ajjf[{ij/ + 0J] sin ij/ (6.87) 

Second Formulation 

^//// = tfos s * n "///W* + 0oc cos "///WO 

+ [^CS s ^ n "///W' + #CC cos ^ // //*A D cos ^ 

+ [0 S5 sin oJ //// iA + 0 5C COS oJ^yi//] sin (6.88) 

where 0 OC , 0 O5 , 0 C5 , 0 CC , 0 55 , and 0 5C arc constant with respect to if/ . 

The First and second derivatives of 0 HH with respect to are as follows. 

* 

~ [^05^//// cos "////A ~ &0(PHH sin "///WO 

+ [{Ocs^HH +&SC) cos "///W* + ( _0 CC"//// +^5s) sin "///WO COS \p 

+ K^SS"//// ~ ^cc) cos "////A - (^.SC"//// + &Cs ) sin "///WO sin <A (6.89) 

* * ^ 

2 2 

d IIH = [ - e 0& n HH sin ^////'A - ^OC"//// cos "///WO 

+ [ _ (^cs("//// + 1 ) + ^SC"////) sin "///W' 
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+ {26 ss <jJnii + 1)} cos <^>////tA D cos \p 

+ [{2# (xPl/H + 0} s > n 

~{20cS O) IIH + ^5^ co //// + * )) cos w ////^] s i n ^ (6.90) 

The two formulations can be made equivalent through the use of a trig- 
onometric identity. 

6 sin(cot j/ + </>)= 6( sin cut j> cos $ + cos an)/ sin 4>) 

Using the above identity, the two formulations can be seen to be equivalent if: 
6 0S= 6 0 COS( t , o d CC= e C 5in( t>c 

0OC = ^O sin 0o Qss= ^5 cos< ^j 

0 = 0 £ cos (f) £ ^ SC ~ ® S 0 s (6.9 1 ) 

For the closed-loop HHC implementation, the control vector 0 will be defined 
as: 

» ={00, «0c l>c, «cc 0,s »,c) T (6.92) 

or an appropriate subset of this vector. 
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6.5.2 Minimization Objectives 


From the calculations described and presented in Ch. 5, rotor hub forces 
and moments are available in three orthogonal directions. In the absence of 
any modeling of the fuselage dynamics, these give the best available indication 
of helicopter vibration levels. Sine and cosine components of the 4/'rev. forces 
and moments are available as well as their peak-to-peak values. Any combi- 
nation of these quantities could be used as an indication of vibration levels and 
therefore be chosen to be minimized. In this study the sine and cosine com- 
ponents of the 4/rev. hub forces in the vertical, lateral, and longitudinal di- 
rections are chosen to be minimized. Decreases in vibration levels will be given 
in terms of the peak-to-peak values of these forces. For closed-loop HHC 
studies, using the notation of Ch. 5, the vector of hub vibratory shears Z will 
be defined as: 

^ = ^v/j4 ^v?c4 ^lts4 ^ltc4 ^lns4 ^lnc4^ (6-93) 

or an appropriate subset of this vector. 
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Chapter VII 

RESULTS SHOWING EFFECTS OF UNSTEADY AERODYNAMICS 
7.1 BASELINE BLADE CONFIGURATIONS 

The results presented in this chapter were calculated using two baseline 
blade configurations. A soft-in-plane and a stiff-in-plane hingeless configura- 
tion identical to those used in Ref. 8 were used. All results in this chapter ex- 
cept those comparing the trim procedures were calculated using the flap trim 
procedure. 

The baseline soft-in-plane blade has fundamental, rotating, uncoupled na- 
tural frequencies in lag, flap, and torsion of 0.735/rev, 1.123/rev, and 3.17/rev 
respectively. The rotor consists of four blades. The Lock number is y = 5.5, 
the rotor solidity is o = 0.07, and the rotor thrust coefficient is C T = 0.005. 
The blade has zero precone (} p , zero root offset e , , and no built in twist. The 
offset between the elastic axis and the center of gravity, x h and the offset be- 
tween the elastic axis and the aerodynamic center, x A , are both zero. All blade 
properties are uniform across the span. The blade chord is c = 0.055 R , the lift 
curve slope is a = 2n, and the profile drag coefficient is C m = 0.01. 

The baseline stiff-in-plane blade is identical to the baseline soft-in-plane 
blade except that the fundamental rotating lag frequency is 1.42/rev. 
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7.2 COMPARISON WITH PREVIOUS UNSTEADY RESULTS 


In Ref. 13 finite-state time-domain unsteady aerodynamics were incorpo- 
rated in a simple rigid blade, offset hinged spring model of a blade with flap 
and lag degrees of freedom. The intent of this study was to simulate the same 
soft-in-plane blade described in the last section, except that the torsional de- 
gree of freedom was not included. Unsteady effects were incorporated using 
two augmented states governing the unsteady aerodynamics at a typical sec- 
tion at 3/4R. These two augmented states were solved for as explicit state 
variables along with those for flap and lag. 

The present study uses an elastic hingcless blade with coupled flap-lag-tor- 
sional degrees of freedom which is a much more realistic model of an actual 
blade. The current model can be made to give results quite similar to the mo- 
del of Ref. 13 by increasing the torsional stiffness of the blade so that it re- 
sembles a torsionally rigid blade. With this in mind, a model identical to the 
baseline soft-in-plane blade, but with a fundamental torsional frequency of 
6.4/rev. was used to determine aeroelastic response and stability of the blade 
using the unsteady aerodynamic formulation of Ref. 13. This blade will be 
referred to as the baseline torsionally stiff blade. 

Figure 7.1 presents the flap and lag blade tip responses as a function of 
azimuth. Plots taken from Ref. 13 for the rigid offset hinged spring restrained 
soft-in-plane blade are also included. The flap responses from Ref. 13 and for 
the torsionally stiff blade are quite similar while the lag responses show the 
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same variation with azimuth but differ considerably in magnitude. These dif- 
ferences are reasonable since the structural models are completely different. 
Decreasing the torsional stiffness to a more realistic level, that of the soft-in- 
plane blade, leads to a considerably different flap response. 

Figure 7.2 presents the stability in the first flap and first lag modes for the 
same three cases. Basic stability trends are the same in all three cases and the 
instability of Ref. 13 is present in the first flap mode for each model. Although 
the first flap eigenvalue branch point for the soft-in-plane blade is at a slightly 
higher advance ratio, all three cases become unstable at the same advance ratio 
of n — 0.45. 

As was explained in Chapt. 4, the instability of Ref. 13 was due to an error 
in formulation which led to a singularity in flap loads at high advance ratios. 
The formulation of this study eliminates this singularity and in addition ap- 
plies the aerodynamics in a more realistic manner. 

Figure 7.3 presents the blade tip response in flap, lag, and torsion, at an 
advance ratio of n = 0.4 , for the soft-in-plane blade using the aerodynamic 
formulation of Ref. 13 and the inproved formulation of this study. Because 
the singularity of Ref. 13 was sharply localized at fi — 0.45, even at an advance 
ratio of n — 0.4 the vertical aerodynamic loads are reasonable. Much of the 
difference in response is due to the new formulation allowing the augmented 
states to vary along the blade. 
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Figure 7.4 shows the variation in first flap and first lag eigenvalues with 
advance ratio. Lag stability differs very little for the two cases. The sudden 
instability at p = 0.45 found using the aerodynamic formulation of Ref. 13 is 
eliminated using the formulation of this study. 


7.3 COMPARISON OF QUASISTEADY AND UNSTEADY RESULTS 
7.3. 1 Response and Stability 

Response and stability calculations were carried out for the baseline soft- 
in-planc and stiff-in-plane blades using the quasisteady aerodynamics of Ref. 
8 and the unsteady formulation of this study. Blade tip response in flap, lag, 
and torsion have been plotted for the two cases at advance ratios of p = 0.2 
and p = 0.4. Figures 7.5 and 7.6 show the responses for the soft-in-plane blade 
at p = 0.2 and p = 0.4. Figures 7.7 and 7.8 show the same responses for the 
stiff-in-plane blade. 

It can bee seen that the unsteady effects produce only moderate differences 
in the flap and lag responses. The response in the presense of the unsteady 
loads tends to lag behind that with quasisteady loads due to the phase lag and 
amplitude modulation associated with the unsteady aerodynamics[13]. As is 
evident from Figs. 7.6 and 7.8, the influence of the unsteady aerodynamics on 
the torsional response is more pronounced. This is due in part to the inclusion 
of several higher order apparent mass terms in Eqn. 4. 17, -the expression for 
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the noncirculatory aerodynamic moment in the unsteady aerodynamic formu- 
lation, which were not included in the quasistcady formulation. 

In order to determine the importance of these terms, response calculations 
were carried out at an advance ratio of n = 0.4 using a quasisteady aerodyna- 
mic formulation which included these apparent mass terms. Figures 7.9 and 
7.10 show the flap, lag, and torsional response using quasisteady aerodynamics 
without the apparent mass terms, quasisteady aerodynamics with the apparent 
mass terms, and unsteady aerodynamics. It can be seen that the apparent 
mass terms effect all three responses to some degree, but that there is a signif- 
icant change in the torsional response. This indicates that these apparent mass 
terms can be important and should not be neglected as they have been in may 
previous aerodynamic formulations. 

For the case of hover, there is no unsteadiness in the aerodynamics and the 
apparent mass terms are zero. The response therefore is identical to that of the 
quasistcady aerodynamics of Ref. 8. There are slight differences in hover sta- 
bility due to the unsteady effects on time derivatives in hover. 

Stability plots for the soft and stiff-in-plane blades are shown in Figs. 7.1 1 
through 7.14. The variation with advance ratio of the real part of the charac- 
teristic exponent associated with each mode is plotted. The influence of time 
domain unsteady aerodynamics on lag and torsional mode stability is small. 
This is reasonable since the primary correction in the unsteady aerodynamics 
theory is to the lift. The influence of the unsteadiness on the drag, which in- 
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flucnccs the lead lag motion, is small. There is however a noticeable effect of 
unsteady aerodynamics on the flap modes, particularly at higher advance ra- 
tios and with the soft-in-plane blade. 


7.3.2 Hub Vibratory Loads 

In Ref. 8, 4/rev. vertical hub shears were calculated using quasisteady 
aerodynamics. The blade structure was then optimized to reduce these shears. 
As described in Ch. 3, Ref. 8 used a less precise convergence criterion than the 
convergence control parameter used in this study when determining whether 
convergence of quasilinearization had occured. Furthermore, less stringent 
local error bounds were used when solving the differential equation associated 
with quasilinearization. The procedure employed in Ref. 8 was adequate for 
determining the low frequency response and stability of the blade, but led to 
response solutions which were not actually converged in the higher harmonics. 
This led to incorrect values of the 4/rev. shears which may have contributed 
to the inconsistent behavior of the optimization results. 

As an example, Fig. 7.15 shows the variation with advance ratio of the 
4/rev. vertical and lateral hub shears for the soft-in-plane blade as calculated 
in Ref. 8 and their comparison to the response solution obtained in this study. 
In both cases the same quasisteady aerodynamic formulation was used. The 
response solution obtained in this study was calculated by carrying out one 
additional iteration of quasilinearization beyond that used in Ref. 8 with in 
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addition much tighter local error bounds. From Fig. 7.15 it can be seen that 
the not completely converged response solution produces substantially differ- 
ent vertical hub shears. With the converged response of this study, the pecul- 
iar dip in vertical hub shears at an advance ratio of n = 0.3 completely 
disappeared. The effect of this convergence problem on the lateral hub shear 
is much less severe, as is evident from Fig. 7.15. 

Hub shear and moment calculations were carried out for the baseline soft- 
in-planc and stiff-in-planc blades using the quasisteady aerodynamics of Ref. 
8 and the unsteady formulation of this study. Figures 7.16 through 7.21 plot 
the variation with advance ratio of the six different 4/rev. hub shears and 
moments. All shears and moments vary somewhat depending on the aero- 
dynamics used. The most significant variation is in vertical shears which show 
a large difference in magnitude which increases with advance ratio. This is 
reasonable since the primary correction in the unsteady aerodynamic theory is 
that due to the lift. 


7.4 COMPARISON OF OLD AND NEW TRIM RESULTS 

The flap trim and full trim procedures described in Ch. 3 were used to de- 
termine the collective pitch 6 0 , cyclic pitch 0 lr and d ic , inflow ratio X, and rotor 
angle of attack. a R which are required to trim the rotor. 

With the full trim procedure, essentially identical results were found using 
quasisteady and unsteady aerodynamics. At hover, the quasisteady and un- 
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steady aerodynamic models give identical air loads and so the trim conditions 
are identical. With the simplified aeroelastic model used in the full trim anal- 
ysis, only the constant and l/'rcv components of each generalized coordinate 
arc retained in forward flight. The high frequency motions which contribute 
to unsteady airloads are therefore not modeled and it is to be expected that 
nearly identical trim states should result. Furthermore it should be noted that 
the flap trim procedure does not have the provision for incorporating the un- 
steady aerodynamics. 

Figures 7.22 and 7.23 depict changes in the trim parameters with advance 
ratio for the flap trim and full trim procedures using the baseline soft-in-plane 
blade. From Fig. 7.22 the value of the collective pitch d 0 can be seen to differ 
most at hover and to become essentially identical at n — 0.3 or above. Con- 
versely the value of the inflow ratio k is nearly identical up to n = 0.1 and then 
the difference increases somewhat at higher advance ratios. The rotor angle 
of attack a R varies considerably above hover, but since a R is not an explicit 
variable in the quasilinearization routine, this does not directly effect the rotor 
response. From Fig. 7.23 the cosine cyclic Q Xc varies only slightly between the 
two procedures while the sine cyclic 6 Xs differs somewhat at advance ratios 
above ^ = 0.2. 

Overall the two trim procedures produce quite similar results for the soft- 
in-planc blade. It is expected that more realistic blades, having for instance 
cross-sectional offsets and non-uniform properties, which cannot be modeled 


118 


by the flap trim procedure, will show more significant differences between the 
two procedures. 

Figure 7.24 shows blade tip response for the flap trim and full trim proce- 
dures applied to the soft-in-plane blade. In each case the trim values deter- 
mined from the trim procedure were used in the full aeroelastic model with 
unsteady aerodynamics at an advance ratio of n =0.3. The nonlinear response 
solution was determined and the nondimensional tip displacements in flap, lag, 
and torsion arc plotted versus blade azimuth. Lag and torsional response vary 
little between the two cases, however the flap response changes considerably. 
The change in magnitude of the flap response is due primarily to the change 
in the inflow ratio X at this advance ratio, while the different value of the sine 
cyclic pitch changes the phase of the response somewhat. 

The two trim procedures led to almost identical stability eigenvalues for the 
blade. The largest difference in value of the real part of the eigenvalue for any 
mode at any advance ratio was less than 2%. Therefore these results are not 
presented here. 

The 4/rev. hub shears and moments were calculated for the soft-in-plane 
blade using the response solution obtained with both the flap trim and the full 
trim. Figures 7.25 through 7.27 dipict of the variation, with advance ratio, of 
the six shears and moments. Application of the coupled trim procedure re- 
sulted in somewhat higher hub shears and somewhat lower hub moments than 
those due to flap trim. The largest difference was in the vertical hub shears. 
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This behavior is reasonable since the vertical hub shears are dependent on the 

flapping response which was shown to differ considerably for the two trim 
formulations. 
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Chapter VIII 

OPEN LOOP RESULTS 


8.1 INTRODUCTION 

This section presents the results obtained during the open-loop application 
of HHC. The primary objective here is to gain an understanding of basic 
trends in rotor response to HHC and to provide initial estimates of the HHC 
transfer matrix T which is required for starting the closed-loop HHC algo- 
rithms. The baseline soft-in-planc blade, as described in Ch. 7, was used in 
all these calculations. Trim was calculated using the full trim procedure. 


8.2 HUB LOAD VARIATION WITH HHC PHASE 

As shown in Ch. 6, the HHC inputs can be expressed as a sine function with 
a phase offset, applied to each of the collective, lateral, and longitudinal con- 
trol channels. A 1/3° HHC pitch input of the form: 

^///y = 1/3° s irt(i// — </>////) (8-1) 

was applied to each of the three control channels with the phase 4> l{ll being 
allowed to vary from 0° to 360° in 90° increments. These calculations were 
carried out for both the quasisteady and the unsteady aerodynamic models. 
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Figures 8.1 through 8.9 present the variation in hub shears and moments 
with HHC phase angle in the three control channels for quasistcady and un- 
steady aerodynamics. Generally for a given shear or moment, if HHC in any 
of the three control channels is phased correctly, a minimum value of the shear 
or moment can be obtained. Except for the case of yawing moments, and 
vertical shears with collective or lateral inputs, this minimum is below the 
baseline value. 

The major difference in responses to HHC between quasisteady and un- 
steady aerodynamics are in the vertical shears. For all the control channels, 
the HHC phase at which minimum vertical shears with unsteady aerodyna- 
mics are obtained, leads the value calculated with quasistcady aerodynamics 
by about 90 . In the case of collective and lateral HHC input, the unsteady 
formulation show's shears above the baseline values for any phasing, w'hile the 
quasistcady formulation indicates a minimum below the baseline value. With 
longitudinal input both the aerodynamic models show' substantial decreases 
below the baseline for vertical shears with correct HHC phasing. In actual 
flights tests[67] also, longitudinal control has been shown to be most effective 
in alleviating vertical shears. 

Rolling and pitching moments show similar behavior of the responses with 
both aerodynamic formulations and for control inputs in all three channels. 
Lateral and longitudinal shears show similar response patterns for the tw'o 
formulations. The responses have similar phases and magnitudes of excursion 
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from the baseline but are offset by the difference in baseline peak-to-peak 
shear values. 


8.3 LINEARITY OF RESPONSE TO HHC 

Linear optimal control theory assumes linear response to control input. It 
is therefore important that actual response of the model be at least fairly linear 
over some reasonable range for the HHC to be most effective. The assumption 
used in the control formulation is that the 4/rev. sine. and cosine harmonics of 
the hub shears vary linearly with the 4/rev. sine and cosine HHC harmonic 
inputs. To verify this assumption, HHC inputs of -1/3° through 3° in 2/3° in- 
crements were applied to the 8 CC HHC degree of freedom, i.e. the cosine com- 
ponent in the longitudinal control channel. The calculations were carried out 
for quasistcady and unsteady aerodynamics at an advance ratio of /r =0.3. 

The variations in the sine and cosine harmonics of the 4/rev. hub shears 
with HHC input pitch magnitude are shown in Figs. 8.10 and 8.11. The re- 
lationship between control harmonics can be seen to be quite linear over this 
range for both aerodynamic formulations. The difference in response between 
the two aerodynamic models can be seen to vary from slight offsets to sub- 
stantially different slopes and intercepts. Thus the transfer matrices needed 
for control implementation can be expected to differ greatly for the two for- 
mulations. The variation of peak-to-peak hub shears with HHC magnitude 
are showm in Fig. 8.12. The peak-to-peak shears are a function of the squares 
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of the sine and cosine harmonics and thus arc not linear with HHC magnitude. 
As can be seen, they show distinct minima, with HHC magnitudes above or 
below this value giving higher shears. 


8.4 TRANSFER MATRIX CALCULATION 

Initial estimates for the HHC transfer matrix at the baseline were calculated 
by incrementing each of the six control degrees of freedom by 1/3°. From Eqn. 
6.1 and 6.2 it can be seen that at the baseline, with no initial HHC input, both 
the local and the global HHC models reduce to: 

AZ=7A0 (8.2) 

From this it can be seen that, assuming low process and measurement noise, 
the columns of the T matrix can be calculated by incrementing each element 
of 9 individually, calculating the changes in hub shears due to this increment, 
and dividing their values by the magnitude of the control increment. For in- 
stance setting the first element of 9 to one and all other elements to zero, and 
multiplying this vector times the T matrix, will result in a vector identical to 
the first column of the T matrix. Baseline transfer matrices were calculated 
using quasisteady and unsteady aerodynamics and are given in tables 8.1 and 
8 . 2 . 

The calculated baseline transfer matrices for the quasisteady and unsteady 
aerodynamic cases differ greatly. In order to see basic trends for the two cases, 
the bar chart given in Fig. 8.13 has been constructed showing the proportional 
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TABLE 8.1 


Quasisteady HHC Transfer Matrix 


0.3145325 

-0.2216286 

-0.0680207 

-0.0569284 

0.0041394 

0.0734329 


-0.2224358 

-0.3188984 

-0.0575412 

0.0649546 

0.0806664 

-0.0100913 


-0.0447006 

0.0572526 

0.0156767 

0.0097563 

-0.0532560 

0.0120377 


0.0600985 

0.0441021 

0.0122215 

-0.0127674 

0.0065953 

0.0480450 


-0.1066384 

-0.2526287 

-0.0879807 

0.0093927 

-0.0125365 

-0.0136280 


-0.251 1 157 
0.1064007 
0.0091254 
0.0814574 
-0.0249821 
0.0000364 


TABLE 8.2 

Unsteady HHC Transfer Matrix 


0.3975075 

0.0018428 

-0.0298565 

-0.0408923 

-0.0094324 

0.0422244 


-0.0000851 

-0.3944278 

-0.0405612 

0.0294869 

0.0419992 

0.0102764 


-0.0326116 

0.0092679 

-0.0071383 

0.0244328 

-0.0516340 

0.0220154 


0.0052673 

0.0358874 

0.0244647 

0.0068263 

0.0210999 

0.0528939 


0.0276221 

-0.2285797 

-0.0822938 

0.0067259 

0.0071792 

-0.0321438 


-0.2291952 

-0.0251755 

0.0075192 

0.0820623 

-0.0322288 

-0.0064157 
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changes in shears due to each of the 1/3° HHC incrementations. From this 
chart some general observations can be made. The vertical shears arc effected 
much more by HHC than the lateral or longitudinal shears. The difference in 
response with quasisteady and unsteady aerodynamics is greatest for vertical 
shears which are most directly dependent on the aerodynamic lift where the 
two formulations differ. Collective and longitudinal HHC inputs have a 
greater relative effect on shears than lateral HHC inputs. 


126 



Chapter IX 

CLOSED-LOOP CONTROL RESULTS 
9.1 INTRODUCTION 

One of the main objectives of this study was to reduce the vertical, lateral, 
and longitudinal 4/rev. hub shears. Sine and cosine components of each of 
these three shears were to be minimized by the control algorithms. In applying 
the control algorithms to the quasilincarization model, a control vector was 
first calculated using the baseline response. The vibratory response of the ro- 
tor with this control was then found and this information was used in calcu- 
lating an improved control vector. Physically this is equivalent to allowing all 
transients to die out between control changes. The response due to the new 
control values was then calculated and the process was continued for a given 
number of iterations or until no further decrease in shears was observed. 

The baseline soft-in-plane blade described in Ch. 7 was used in all these 
studies. Unless otherwise stated, all the control calculations were done at an 
advance ratio of ^ = 0.3. Trim was calculated using the full trim procedure. 
The initial baseline transfer matrices for the quasistcady and unsteady aero- 
dynamic models at this advance ratio were those calculated in Ch. 8. 
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9.2 FIXED GAIN CONTROL 


The baseline HHC transfer matrices were used to construct Fixed gain con- 
trollers. For these controllers, there was no identification of control parame- 
ters and the hub shears were assumed known exactly. Control was therefore 
determined by a simple inverse relation: 

W+ D= T~ l Z(i) (9.1) 

This assumes the local HHC model for multistep application. The first step is 
the same as a global controller, which is a one step process when parameters 
cannot change through identification. 

With quasisteady aerodynamic modeling the HHC input angles, as shown 
in Fig. 9.1, moved smoothly to their optimums with almost no overshoot or 
oscillation. The peak-to-peak values of the three shears were reduced smoothly 
to less than 2% of their baseline values within five iterations as shown in Fig. 
9.2. With unsteady aerodynamic modeling the HHC input angles also moved 
smoothly to their optimums, as shown in Fig. 9.3. With unsteady aerodynamic 
modeling the hub shears were reduced to less than 5% of their baseline values 
in two iterations but the vertical shear had risen to 10% of its baseline value 
by the fifth iteration. 

These results indicate that both transfer matrices were quite good estimates 
of the actual relationship between HHC input and hub shear harmonics at the 
optimal HHC angles. The increase in vertical shear with continued iterations 
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for the unsteady ease indicates that this transfer matrix was not quite as good 
an estimate as the quasisteady transfer matrix. 


9.3 ADAPTIVE CONTROL AT DESIGN FLIGHT CONDITION 

The minimum variance controller described in Ch. 6 provides for one local 
and three global controller designs using different Kalman identification 
schemes. These were referred to as controllers I, 2, 3, and 4. Each of these 
controllers may be implemented as a deterministic, cautious, or dual controller. 
Throughout this study equal weighting was given to the different shear com- 
ponents and no weighting was applied to the control input magnitudes or rates 
of change. It should be stressed that in these simulations a change in step in- 
dex from i to i + 1 indicates an iteration of quasilinearization rather than any 
real time increment. Also the Kalman filter identifies computer code output 
parameters which move in an unpredictable manner but does not have to deal 
with any articficially introduced noise. 

In the initial phase of this investigation the deterministic versions of the four 
controllers were were implemented for the quasisteady and unsteady aeroelas- 
tie models. For these calculations the Kalman filters were initialized by as- 
suming the variances IV 0 and W T , and the initial values of the elements of the 
covariance matrix P, to be equal with a value of 0.5* 10~ 2 , and assuming the 
variance V to have a value of 0.3T0 -7 . 
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Figures 9.5 through 9.7 give the iteration history for the three hub shear 
conponcnts for the four controllers applied to the quasisteady model. In each 
case the First step is identical for deterministic control. The results for cases I 
and 2 arc quite similar as are the results for cases 3 and 4. When these con- 
trollers were applied to the unsteady model, all four controllers produced al- 
most identical hub shear histories with only vertical shear showing some 
spread in the responses. All four controllers gave staisfactory results and of 
the global controllers, controller 4 appeared to give the best results after several 
iterations. It was decided therefore to continue with more extensive studies 
using the deterministic and cautious versions of just two controllers, controller 
1, a local controller with identification of only the HHC transfer matrix, and 
controller 4, a global controller with identification of the baseline vibrations 
and the HHC transfer matrix. In these additional studies the Kalman filter 
was initialized by assuming the variances V, W, and W 0 to be equal with a 
value of 0.3x10“ 7 and by initializing the elements of the covariance matrix P 
as 0.3x 10“ 5 . 

With quasisteady aerodynamics the difference in control sequence between 
the deterministic and cautious controllers was relatively small, reflecting the 
excellent initial transfer matrix. Figures 9.8, 9.9, and 9.10 show the iteration 
history of the six HHC inputs and the three hub shears for deterministic and 
cautious versions of the local controller. Figures 9.11 through 9.13 show the 
same iteration histories for the global controller. All four controllers led to the 
same control solution for hub shear minimization. 
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With unsteady aerodynamics the cautious and deterministic controllers 
show more difference in control values initially but all four quickly converge 
to the same control values. The deterministic controller tends to slightly over- 
shoot the optimal control on the first iteration. Figures 9.14, through 9.19 
show the iteration histories of the HHC inputs and hub shears for the local 
and global controllers with unsteady aerodynamics. 

When applied to either the quasisteady or unsteady model, HHC was able 
to eliminate essentially all the 4/rev. hub shears. With local control applied to 
the unsteady model, the vertical shears rise slowly after the second iteration 
indicating that the transfer matrix identification has not been ideal. Since the 
control input is changing only slowly and identification of the local T matrix 
is based on changes in control, this result can be easily anticipated. 

Figures 9.20 and 9.21 compare the final hub shear values after five iter- 
ations for the four control cases as a percentage of their baseline values. When 
unsteady aerodynamic effects are included the vertical shear becomes the crit- 
ical component, being decreased much less than the lateral and longitudinal 
shears. With quasisteady aerodynamics there is no discernible pattern of one 
shear being reduced more than another. The unsteady baseline vertical shear 
started at 1/3 the level of the quasisteady baseline, so in absolute terms rather 
than as percentages of the baseline, the global controllers produced similar 
vertical shear magnitudes using the two aerodynamic formulations. The local 
controllers however were much less successful when used with the unsteady 
model than they were with the quasisteady model. Figure 9.22 compares the 
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optimal HHC inputs after five iterations of the local controller using quasis- 
tcady versus unsteady aerodynamics. Completely different, but similarly ef- 
fective, control inputs are needed for hub shear reductions depending on the 
type of aerodynamic representation used in the aeroelastic model. 

The LMS Direct Inverse Control algorithm described in Ch. 6 provides a 
simple alternative method of determining HHC inputs for hub shear minimi- 
zation. This algorithm was implemented using the unsteady aerodynamic 
model. In Ref. 32, values of the gain term K s of between 0 and 5 were used. 
A very simple matrix representation of rotor response to HHC was used, and 
often on the order of 100 or more iterations were carried out for convergence. 
When a value of K s = 0.3 was used to implement this control, the results in 
terms of HHC angle and shear minimization iteration histories were essentially 
identical to those results already presented for fixed gain HHC calculations 
and shown in Figs. 9.3 and 9.4. This was because changes in the identified C 
matrix were too slow to be detectable within five iterations. It was necessary 
to increase K s far beyond the range in which stability problems were found in 
Ref. 32 before appreciable changes in convergence within five iterations were 
found. Unfortunately the high computational cost of the aeroelastic analysis 
used here precludes the investigation of convergence characteristics much be- 
yond the fifth HHC iteration. 

In order to determine the effect of HHC on overall blade response, the 
clastic tip deformations in flap, lag, and torsion were plotted for the baseline 
response and for the optimal response for the case of controller 4, the global 
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cautious controller. From Fig. 9.23 it can be seen that the principal differences 
in elastic response are in the flap and torsional degrees of freedom. The tor- 
sional response is somewhat missleading in that the actual tip angular de- 
flection includes also a rigid body rotation due to the total blade root angle 
associated with trim and HHC, and a rotation of the form due to the 

bending of the blade. When these geometric rotations are added to the tip 
elastic torsional deflection, the total tip angular displacement shown in Fig. 
9.24 results. The term accounts for less than 2% of this total angular 
displacement. This plot indicates that vibration reduction through HHC is 
achieved not by reducing overall blade response, but by introducing modified 
airloads on the blade which cancel out the vibratory hub loads. 


9.4 CONTROL WITH STEP CHANGE IN FLIGHT CONDITION 

A practical controller must be able to adapt to changing flight conditions. 
As a test of this ability a step change in flight condition from /i = 0.3 to 
H = 0.35 was applied to the deterministic and cautious versions of the local and 
global controllers. This was done by starting with the converged optimal sol- 
ution and its response at fi = 0.3, changing the propulsive trim values to those 
for /r = 0.35, and proceeding with iterative control calculations and quasili- 
nearization solutions. 

When this procedure was applied to the local controller there were large 
oscillations in the calculated control inputs and the resultant hub shears. As 
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can be seen from Figs. 9.25 and 9.26, these oscillations had damped out con- 
siderably by the sixth iteration and all the shears were considerably below the 
baseline values for this advance ratio. 

When this step change in fx was applied to the global controller, there was 
little oscillation in the required control inputs and the controller moved fairly 
smoothly from an initial large increase in shears toward a minimum. Figures 
9.27 and 9.28 show the iteration history of the HHC inputs and hub shears for 
this case. The procedure was stopped at five iterations but the shears were 
continuing to decrease at that point. 

A comparison of the three shear components and their baseline values for 
the local and global controllers is given in Fig. 9.29. As can be seen, the global 
controller has been more successful in reducing shears. Additional iterations 
of the analysis would be expected to reduce the shears to lower levels. The less 
effective performance exhibited by the local controllers is attributed to the 
method employed for the calculation of the transfer matrices, which is based 
on changes in control input rather than the total control magnitudes. Since the 
controllers quickly reached minimum shear levels from the baseline, changes 
in control became very small and the transfer matrices were not completely 
identified. When a sudden change in flight condition occurred this uncertainty 
degraded the controller performance. 
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Chapter X 

HHC FOR ARTICULATED VERSUS HINGELESS ROTORS 
10.1 INTRODUCTION 

In this chapter the effects of HHC on a hingclcss rotor and a roughly 
equivalent articulated rotor are investigated. This includes comparisons of re- 
quired control angles, shear minimization, accompanying changes in hub mo- 
ments, attempts at shear and moment minimization, changes in rotating blade 
loads, changes in blade stability, and changes in power requirements. 

The hingeless rotor used in the study is the same as the baseline soft-in- 
plane rotor described in Ch. 4 except that a blade root offset from the center 
of rotation of 5% of the blade elastic length has been introduced. The artic- 
ulated rotor used was produced from this rotor by introducing hinges in flap 
and lag at the same 5% blade root offset. Additionally a rotational damper 
was applied to the articulated blade root lag hinge as is generally required to 
maintain stability in articulated rotors. This was a viscous damper and gave 
a nondimcnsional damping coefficient of: 

C = — ^—=0.02 (10.1) 

m (file 

The value 0.02 was chosen as being approximately the value of the damper on 
an OH6A helicopter rotor blade. 
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The hingelcss blade had fundamental uncoupled rotating natural frequen- 
cies in lag, flap, and torsion of 0.783/rev., 1. 158/rev., and 3.176/rev. The ar- 
ticulated blade had natural frequencies in lag, flap, and torsion of 0.274/rev., 
1.037/ rev., and 3. 176/rev. Mode shapes for the first six modes were very sim- 
ilar for the two blades except in the area of the blade root where the two 
structural boundary conditions were different. 

Trim calculations were carried out for the hingeless and articulated blades 
at an advance ratio of n — 0.3 using the full trim procedure. The required col- 
lective and sine cyclic angles were almost identical for the two rotors. The 
hingeless rotor required a rotor shaft angle 1% higher and a cosine cyclic angle 
10% lower than that of the articulated rotor. 

Using these trim values, baseline response and stability were calculated for 
the articulated and hingeless blades at an advance ratio of fx =0.3. Subsequent 
HHC calculations were carried out from these baselines. All calculations in 
this chapter were carried using unsteady aerodynamics. 


10.2 REDUCTION OF HUB SHEARS 

From the baseline flight condition at an advance ratio of /r=0.3, initial es- 
timates of the HHC transfer matrices for the hingeless and articulated rotors 
were calculated using the method described in Ch. 8. Using these initial 
transfer matrices, controller 4, the cautious global controller of Ch. 6, was ap- 
plied to reduce 4/rev. vibratory hub shears. Five iterations of the shear min- 
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imization procedure using HHC were performed for each rotor. All results in 
this chapter used this same HHC algorithm. 

Figure 10.1 shows the iteration history of the HHC angles required to re- 
duce hub shears for the hingelcss rotor. The maximum HHC angle required 
was 4°. By the second iteration the HHC angles had essentially reached their 
optimum values. Figure 10.2 shows the iteration history of the 4/rev. hub 
shears being minimized. All three hub shears were greatly reduced. The lat- 
eral and longitudinal hub shears were reduced to less than 2% of their baseline 
values by the fifth iteration. The vertical shear, which started at one fifth the 
level of the inplane shears, was reduced to 13% of its baseline value. While 
the shears were being minimized, the hub vibratory moments increased greatly. 
Figure 10.3 shows the iteration history of these moments. The rolling moment 
and pitching moment rose to respectively 6 times and 4 times their baseline 
values when HHC was applied. The yawing moment started at a very small 
value and rose to 250 times its baseline value. The yawing moment does not 
act directly on the fuselage but acts indirectly through the helicopter power- 
plant. Typically rotor torque is developed through complicated fluid forces 
acting on the power turbine of a gas turbine powcrplant. Therefore the sig- 
nificance of this large increase in yawing vibratory moment is not at all obvi- 
ous. 

Figure 10.4 shows the iteration history of the HHC angles required to re- 
duce hub shears for the articulated rotor. The maximum HHC angle required 
to minimize the hub shears was 0.9° but intermediate values of up to 1.1° were 
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used. By the second iteration the HHC angles had essentially reached their 
optimum values. The HHC angles required for hub shear reduction with the 
articulated blade were much lower than those required with the hingelcss 
blade. Figure 10.5 shows the iteration history of the 4/rev. hub shears being 
minimized for the articulated blade. The lateral and longitudinal shears were 
reduced to less than 0.4% of their baseline values. The vertical shear was ini- 
tially much lower and was reduced to 9% of its baseline value. 

Figure 10.6 shows the iteration history of the 4/rev. hub vibratorv moments 
as the shears were being minimized. For the articulated rotor, the zero mo- 
ment condition at the blade root flap hinge and the low vibratory moment 
transmitted by the rotational damper at the lag hinge kept the vibratory 
pitching and rolling hub moments from rising greatly as they had for the 
hingeless blade. The greatest increases in moment were when HHC was first 
applied. At the first iteration of HHC the hub rolling, pitching, and yawing 
moments rose respectively to 1.6, 3.4, and 4.0 times their baseline values. 
These values for the pitching and rolling moments were still only slightly above 
the baseline values for the hingeless rotor. From the second HHC iteration 
onward the hub rolling, pitching, and yawing moments were 1.3, 0.8, and 2.2 
times their baseline values. These values were an order of magnitude lower 
that those for the hingeless rotor. The moderate peaks in moment on the first 
HHC iteration were associated with the first, poor approximation to the opti- 
mum HHC angles and it is felt they could be avoided by tuning the controller 
to give a more cautious first estimate of HHC requirements. 
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10.3 REDUCTION OF SHEARS AND MOMENTS 


Vibrations at a point in a helicopter fuselage will depend on both the hub 
shears and the hub moments being fed from the rotor into the fuselage. In 
actual implcmcntation[68], the quantities being minimized by the HHC algo- 
rithm would normally be three linear orthogonal acceleration components at 
some point in the fuselage such as the pilot's scat. Due to the offset of this 
point from the rotor hub, angular accelerations are transformed into linear 
accelerations, and therefore hub vibratory moments also contribute to the ac- 
celerations experienced at the pilot's seat. The dynamic properties of the fu- 
selage will also act as a filter changing the phase and magnitude of the forces 
and moments felt at this point due to the hub loads. Because of these factors, 
minimizing vibrations at a specific point in the fuselage may actually increase 
some loads at the rotor hub. Conversely, trying to reduce all hub loads may 
not be the most effective way to reduce vibrations at a specific point in the 
fuselage. The model used in this study is only capable of modeling hub shears 
and moments and has no provision for calculating fuselage vibrations. 

In the last section it was seen that minimizing hub shears could lead to 
greatly increased hub moments, particularly for the hingeless blade rotor. It 
is therefore interesting to attempt to minimize both the hub shears and the hub 
moments simultaneously. The same cautious global HHC algorithm is used 
to reduce the hub loads, however now the problem is slightly redefined: 

Z(i + 1 ) = Z 0 + 710(0 ( 1 0.2) 
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Where the vibration vector Z is now of length 12 and consists of the sine and 
cosine harmonics of the three 4/rev. hub shears and the three 4/rev. hub mo- 
ments. The transfer matrix T is therefore redefined to be a 12x6 matrix re- 
lating these 12 hub load harmonics to the 6 HHC input harmonics. There arc 
now 12 hub load harmonics being minimized using only 6 HHC harmonics as 
input, so the number of quantities being minimized is larger than the number 
of input variables. 

Figure 10.7 shows the iteration history of the HHC angles required by the 
12x6 HHC algorithm applied to the hingeless blade.' As with the 6x6 algo- 
rithm, the optimum HHC angles were essentially reached by the second iter- 
ation. The maximum HHC angle was 0.7°, which was much smaller than the 
4° commanded by the 6x6 algorithm. Figure 10.8 shows the iteration history 
of the hub vibratory shears. With the 6x12 algorithm, there was no minimi- 
zation of the vertical shear and the lateral and longitudinal shears were only 
reduced by 20% to 30%. Figure 10.9 shows the iteration history of the hub 
vibratory moments. The pitching and rolling moments were reduced to about 
40% of their baseline values and the yawing moment increased to about 24 
times its baseline value. 

The results for the 6x6 and 6x12 HHC algorithms are summarized in Fig. 
10.10 where they are compared to the baseline values for the hingeless blade. 
It can be seen that including the requirement for hub moment minimizations 
in the HHC algorithm eliminated the high hub moments of the 6x6 controller, 
but at the expense of very poor hub shear reductions. As was said earlier, this 
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docs not necessarily mean vibrations cannot be minimized in some area of the 
fuselage, but the implication is that this may be more difficult with a hingelcss 
rotor than with an articulated rotor for which vibratory hub moments are in- 
herently smaller. 

Figure 10.1 1 shows the iteration history of the HHC angles commanded by 
the 6x12 control algorithm for the articulated blade. By the second iteration 
the angles had essentially reached their optimum values. The maximum HHC 
angle commanded was 0.8°. Comparing Figs. 10.4 and 10.11 it can be seen 
that almost the same HHC input was commanded by the 6x6 and 6x12 con- 
trollers. Figure 10.12 shows the iteration history of the hub vibratory shears. 
The lateral and longitudinal shears were reduced to about 5% of their baseline 
values while the vertical shear was reduced to a similar numerical value but 
only 30% of its baseline value. Figure 10.13 shows the iteration history of the 
hub moments. Again there were moderate peaks in the pitching and yawing 
moments at the first HHC iteration. The hub moments then settled at nearly 
the same values as with the 6x6 algorithm. The pitching moment w r as 60%, 
the rolling moment 130%, and the yawing 170% of its baseline value. 

The results for the 6x6 and 6x12 HHC algorithms applied to the articulated 
rotor are summarized in Fig. 10.14. For the articulated rotor, the increases in 
hub moments associated with the 6x6 algorithm were small. Including mo- 
ment minimization in the control algorithm had relatively small effects, leading 
to slightly poorer shear reductions but slightly lower hub moments. For the 
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articulated rotor, the large trade off between shear and moment reductions 
observed in the hingeless rotor was not evident. 


10.4 ROTATING BLADE ROOT LOADS 

In the last two sections, changes in vibratory hub loads in a fixed, non-ro- 
tating reference frame were investigated as HHC was applied. It is important 
to determine whether minimization of non-rotating loads was attained at the 
expense of higher vibratory loads in the rotating system. Higher rotating 
peak-to-peak loads could lead to fatigue problems while higher maximum ro- 
tating loads could cause blade failure. With an articulated blade, the flap and 
lag hinges reduce blade bending moments to zero at the blade root. With a 
hingeless blade however, these moments could become undesirably large. 
Large increases in peak-to-peak blade root shears or feathering moment would 
also be undesirable. 

In order to investigate rotating blade root loads, the procedure described in 
Ch. 5 for calculating hub forces and moments was modified. A Fourier anal- 
ysis of blade root loads was carried out in the rotating system calculating forces 
and moments at the blade root offset rather than at the center of rotation, and 
the summation process over the four blades was deleted. This analysis pro- 
vided a constant component and sine and cosine components in the first five 
harmonics for axial, cordwise, and vertical blade root shear and feathering, 
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flapping, and lagging blade root moments. The maximum pcak-to-peak value 
and maximum absolute value of each load were also calculated. 

Figure 10.15 shows the variation with iteration number of peak-to-peak vi- 
bratory blade root shears for the hingeless blade. The vertical peak-to-peak 
shear increased by 40% as HHC was applied while the cordwise and axial 
shears increased by about 25%. Figure 10.16 shows the variation in blade root 
moments as HHC was applied. The peak-to-peak feathering moment in- 
creased to 10 times its baseline value while the flapping moment increased to 
1.7 times its baseline value and the lagging moment to 1.5 times its baseline 
value. 

Similar calculations were carried out for the articulated blade. For this case 
blade root shears increased by less than 5%. The blade root feathering mo- 
ment increased to 2.5 times its baseline value. For the hingelcss and to a lesser 
degree the articulated blade, the increases in peak-to-peak blade root feather- 
ing moment would need to be accounted for when determining the fatigue life 
estimates of the rotor blade on which HHC is applied. Furthermore, the 
moderate increases in peak-to-peak blade root bending moments and shears 
encountered for the hingeless blade could lead to fatigue problems. 

Figure 10.17 shows the variation in the maximum absolute value of the 
blade root moments as HHC was applied to the hingelcss blade. The maxi- 
mum feathering moment increased to 4 times its baseline value while the lag- 
ging moment increased by 50% and the flapping moment increased by 15%. 
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For the articulated blade again there were essentially no increases in maximum 
blade root shears while the feathering moment increased by 50%. These in- 
creased maximum absolute blade root loads, particularly for the hingeless 
blade, may imply the need for increased blade strength and cause weight pen- 
alties. 


10.5 POWER REQUIREMENTS AND STABILITY 

Operating the electro-hydraulic actuators needed to implement HHC will 
of course require power from the helicopter powerplant. In addition, the heli- 
copter rotor may require more or less power than at the baseline condition 
because of the additional aerodynamic loads which are imposed on it by the 
HHC inputs. In the flight tests of Ref. 67 a small decrease in rotor power re- 
quired was observed when HHC was applied. If any such decrease in required 
power could be relied on, it could compensate somewhat for the additional 
power required to run the HHC servo-actuators. 

When HHC was applied to the hingeless rotor, power required increased 
by 0.58% at the first iteration and remained between 0.90% and 1.03% above 
the baseline from the second through the fifth HHC iteration. This increase 
in required power was accompanied by a 0.2% increase in rotor thrust. In 
actual flight the pilot would retrim the helicopter after the application of HHC 
by decreasing the collective angle slightly to return to the baseline thrust value. 
This would alleviate part of the rise in required power seen in this simulation. 
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When HHC was applied to the articulated rotor, the rotor power required in- 
creased by 0.19% on the first iteration and then remained constant at 0.15% 
above the baseline for iterations two through five. The increase in required 
power was accompanied by a 0.06% decrease in rotor thrust, meaning that 
when the rotor was retrimmed the increase in required power would be some- 
what greater. The hingclcss rotor therefore had a six times greater increase in 
power required with application of HHC than the articulated rotor, but this 
difference was somewhat mitigated by the increase in thrust of the hingclcss 
rotor and decrease in thrust for the articulated rotor. In both cases the in- 
crease in required power was less than 1%. There was no indication of any 
decrease in power requirements when HHC w'as applied. 

Another possible problem area with the application of HHC would be any 
substantial decrease in the blade aeroelastic stability margins. With this in 
mind, the real parts of the characteristic exponents for the six generalized co- 
ordinates, which determine the stability of the six blade modes, were calculated 
at the baseline and at each iteration of the shear minimization process. The 
characteristic exponents for the articulated blade were all very insensitive to 
HHC application. The maximum changes in the real part of the characteristic 
exponent values were from -0.05% to +0.04%. For the hingclcss blade all but 
the second lag mode characteristic exponent had changes of between -0.34% 
and +0.58%. The real part of the characteristic exponent for the second lag 
mode, which was the closest to its stability margin, had an increase in stability 
of between 27% and 30%. Thus one can conclude that overall stability mar- 
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gins were not adversely effected by the application of HHC to either hingeless 
or articulated blades. 
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Chapter XI 

MODELING THE OH6A HELICOPTER 
11.1 INTRODUCTION 

This chapter presents the results of a simulation aimed at modeling the re- 
sponse to HHC of an OH6A light helicopter rotor. This helicopter was used 
for extensive flight testing of open-loop and closed-loop HHC by the McDon- 
nell Douglas Helicopter Company (MDHC, then Hughes Helicopters). The 
acroclastic simulation employed in this study has certain limitations in its 
ability to model the OH6A rotor as will be discussed later. However another 
a fundamental limitation of this aeroelastic analysis is the lack of any repre- 
sentation of the fuselage dynamics. Vibration data for the OH6A flight tests 
were recorded as accelerations at a particular point in the fuselage, namely the 
pilot's seat, therefore finding close correlation between the flight test data and 
the hub loads calculated in this study is not anticipated. The main objective 
of the simulations conducted was to identify similar trends, to those observed 
during the flight tests, when HHC was used. 

The OH6A is a light turbine powered helicopter with a fully articulated four 
bladcd rotor. Basic descriptive data for the helicopter and its rotor are given 
in Table 11.1. Geometric, structural, and mass properties for the OH6A rotor 
were provided by MDHC along with information on usual methods used in 
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TABLE 11.1 


OH6A Basic Data 


Gross Weight 

2550 lbs 

Equivalent Flat Plate Area 

5.0 ft 2 

Rotor Diameter 

26.33 ft 

Blade Chord 

6.75 in 

Blade Twist 

-9° (linear) 

Flap Hinge Offset 

5.5 in 

Lag Hinge Offset 

16.19 in 

Rotor RPM 

483 

Lock Number 

4.919 

Precone 

0° 


modeling the rotor. An illustration of the OH6A rotor hub is given in Fig. 
11.1. Plots were provided by MDHC for the cordwise position of the blade 
center of mass, shear center, and feathering axis. Additional data in the form 
of plots for the values of the blade chordwise, flapwise, and torsional stiff- 
nesses as a function of blade span together with the radial distribution of 
lumped blade mass and lumped blade chordwise mass moment, were also 
provided by MDHC. Several important aspects of the modeling process for 
this rotor configuration are presented below: 
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1. Control system stiffness is modeled by a torsional spring acting on the 
pitch case about the feathering axis. 

2. The lead-lag friction damper is usually modeled by a torsional dashpot 
acting between the blade root and the pitch ease about the vertical axis 
with an equivalent viscous damping coefficient of 1016 in-lbf- 
scc/radian. 

3. The pitchcase/strap pack chordwise clastic bending stiffness is modeled 
as an 884,800 in-lbf/radian torsional spring acting about the vertical 
axis between the root and flap/feather bearing. 

4. The pitch case flapwise elastic bending stiffness is considered infinite. 

5. The masses of the strap pack, pitch link, pitchcase, and damper are 
lumped at the lag hinge to give an equivalent first mass moment about 
the center of rotation. 

Natural rotating frequencies and mode shapes were provided for model 
verification. These were calculated at a collective angle such that the 3/4 ra- 
dius collective angle was zero. 

The aeroelastic simulation used in this study has two main limitations in its 
ability to model this rotor: 

1. The feathering and elastic axis are required to be coincident whereas for 
the OH6A blade the elastic axis is 1.2 inches (0.178xchord) ahead of the 
feathering axis. 

2. Non-coincident flap and lag hinges cannot be modeled whereas the lag 
hinge is 10.69 inches outboard of the flap hinge for the OH6A. 
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1 1.2 SIMPLIFIED OH6A BLADE 


In this study, a simplified model of the OH6A blade was developed. The 
objective here was to model the basic design characteristics of the blade as 
closely as possible within the constraints imposed by limitations of the analysis. 
In general, blade properties and their distributions were used as given. The 
major compromises made in the modeling are listed below: 

1. Since non-coincident elastic and feathering axes could not be modeled, 
these were assumed to be coincident at the elastic axis. 

2. Flap and lag hinges were assumed to be coincident at the flap hinge 
offset of 5.5 inches as given. Changes in the lag mode characteristics 
were minimized by tuning the torsional spring representing the strap 
pack cordwise flexibility so as to give the best possible first lag fre- 
quency. 

3. The torsional spring representing control system stiffness was tuned to 
give the best possible first torsional frequency. 

Other than these compromises, the data on the OH6A blade was used as 
given. Using this simplified model of the OH6A blade, natural rotating fre- 
quencies and mode shapes were calculated. Table 11.2 compares the first six 
rotating natural frequencies obtained from model of this study and the actual 
OH6A blade. 

The first mode in each of flap, lag, and torsion for the simplified model 
closely matched the frequencies and mode shapes provided for the OH6A. The 
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TABLE 11.2 


OH6A Rotating Natural Frequencies 


OH6A 

Simplified Model 

1 L, 0.442/rev. 

1 L, 0.441 /rev. 

1 F, 1 .028/rcv. 

IF, 1.028/rev. 

2F, 2.701 /rev. 

2F, 2.389/rev. 

IT, 3. 143/rev. 

IT, 3. 144/ rev. 

3F, 4.645/rev. 

2L, 4.442/rev. 

2L, 4.809/rcv. 

3F, 4.687/rcv. 


torsional components of the last three modes differed substantially from those 
given in the OH6A modal data. Furthermore, the last two modal natural fre- 
quencies appear in a reversed order. Within the constraints of the aeroclastic 
model used, it was felt that the representation of the modal characteristics 
provided above, was adequate. 


1 1.3 OPEN-LOOP HHC FOR THE OH6A 

In Ref. 67 open and closed-loop HHC flight tests were carried out at flight 
speeds of 60, 70, 80, and 100 knots. It was decided to do simulations at the 
highest of these speeds, where vibrations were most severe. The speed of 100 
knots is equivalent for the OH6A to an advance ratio of n =0.253. Plots of the 
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variation in vertical and lateral vibrations at the pilots scat, given in g's, with 
the phasing of a 1/3° amplitude lateral HHC input were given in Ref. 67. The 
plot of vertical and lateral vibrations for the 100 knot flight condition is re- 
produced in Fig. 1 1.2. Calculations were made for the model of this study to 
produce similar plots of the changes in vertical, lateral, and longitudinal hub 
shears with changes in the phase of collective, lateral, and longitudinal HHC 
inputs. These plots are presented in Figs. 11.3, 11.4, and 11.5 

These plots show that lateral and longitudinal shears respond fairly well to 
reduction through HHC inputs. Some phasing of HHC input in any one of the 
three control channels is capable of reducing either the lateral or longitudinal 
shear considerably below its baseline value with just this 1/3° of HHC ampli- 
tude. Vertical shear on the other hand, responds very well to lateral or longi- 
tudinal HHC but for collective HHC gives values well above the baseline for 
any phasing of a 1/3° input. Comparing the results for lateral HHC input 
from Ref. 67 and this study, the phasing for minimum vertical and lateral 
shears are coincident at 320° for Fig. 1 1.2, the plot from Ref. 67, while from 
the plot of this study, Fig. 1 1.4, the minimum for lateral shears is at 140° and 
that for vertical shears is at 85°. For both the lateral and the longitudinal 
control channels, there is an HHC phasing which reduces all three shears be- 
low the baseline for this 1/3° input. 

There is little correlation between the fuselage vibrations reported in Ref. 
67 and the hub loads determined in this study. For instance Ref. 67 reported 
the baseline vertical vibration level as being 2.5 times the baseline lateral vi- 
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bration level. Conversely the vertical hub shear of this study was found to be 
only 30% of the lateral or longitudinal hub shear. In fact this is not unrea- 
sonable. In Ref. 62 an analytical analysis of HHC was done on a model which 
provided both hub loads and loads at the pilot's scat. Exactly this same in- 
version was seen with lateral and longitudinal vibrations being predominant 
at the hub and vertical vibrations being predominant at the pilot's seat. In an 
actual helicopter, the vibrations at a point in the fuselage due to rotor hub 
loads will be changed in phase and magnitude due to the dynamic character- 
istics of the fuselage. In addition, vibrations due to the powerplant, the tail 
rotor, impingement of the rotor wake on the empennage, and other factors will 
contribute to fuselage vibrations. Therefore it remains to be shown what lo- 
cations for vibration sensors will be most effective in minimizing overall vi- 
brations. 


1 1.4 CLOSED-LOOP HHC FOR THE OH6A 

The HHC algorithm used for closed-loop flight tests in Ref. 67 was the 
cautious, global controller referred to as controller 4 in Ch. 6. This same con- 
troller was used to apply HHC to the simplified model of the OH6A rotor at 
a 100 knot or n = 0.253 flight condition. Five iterations of HHC were carried 
out. Figure 11.6 shows the iteration history of the HHC input angles com- 
manded by the controller. By the second iteration HHC angles are very close 
to their optimums. The collective channel was shown in the last section to be 
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quite ineffective in reducing vertical hub shears, and this is reflected in the two 
collective HHC inputs being essentially zero. 

Figure 11.7 shows the iteration history of the hub shears being minimized 
by HHC. The shears are essentially totally suppressed, being less than 1% of 
their baseline values, in all three directions. In Ref. 67 the vertical, lateral, and 
longitudinal pilot's seat vibrations were reduced to respectively 7%, 23%, and 
77% of their baseline values at 100 knots. Figure 11.8 shows the iteration 
history of the hub vibratory moments as HHC is applied. There is no increase 
in the yawing moment but substantial changes in the pitching and rolling mo- 
ments which increase respectively to 5.3 and 2.4 times their baseline values. 
Again the significance of any increase in moments as shears are reduced is 
dependent on actual fuselage dynamics, which are not included in this model. 


1 1.5 POWER REQUIREMENTS AND STABILITY 

When HHC was applied to the simplified model of the OH6A rotor, the 
required rotor power increased by 'A % . There was no indication of a decrease 
in required rotor power as was reported in Ref. 67. This small increase in re- 
quired power would present no problem in itself, but in any actual prediction 
of required power for HHC implementation, the power requirements of the 
servo-actuators used to produce the HHC pitch changes would also have to 
be taken into account. There were no appreciable negative effects on aeroe- 
lastic stability when HHC was applied to this rotor. The real part of the 
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characteristic exponent associated with the second lag mode, which was closest 
to its stability margin, increased in stability by 18% while the other modes had 
changes in stability of 4% or less. 
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Chapter XII 

CONCLUDING REMARKS 


This study has presented a rotor aeroelastic analysis capable of modeling 
the effects of HHC in forward flight. This is a major extension of an existing 
analysis which includes several features which make it a superior model for 
HHC studies. The original analysis incorporates structural nonlinearities in a 
consistent manner and provides direct stability information. An improved, 
implicitly derived, implementation of a finite-state, time-domain unsteady 
aerodynamic formulation has been incorporated in the model to capture high 
frequency aerodynamic effects due to HHC pitch changes. Trim values re- 
quired for the analysis have been calculated using an improved helicopter trim 
procedure which accounts for flap, lag, and torsional blade deformations. This 
complete analysis has then been coupled with a procedure for determining the 
rotor vibratory hub loads and a set of control routines used to calculate the 
HHC inputs needed to minimize these vibratory loads. 

Using this model, a study of the effects of unsteady, as opposed to quasis- 
teady, aerodynamics on a hingeless rotor was carried out. The unsteady 
aerodynamic formulation, which also included several apparent mass terms 
which had been neglected in the quasisteady formulation, was found to have 
moderate effects on the stability and low frequency response of the blade. In 
order to compare high frequency response, the 4/'rev. rotor hub shears were 
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compared for the two aerodynamic formulations. It was found that there can 
be large differences in high frequency response between the two aerodynamic 
formulations especially at higher advance ratios. This difference is noticeable 
particularly in the vertical hub shear response which is primarily affected by 
the changes in lift due to the unsteady aerodynamics. The differences were 
pronounced when the blade was excited at high frequencies in pitch, as is the 
case with HHC. These results indicate that high frequency unsteady aero- 
dynamic effects should be included when modeling HHC. A flap trim and an 
improved full trim procedure were compared. The principal differences in trim 
values were in the collective at low advance ratios and in the cosine cyclic and 
angle of attack at high advance ratios. The new trim values gave a consider- 
able difference in low frequency flap response and also led to moderate 
changes in vertical hub shears at advance ratios above ^ = 0.3. Blade stability 
was not appreciably affected by the new trim solution. The convergence cri- 
terion used for quasilinearization in this study was more precise than that of 
Ref. 8 and results were shown indicating that an improper convergence crite- 
rion could lead to poorly converged higher harmonics of the response solution. 
This led to erroneous values for the vertical hub shears which may have con- 
tributed to inconsistent optimization results in Ref. 8. 

A preliminary study of the effects of open-loop HHC was carried out on a 
hingeless rotor. HHC inputs of 1/3° in the collective, lateral, and longitudinal 
control channels were phased at 90° intervals between 0° and 360° and the 
corresponding changes in vibratory hub shears and moments were calculated. 
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With appropriate phasing, lateral and longitudinal shears could be driven be- 
low the baseline with control input in any channel using either aerodynamic 
model. Vertical shears could only be reduced below the baseline by longitudi- 
nal input for the unsteady aerodynamic model however they could be reduced 
below the baseline by airy channel of input with the quasistcady aerodynamic 
model. Similar phasing for minimums of lateral and longitudinal shears were 
found for the two aerodynamic models but the phases for vertical shear min- 
imization differed greatly. Calculations were made to show the variation in 
4/'rev. harmonics of hub shears as the HHC input angle magnitude was varied 
for a set phase angle. The relationships were found to be quite linear over the 
range 0° to 3° of amplitude for both the quasistcady and the unsteady models. 

In order to initiate closed-loop studies of HHC, initial estimates of HHC 
transfer matrices relating HHC inputs to changes in hub vibratory shears were 
calculated at an advance ratio of ^=0.3. This was done by individually in- 
crementing each HHC degree of freedom by 1/3° and measuring the resultant 
changes in hub shears. These transfer matrices were then used to implement 
fixed gain closed-loop control on both the quasisteady and unsteady aero- 
dynamic models. In both cases shears were very successfully reduced, indi- 
cating that the off-line estimates of the transfer matrices were quite accurate. 
With the unsteady model there was a tendency for the vertical shears to drift 
upward after minimization. The quasisteady model was next used to compare 
the performance of a local and three global HHC algorithms. Based on these 
results, it was decided to conduct the rest of the investigation using two con- 
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trollcrs. a local controller using Kalman filter identification of only the transfer 
matrix, and a global controller using Kalman identification of the transfer 
matrix and the baseline vibration level. 

Next cautious and deterministic versions of these two controllers were im- 
plemented for the quasisteady and unsteady aerodynamic models. With qua- 
sistcady aerodynamics there were very small differences in control sequence 
and shear minimizations. With unsteady aerodynamics the deterministic con- 
troller tended to slightly overshoot the optimum control values and the local 
controller was not as effective in holding the vertical shear to its minimum. 
This was attributed to the identification process of the transfer matrix for the 
local controller which was based on changes in the HHC inputs, rather than 
the total magnitude of the HHC inputs in the global case. Thus when the 
minimum was quickly reached changes in HHC input became essentially zero 
and identification to improve the transfer matrix estimate stopped. In order 
to gauge the ability of the controllers to respond to changing flight conditions, 
calculations were carried out wherein a step change in flight condition from 
fi = 0.3 to n =0.35 was applied from the optimal HHC solution at /x =0.3. The 
local controller responded with large oscillations in commanded control inputs 
and resultant shear minimizations however by the sixth iteration the hub 
shears were reduced considerably below the baseline for this advance ratio. 
Under the same conditions, the global controller responded quite smoothly and 
only small oscillations were observed. The reduced performance of the local 
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controller was again attributed to poor transfer matrix identification at the 
original baseline. 

In order to compare the differences in effects of HHC on articulated and 
hingeless rotors, HHC simulations were carried out on two roughly equivalent 
rotors, one with a hingeless and one with an articulated hub. It was found that 
in order to achieve a similar reduction in hub shears, the hingclcss rotor re- 
quired HHC angle magnitudes which were four times larger than those for the 
articulated rotor. This implies that considerably more power might be needed 
to drive the HHC servo-actuators for a hingeless rotor. While hub shears were 
being reduced, hub pitching and rolling moments for the hingclcss blade rose 
to four and six times their baseline values while for the articulated rotor they 
increase to only 1.3 and 0.8 times their baseline values. Moments for the ar- 
ticulated rotor remained low due to the flap and lag hinges which keep mo- 
ment transfer from the rotor to the non-rotating system to a minimum, and 
were an order of magnitude less than those for the hingeless blade. Subse- 
quently HHC was used to attempt to suppress simultaneously both the hub 
shears and the hub moments. For the hingeless blade this led to pitching and 
rolling moments which became half the baseline values, however this decrease 
was accomplished only at the expense of shear reductions of less than 30%. 
For the articulated blade, a small decrease in the moments, which were already 
small, was accompanied by a minor increase in shears. This shows that overall 
vibration reductions may be more difficult to achieve with HHC when applied 
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to hingclcss rotors because the simultaneous minimization of both hub shears 
and moments could be difficult to accomplish. 

Rotating blade root loads were also calculated for the hingeless and artic- 
ulated blades as HHC was applied. Both pcak-to-pcak and absolute blade 
root bending moments increased significantly with HHC application for the 
hingclcss blade, and feathering moment increased substantially for both 
blades. These increases need to be taken into accounted when designing the 
strength and fatigue properties of rotor blades on which one intends to use 
HHC. Application of HHC led to increases in required rotor power of less 
than 1% for both blades. Acroclastic stability margins for the two blades 
changed by less than 1% except for the stability of the second lag mode of the 
hingeless blade, for which the stability margin increased by 30%. 

A simulation of application of HHC to the rotor of an OH6A helicopter 
was conducted and the results were compared to actual flight test data. The 
model used was somewhat limited due to the inability to model the noncoinci- 
dent flap and lag hinges and noncoincident elastic and feathering blade axes 
of the OH6A rotor blade. Within these limitations, the model was designed 
so as to have similar geometric, structural, and mass properties to the actual 
rotor blade. Using this model, open-loop HHC simulations were carried out. 
Phasing or the HHC for shear minimization did not correlate well between the 
model and actual flight tests, but this may have been caused by the differences 
in the objective functions being minimized. In the flight tests vibratory accel- 
erations at the pilot's seat were minimized whereas in the simulation hub vi- 
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bratory shears were minimized. Due to fuselage dynamics and vibrations from 
sources other than the rotor, vibration minimization at a given point in the 
fuselage need not correspond to hub shear minimization, and in fact might lead 
to increases in some hub load components. Lateral and longitudinal HHC in- 
puts were found to be the most effective in reducing hub shears. When 
closed-loop HHC was applied to this model, all shears were suppressed to less 
than 1% of their baseline values at a speed of 100 knots. Concurrently, 
pitching and rolling hub vibratory moments increased to respectively 5.3 and 
2.4 times their baseline values. As HHC was applied, rotor required power 
increased by V 2 %. The small decrease in rotor power reported in Ref. 67 was 
not reproduced by the simulation. Acroclastic stability was not adversely af- 
fected by the application of HHC. Because of the lack of fuselage dynamics 
in the model, these results cannot be construed as any sort of prediction of 
actual OH6A response to HHC, but serve to show the ability of the analysis 
to model more realistic rotor configurations. 
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Figure 2.2: Geometry of the Undeformed Blade in the Rotating System. 
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Figure 2.3: Geometry of the Elastic Axis of the Deformed Blade. 
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Figure 2.4: Deformed Blade Cross Sectional Geometry. 
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Figure 3.1: Forces and moments on the helicopter in straight and level 

flight. 
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Figure 3.2: Flow chart showing procedure used to solve the full trim 

equations. 
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Figure 4.1: Geometry of motion of a typical airfoil. 
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Figure 5.1: Flow chart showing procedure for hub force and moment calcu- 

lation. 
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Figure 6.1: Schematic of helicopter plant model. 
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Figure 6.2: Block chart showing generation of the adaptive error vector for 

the extended LMS identification algorithm. 
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Figure 7.1: Flap and lag response, advance ratio /i = 0.4, results of Ref. 13, 

torsionally stiff blade, and baseline soft-in-plane blade. 
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Figure 7.2: 
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Figure 7.3: Flap, lag, and torsional response, soft-in-plane blade, advance 

ratio n = 0.4, previous and current unsteady aerodynamics. 
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Figure 7.4: Stability eigenvalues for flap and lag, soft-in-plane blade, previ- 

ous and current unsteady aerodynamics. 
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Figure 7.5: Flap, lag, and torsional response, soft-in-plane blade, advance 

ratio n=0.2 , quasisteady and unsteady aerodynamics. 
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Figure 7.6: Flap, lag, and torsional response, soft-in-plane blade, advance 

ratio ^ = 0.4, quasisteady and unsteady aerodynamics. 
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Figure 7.7: Flap, lag, and torsional response, stiff-in-plane blade, advance 

ratio p. =0.2, quasisteady and unsteady aerodynamics. 


187 





BLADE 


Figure 7.9: Flap response, soft-ir 

teady, quasisteady w 
mics. 


quasisteady 


“ quasisteady with apparent mass 


“ unsteady 




quasisteady 



Figure 7.10: Lag and torsional response, soft-in-plane blade, advance ratio 

//=0.4, quasisteady, quasisteady with apparent mass, and un- 
steady aerodynamics. 
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Figure 7.1 1: Stability eigenvalues for the first lag, first flap, and third flap 

modes, soft-in-plane blade, quasistcady and unsteady aero- 
dynamics. 
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Figure 7.12: Stability eigenvalues for the second lag, second flap, and first 

torsion modes, soft-in-plane blade, quasisteady and unsteady 
aerodynamics. 
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Figure 7.13: Stability eigenvalues for the second lag, first flap, and second 

flap modes, stiff-in-plane blade, quasisteady and unsteady 
aerodynamics. 
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Figure 7.14: Stability eigenvalues for the first lag, third flap, and first torsion 

modes, stiff-in-plane blade, quasisteady and unsteady aero- 
dynamics. 




Figure 7.15: Variation of 4/rev. vertical and lateral hub shears with advance 

ratio, soft-in-plane blade, quasisteady aerodynamics, calculated 
from converged and unconverged response solutions. 
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Figure 7.16: Variation of 4/rev. vertical shear and yawing moment with ad 

vance ratio, soft-in-plane blade, quasisteady and unsteady 
aerodynamics. 



Figure 7.17: Variation of 4/rev. lateral shear and pitching moment with ad- 

vance ratio, soft-in-plane blade, quasisteady and unsteady 
aerodynamics. 
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Figure 7.18: Variation of 4/rev. longitudinal shear and rolling moment with 

advance ratio, soft-in-plane blade, quasistcady and unsteady 
aerodynamics. 



Figure 7.19: Variation of 4 rev. vertical shear and yawing moment with ad- 

vance ratio, stiff-in-plane blade, quasisteady and unsteady 
aerodynamics. 
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Figure 7.20: Variation of 4/rev. lateral shear and pitching moment with ad- 

vance ratio, stiff-in-plane blade, quasisteady and unsteady 
aerodynamics. 
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Figure 7.21: Variation of 4/rev. longitudinal shear and rolling moment with 

advance ratio, stiff-in-plane blade, quasisteady and unsteady 
aerodynamics. 
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Figure 7.24: Flap, lag, and torsional response, soft-in-plane blade, advance 

ratio pL- 0.3, unsteady aerodynamics, flap trim and full trim. 
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Figure 7.25: Variation of 4/rev. vertical shear and yawing moment with ad 

vance ratio, soft-in-plane blade, unsteady aerodynamics, flap 
trim and full trim. 
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Figure 7.26: Variation of 4/rev. lateral shear and pitching moment with ad 

vance ratio, soft-in-plane blade, unsteady aerodynamics, flap 
trim and full trim. 





Figure 7.27: Variation of 4/rev. longitudinal shear and rolling moment with 

advance ratio, soft-in-plane blade, unsteady aerodynamics, flap 
trim and full trim. 
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Figure 8.1: Variation of 4/rev. vertical shear and yawing moment with HHC 

phase, collective input, soft-in-plane blade, quasisteady and un- 
steady aerodynamics, /i =0.3. 
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Figure 8.2: Variation of 4/rev. lateral shear and pitching moment with HHC 

phase, collective input, soft-in-plane blade, quasisteady and un- 
steady aerodynamics, // =0.3. 
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Figure 8.3: Variation of 4/ rev. longitudinal shear and rolling moment with 

HHC phase, collective input, soft-in-plane blade, quasisteady 
and unsteady aerodynamics, =0.3. 
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Figure 8.4: Variation of 4/'rev. vertical shear and yawing moment with HHC 

phase, longitudinal input, soft-in-plane blade, quasistcady and 
unsteady aerodynamics, n =0.3. 
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Figure 8.5: Variation of 4/rev. lateral shear and pitching moment with HHC 

phase, longitudinal input, soft-in-plane blade, quasisteady and 
unsteady aerodynamics, n =0.3. 
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Figure 8.6: Variation of 4/rev. longitudinal shear and rolling moment with 

HHC phase, longitudinal input, soft-in-plane blade, quasisteady 
and unsteady aerodynamics, /i = 0.3. 
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Figure 8.7: Variation of 4/rev. vertical shear and yawing moment with HHC 

phase, lateral input, soft-in-plane blade, quasisteady and un- 
steady aerodynamics, n =0.3. 




Figure 8.8: Variation of 4/rev. lateral shear and pitching moment with HHC 

phase, lateral input, soft-in-plane blade, quasisteady and un- 
steady aerodynamics, u = 0.3. 




Figure 8.9: Variation of 4/rev. longitudinal shear and rolling moment with 

HHC phase, lateral input, soft-in-plane blade, quasisteady and 
unsteady aerodynamics, n =0.3. 
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Figure 8.10: Variation of hub shear harmonics with HHC amplitude, qua- 

sisteady and unsteady aerodynamics, longitudinal cosine HHC 
applied, /i = .0.3. 
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Figure 8.1 1: Variation of hub shear harmonics with HHC amplitude, qua- 

sisteady and unsteady aerodynamics, longitudinal cosine HHC 
applied, u =0.3. 
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Figure 8.12: Variation of peak-to-peak hub shears with HHC amplitude, 

quasisteady and unsteady aerodynamics, longitudinal cosine 
HHC applied, n = 0.3. 
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Figure 8.13: Proportional change in hub shears with different 1/3° HHC in- 

crementations, quasisteady and unsteady aerodynamics, soft- 
in-plane blade, advance ratio fi= 0.3. 
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Figure 9.2: Iteration history of vertical, lateral, and longitudinal shears for 

fixed gain control, quasisteady aerodynamics, u =0.3. 
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Figure 9.3: Iteration history of HHC input angles for Fixed gain control, 

unsteady aerodynamics, n =0.3. 
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Figure 9.4: Iteration history of vertical, lateral, and longitudinal shears for 

fixed gain control, unsteady aerodynamics, /i = 0.3. 
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Figure 9.5: Iteration history of vertical shears, four different deterministic 

controllers, quasisteady aerodynamics, u =0.3. 



Figure 9.6: Iteration history of lateral shears, four different deterministic 

controllers, quasisteady aerodynamics, u =0.3. 
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Figure 9.7: Iteration history of longitudinal shears, four different determin- 

istic controllers, quasisteady aerodynamics, n = 0.3. 
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Figure 9.8: Iteration history of HHC input angles for deterministic versus 

cautious local control, quasisteady aerodynamics, u =0.3. 
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Figure 9.9: Iteration history of vertical and lateral hub shears for dctcrmin 

istic versus cautious local control, quasisteady aerodynamics, 
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Figure 9.10: Iteration history of longitudinal hub shear for deterministic 
VerS Q j Caut ' ous * oca * contro, » quasisteady aerodynamics, 
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Figure 9.1 1: Iteration history of HHC input angles for deterministic versus 

cautious global control, quasisteady aerodynamics, /i =0.3. 
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Figure 9.12: Iteration history of vertical and lateral hub shears for deter 

ministic versus cautious global control, ouasisteady aerodvn 
mics, u=0.3. 
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Figure 9.13: Iteration history of longitudinal hub shear for deterministic 

versus cautious global control, quasisteady aerodynamics, 

= 0 . 
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Figure 9.14: Iteration history of HHC input angles for deterministic versus 

cautious local control, unsteady aerodynamics, n =0.3. 
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Figure 9.15: Iteration history of vertical and lateral hub shears for deter- 

ministic versus cautious local control, unsteady aerodynamics, 
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Figure 9.16: Iteration history of longitudinal hub shear for deterministic 

versus cautious local control, unsteady aerodynamics, tt =0.3. 



Figure 9.17: Iteration history of HHC input angles for deterministic versus 

cautious global control, unsteady aerodynamics, \i = 0.3. 
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Figure 9.18: Iteration history of vertical and lateral hub shears for deter- 

~ versus cautious global control, unsteady aerodynam 
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Figure 9.19: Iteration history of longitudinal hub shears for deterministic 

versus cautious global control, unsteady aerodynamics, n =0.3. 
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Figure 9.20: Optimal shears for deterministic and cautious versions of local 
and global controllers using quasisteady aerodynamics, n =0.3. 
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Figure 9.21: Optimal shears for deterministic and cautious versions of local 

and global controllers using unsteady aerodynamics, =0.3. 
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Figure 9.22: Comparison of optimal HHC inputs for the cautious version of 
the local controller using quasisteady versus unsteady aero- 
dynamics, pi =0.3. 
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Flap, lag, and torsional elastic tip response, baseline response 
and optimal response due to cautious global HHC, unsteady 
aerodynamics, soft-in-plane blade, n =0.3. 
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Figure 9.24: Total torsional geometric and elastic tip displacement, baseline 
response and optimal response due to cautious global HHC 
unsteady aerodynamics, soft-in-plane blade, // =0.3 
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Figure 9.25: Iteration history of HHC input angles for deterministic versus 

cautious local control, unsteady aerodynamics, step change 
from n =0.3 to n = .35. 
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Figure 9.26: Iteration history of hub shears for deterministic versus cautious 

local control, unsteady aerodynamics, step change from u =0.3 
to u= 0.35. 
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Figure 9.27: Iteration history of HHC input angles for deterministic versus 

cautious global control, unsteady aerodynamics, step change 
from n = 0.3 to p = .35. 
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Figure 9.28: Iteration history of hub shears for deterministic versus cautious 

global control, unsteady aerodynamics, step change from 
H =0.3 to u=0.35. 
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Figure 9.29: Baseline shears for ^ = 0.35, and shears five iterations after a 

step change from =0.3 to fi= 0.35, local and global control 
lers, unsteady aerodynamics. 
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Figure 10.2: Iteration history of vertical, lateral and longitudinal peak-to- 

peak hub shears, hingeless blade, 5% blade root offset, cautious 
global control, = 0.3. 
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Figure 10.4: Iteration history of HHC input angles for articulated blade, 5% 

blade root offset, cautious global control, fi = 0.3. 
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Figure 10.5: Iteration history of vertical, lateral and longitudinal peak-to- 

peak hub shears, articulated blade, 5% blade root offset, cau- 
tious global control, pi =0.3. 
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Figure 10.6: Iteration history of pitching, rolling and yawing peak-to-peak 

hub moments, articulated blade, 5% blade root offset, cautious 
global control, n = 0.3. 
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Figure 10.7: Iteration history of HHC input angles for hingeless blade, 5% 

blade root offset, 6x12 cautious global control, n =0.3. 
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Figure 10.8: Iteration history of vertical, lateral and longitudinal peak-to- 

peak hub shears, hingeless blade, 5% blade root offset, 6x 1 2 
cautious global control, fx =0.3. 



ZOZQ-5W2I/1 — OZ<j > < j 3 



Figure 10.9: Iteration history of pitching, rolling and yawing peak-to-peak 

hub moments, hingeless blade, 5% blade root offset, 6x12 cau 
tious global control, u = 0.3. 


Hub Load Component 


Figure 10.10: Comparison of peak-to-peak hub loads for baseline, 6x6 opti- 

mum, and 6x12 optimum, hingeless blade, 5% blade root off 
set, cautious global control, // = 0.3. 









Figure 10.1 1: Iteration history of HHC input angles for articulated blade, 

5% blade root offset, 6x12 cautious global control, fx =0.3. 
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Figure 10.12: Iteration history of vertical, lateral and longitudinal peak-to 

peak hub shears, articulated blade, 5% blade root offset, 
6x12 cautious global control, n =0.3. 
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Figure 10.13: Iteration history of pitching, rolling and yawing peak-to-peak 
hub moments, articulated blade, 5% blade root offset, 6x12 
cautious global control, u = 0.3. 



Hub Load Component 


Figure 10.14: Comparison of peak-to-peak hub loads for baseline, 6x6 opti 
mum, and 6x12 optimum, articulated blade, 5% blade root 
offset, cautious global control, n=0.3. 
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Figure 10.15: Iteration history of axial, cordwise, and vertical peak-to-peak 
rotating blade root shears, hingeless blade, 5% blade root off- 
set, cautious global control, u =0.3. 
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Figure 10.16: 


Iteration history of feathering, flapping, and lagging peak- 
peak rotating blade root moments, hingeless blade, 5% bU 
root offset, cautious global control, u =0.3. 
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Figure 10.17: Iteration history of maximum values of feathering, flapping, 

and lagging rotating blade root moments, hingeless blade, 5% 
blade root offset, cautious global control, u = 0.3. 
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Figure 11.1: OH6A main rotor hub. 
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Figure 11.2: 


Variation of 4/rev. vertical and lateral pilot seat vibrations witl 
Phaseo 1/3° lateral HHC input, flight test results of Ref. 67 fo 
the OH6A helicopter, 100 knots. 
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Figure 1 1.3: Variation of 4/rev. vertical, lateral, and longitudinal hub shears 

with HHC phase, 1/3° collective HHC input, OH6A rotor, un- 
steady aerodynamics, /r =0.253. 



Figure 1 1.4: Variation of 4/rev. vertical, lateral, and longitudinal hub shears 
with HHC phase, 1/3 lateral HHC input, OH6A rotor, un- 
steady aerodynamics, n =0.253. 
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Figure 1 1.5: Variation of 4/rev. vertical, lateral, and longitudinal hub shears 

with HHC phase, 1/3° longitudinal HHC input, OH6A rotor, 
unsteady aerodynamics, n= 0.253. 
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Figure 1 1.7: Iteration history of vertical, lateral, and longitudinal 4/rev. hub 

shears for OH6A rotor, cautious global controller, \i = 0.253, 
unsteady aerodynamics. 
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Figure 1 1.8: Iteration history of pitching, rolling, and yawing 4/rev. hub 

moments for OH6A rotor, cautious global controller, n =0.253, 
unsteady aerodynamics. 
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